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Preface 


I am highly honored to welcome you all to the 7th GNSS Workshop — International Symposium on 
GPS/GNSS on behalf of Korea GPS Technology Council. Thank you for joining us. Your participation 


is vital to the conference, and your satisfaction is our goal and the best measure of our success. 


As many of you may agree, GPS/GNSS is a revolution of navigation which gives us extremely 
accurate position and timing information than ever obtained. Together with the mobility technology, 


GPS/GNSS is now creating new challenging technology areas and market, namely L-commerce. 


GNSS Workshop — International Symposium on GPS/GNSS is the only navigation conference in Asia 
sponsored by the ION, the world’s leading professional society. The purpose of this International 
Symposium is to provide the current status and to make clear what should be done next as Steps 
Forward concerning GPS-related technology through the international cooperative efforts of exchanging 


knowledge and experience. 


There are many people to thank for this symposium. We have had strong support from Ministry of 
Information and Communication, which sponsored this symposium. The Organizing Committee has 
provided excellent leadership and an enormous amount of work. They have organized an outstanding 
technical program, and the symposium would not be possible without the efforts of all the authors. 
Those who present the tutorials and those who staff the exhibits also deserve great thanks. I am 
deeply grateful to all the distinguished speakers in the plenary. Finally, I want to acknowledge the 
support and participation of all the organizations including the commercial ones which permitted and 
encouraged their employees to spend the time, efforts and resources needed to make this symposium 


most educational, idea stimulating and contact producing. 


Sang Jeong Lee 
General Chairman, the 7th GNSS Workshop — International Symposium on GPS/GNSS 
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Abstract: This paper presents the technologies to make up a GPS/INS integrated system. Based on the 


fundamental backgrounds of INS and GPS systems, such as the configurations, fundamental principles and 


error characteristics, it describes how to configure a GPS/INS integrated system. And the performance of the 


integrated system is analyzed by the simulation. In addition the technology trends and processes of a 


GPS/INS integrated navigation system are also presented. 


Keywords: INS, GPS, Integrated Navigation 
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Dead Reckoning Navigation System for Autonomous 


Mobile Robot Trajectory Following 
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Abstract: Dead reckoning navigation system is designed to make the autonomous mobile robot 
follow the desired trajectory. Sensor-based autonomous navigation architecture is composed of 
gyroscope, and odometer using digital rotary encoder. Steering controller and the speed controller are 
used to change the direction and the speed of the mobile robot. One embeded computer is used as a 
main processor, and radio frequency modem is used to transmit the command as well as to monitor 
the status of the mobile robot. Two-dimensional planar dynamics is considered in this study. The 
mobile robot is considered as a rigid body, and the dynamics of the actuator is approximated as first 
delay filter. System identification experiments are performed to obtain the dynamic model of the 
mobile robot system, and the identified model is analyzed through time domain analysis with 
feed-forward control law. 


Keywords: Dead reckoning, autonomous mobile robot, system identification 
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ABSTRACT 

This paper investigates observability properties of 
strapdown inertial navigation system (INS) aided by the 
Global Positioning System (GPS). Observability 
analysis is made on a linear error dynamics model of 
navigation equations in the Earth-centered Earth-fixed 
(ECEF) frame for both loosely-coupled and tightly- 
coupled integrations. Position and velocity 
measurements are given from multiple GPS antennas 
mounted on the vehicle. The motivation to consider 
multiple GPS antennas is that the lever-arms of them 
play an important role in observing vehicle attitude and 
biases of inertial motion unit (IMU) during navigation. 
Biases of IMU are assumed to be constant. Necessary 
and sufficient conditions for the instantaneous 
observability on the GPS/INS integrations are given. 

A numerical simulation is given to show the 
effectiveness of multiple GPS antennas in estimating 
vehicle attitude and biases of IMU when IMU has 
considerable magnitude of biases. 


INTRODUCTION 

Improved navigation can be realized by the integration 
of GPS and INS. GPS receivers provide position and 
velocity of vehicle with bounded accuracy. The accuracy 
is independent of elapsed time from the start of 
measurement. However, GPS receivers can be 
considered as discrete-time sensors; in many cases, their 
sampling period is about one second. Occasionally, 
measurement is not available during loss of lock on 
satellites due to shading of GPS antennas or radio- 
frequency (RF) interference. The INS is a continuous- 
time measurement system. It offers short-term stability, 
but has poor long-term stability because of bias of IMU 
which consists of gyros and accelerometers. Using the 
above complementary properties of GPS and INS, 
various integration schemes are under development to 
overcome shortcomings of each sensor system. 

An attractive scheme of the integration is to estimate 
biases of IMU in strapdown INS during navigation with 
GPS measurement. With this integration scheme, an 
accurate, low-cost navigation system that provides long- 
term stability and continuous measurement can be 
constructed with relatively low-accurate inertial sensors. 


There have been several approaches to the estimation of 
biases of IMU using GPS measurements. A technique to 
determine attitude of a vehicle by GPS antennas without 
IMU was developed in [1]. The attitude measured by 
GPS in [1] was employed in [2] and [3] to estimate gyro 
bias. A nonlinear observer for attitude and IMU bias 
was suggested in [4]. The observer was proven to be 
exponentially convergent for the IMU biases that were 
modeled as Markov processes. A_ sensor fusion 
technique is introduced in [5] to estimate attitude and 
IMU biases with GPS, IMU, and air-data sensors. 

In this paper, observability properties of strapdown INS 
aided by GPS are investigated. The biases of IMU are 
assumed to be constant. The time-constant model could 
be useful for biases that change very slowly compared 
with vehicle Velocity and position 
measurements are given from multiple GPS antennas on 
the vehicle. Observability analysis is made on a linear 
error dynamics model in the Earth-centered Earth-fixed 
(ECEF) frame. The state in the error model consists of 
3-dimensional biases of gyro and accelerometer, and 
errors for position, velocity, and attitude. Analysis is 
made on both loosely-coupled and tightly-coupled 
GPS/INS integrations. 


dynamics. 


NOMENCLATURE 


Throughout this paper, the following notations are given: 
o;,: column vector of angular velocity of frame b 


relative to frame a, decomposed in frame c. 
: position vector decomposed in frame a. _ 
: velocity vector decomposed in frame a. 


: rotation matrix from frame a to frame b. 


Qi, : Skew-symmetric cross product matrix of @‘,. 
O : estimated value of (). 

6(): estimation error of (). 

6) : time derivative of (). 


(|: absolute value of (). 


( )x(_ ): cross product of two vectors 


The navigation frames used in the paper are : 
i-frame : Earth-centered inertial (ECI) frame. 
e-frame : ECEF frame. 
b-frame : Body-frame. 


NAVIGATION ERROR PROPAGATION MODEL 
The navigation equations in the ECEF frame are [6] 


Peay (1.1) 
Vo = Rf? -200 xV +g" (1.2) 
Ry = R,Qh, (1.3) 


where f’ is specific force in body frame and g“ is 


gravity in the ECEF frame. The corresponding INS 
mechanization differential equations are 


pave (1.4) 
Vi = Rf? Qo x Pee ge (1.5) 
Re = RO, (1.6) 
= OF - Rot (1.7) 
where f/f” and @% are measurements from 


accelerometer and gyro, respectively. Let the 
mechanization errors are modeled as 


P' = P+ 5P (1.8) 
Vio=V'+6V (1.9) 
Re = RU + [Pd (1.10) 
fraf? ret, (1.11) 
Oj, = Wy, +E, + W, (1.12) 


where y+ is attitude error, [yx ] is the cross product matrix 


of y . «¢, is accelerometer bias vector, @, is 


i a 


accelerometer noise, ¢€, is gyro bias vector, and 0, 1s 


gyro noise. Bias vectors é, and €, are assumed to be 
constant. Then the linearized error propagation equations 
are 

dP = OV (1.13) 


dV = GSP-Q.5V - RFy + Re, (1.14) 
y =-Oy te, +1, (1.15) 
é,=0 (1.16) 
o=0 (1.17) 
where 6-= ,F is the cross product matrix of /’, 


Q,=2Q°, R and Q are the simplified notations of 


t 


R’ and Q"., respectively. 


GPS MEASUREMENT ERROR MODEL 
LOOSELY-COUPLED GPS/INS SYSTEMS 


Measurements from GPS receivers are modeled as 


Pra P'+Ri+y, (2.1) 


Vi =V' + ROL +y,, fel 2yom (2.2) 


y? 


where P° and VY; are position and _ velocity 


measurements from the jth GPS receiver antenna, 
respectively, J; , which is called level arm of jth GPS 
antenna, is the position of jth GPS receiver antenna 
relative to that of IMU decomposed in the body frame, 


Q, is the simplified notation of 9’,, v,, andy, are 


position and velocity measurement error of the jth GPS 
receiver antenna, respectively, and m is the number of 
GPS receiver antennas. Estimation for measurements is 
given as 


(2.3) 


Vom ROR, Gd (2.4) 
The estimation errors for measurements are defined as 


Pe = P+ SP (2.5) 


! 


VeaV! +6V.. (2.6) 


Then, the linearized measurement estimation errors can 
be shown as 


OP’ =5P-RLy-v, (2.7) 


6V! = BV + RL Q-QL, Wy -Le,} — RL, wy ~ v, (2.8) 


where L; is the cross product matrix of /,. 


TIGHTLY-COUPLED GPS/INS SYSTEMS 
Measurements from GPS receivers are modeled as 


ya (2) 
p= 


el) e 
P -P tet, +77; 


(2.9) 


yi) 
P, = 


(2.10) 


ye ; Ope 
ye —Vil+ ci, +60, i=1,2,3,4 
where. pis the pseudorange of the ith satellite from 


the jth GPS receiver antenna, P“” is position of the ith 
satellite decomposed in the ECEF frame, c is the speed 


of light, 1, is clock bias of the jth GPS receiver, a 


J 


is the composite of errors produced by atmospheric 


delays, satellite ephemeris mismodeling, receiver 


s() 


tracking error, etc., \” is the pseudorange rate of the ith 


satellite from jth GPS receiver antenna, V“” is velocity 


of the ith satellite in the ECEF frame, ¢ is clock drift of 


the jth GPS receiver, ¢‘ is the measurement error. 


Estimation for measurements is given as 


pd =|PO — Pela ch (2.11) 
Ai = pe Vil +c, (2.12) 
The estimation errors for measurements are defined as 
BY = p+ 6p (2.13) 
p= p® + 6p (2.14) 
f=t,+6t, (2.15) 
i, =1,+61, (2.16) 
Let 
so) faa] fa] fee 
dp dp 7 Gt 
ry ’ |, OL S| a dyes lee AQT 
Pp; 6p" Pp, op’? i hi G: he ( ) 
op” | op” a o° 


Then, the measurement estimation errors can be shown 


as 


(2.18) 


(2.20) 


ov’ 
6p,=H,| / |= 
eA , 
where 
[ api” 
oP; 
ap 
i 
oP? 
A, = (3) 
ap, 
oP’ 
dp 
j 
| 
OBSERVABILITY PROPERTIES OF GPS/INS 
LOOSELY-COUPLED GPS/INS SYSTEMS 
Let 
x=[6P 8V’ y et a] (3.1) 
yu= [BR -- OR ays ors] (3.2) 
c,=[/ 0 ~ RL, 0 0] (3.3) 
C,=[0 7 -R(Q,L,-£,0) -RL, 0] (3.4) 


where ( ) denotes transpose of a vector or a matrix. 


Then, the equations of errors for INS mechanization and 
measurement estimation are 


xX=Ax+w 


(3.5) 


=C x+v 


Ym aN (3.6) 


where y,, is estimation error for measurement from m 
GPS receiver antennas, wand v are the first order 
approximation errors, and 


i orl 
OF. <0 00 i 
G -Q, -RF 0 R a 
A=/0 0 -Q FF OLE,=] (3.7) 
0 0 0 0 0 
0 0 0 00 
eee 
Let 
Pid. 020-40 
07000 
T=|0 0 7 0 0 (3.8) 
0.0 Qi d 
O00 F OF 


iL. 


(2.19) 


I 0 RL, 0 0 
0 ] 0 RL, 0 
Le 6 0 I 0 0} (3.9) 
0 0 Q i 0 
-R'G R'Q. F-R'GRL, R'Q.RL,-LQ 1 
0 | 0| 0 | 0 
0 0 0 0 
eet Lael tees bh lex =|. 0 (3.10) 
0 0 0 1-1, 
| 0 0] 0 0 
Xp Sit tye er hi = ben te bee CAD 


Then, we have the following theorem: 
Theorem I: Suppose 1; , 1,, and J; 


are linearly 
independent. Then: 


i) (C,,,4) is observable for m>3 


ii) If ’, and (/,x/,) are linearly independent 


andw’, #0, then (C,,A4) is observable. Otherwise, 


x,, is an unobservable mode of (C,, A). 
iii) The following three conditions, 


(1) o!, is parallel with J,x/,or w/,=0, 
(2) w) is parallel with /,-/, or ow) =0, 
(3) GRI,xl,=RF(l,-|), 
are the only conditions for (C,,A) to have two 


unobservable modes. In this case,x,, and x,, are 


the unobservable modes. 


iv) x,, is an unobservable mode of (C,, 4) ). 


v) If /,=0, then (C,,A) has three unobservable 


modes, x,,. 
Proof: See Appendix. 
TIGHTLY-COUPLED GPS/INS SYSTEMS 
For simplicity of expression, it is assumed that the 


multiantenna measurement system is implemented by 
means of a single multiantenna receiver, so that 


Lai 


_ 14- 


Let 


x,=[x" cét cdi] (3.12) 
Yom = [80 6a Ba BD] (3.13) 
re ee CH 0. [04 3.14 
te, 0 0 ’ a ’ eo =[ ] ( 7 ) 
Cc, oO Ce 20 
Cm] zl ma “| (3.15) 
Then, 
X,= A,X, +6, (3.16) 
Vom = ConXp (3.17) 
where 
A,= Ae 3.18 
BO 0 A (3. ) 
Ors (Ca Co. Cr Cia : (3.19) 
Let 


Xx x, x2 Xoo 
X pyo -| ee -| Al X p72 -| 0 | Xp? -| 0 } (3.20) 


Then, the following theorem on the observability of 
tightly coupled GPS/INS system can be inferred from 
the above Theorem 1. 

Theorem 2: Suppose /,, /,and /; are linearly independent 
and H; is of full rank. Then, the observability conditions 


for (C 


Zink A) with Kins Ries 


A,) are the same as (C 


me? 


x and x,., instead of x, 4%, 4%, and x 


py2? yl? yo? x2 


respectively. 

Remark J: Let the level arms /, , />, and /; are linearly 
independent. Then, it can be shown that GPS/INS 
system is instantaneously observable with only position 
measurement from the three GPS antennas for both 
loosely-coupled and tightly-coupled systems. 


SIMULATION RESULTS 

A numerical example is given to demonstrate the 
behavior of multiantenna GPS measurement system in 
estimating attitude of vehicle and biases of IMU. 
Responses of multi-rate extended Kalman filter is 
simulated on a simple vehicle trajectory. Estimation of 
navigation state is updated with the IMU measurement 
at 10Hz. The error of the estimation is corrected with 
position measurement from four GPS antennas at 1Hz. It 


is assumed that positions of GPS antennas are obtained 
from the double differenced carrier phase measurement. 
Integer ambiguity problem is assumed solved so that 
GPS position measurement has cm level accuracy. 
Length of each level arm of GPS antennas is about 1.7m. 
All the noises in the GPS receiver and IMU are assumed 
to be Gaussian white. Standard deviation of each 
component of GPS position measurement noise is set to 
5cm. Bias and standard deviation of each component of 
accelerometer noise are approximately [0.06 -0.03 0.02] 
and 0.01 in m/s’. Bias and standard deviation of each 
component of Gyro noise are [-0.06 0.1 -0.15] and 0.036 
in degree/s. Initial estimation errors for roll, pitch, and 
yaw in degrees are 0.5, -0.5, and 1.0, respectively. 

The path and attitude of vehicle in the simulation are 
given in Fig. 4.1 and Fig 4.2. Estimation errors for 
attitude and biases are given in Fig 4.3, Fig 4.4, and Fig 
4.5. These figures show that the GPS/INS system 
estimates attitude of vehicle and biases of IMU quite 
accurately. 
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Fig 4.2: Vehicle attitude 
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Fig 4.3: Attitude estimation error 
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Fig 4.4: Accelerometer bias estimation error 
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Fig 4.5: Gyro bias estimation error 
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CONCLUSIONS 

In this paper conditions for the instantaneous 
observability of linear error model of navigation 
equation are given. It is shown that the error model is 
instantaneously observable for all conditions of motion 
with measurement from at least three GPS antennas. 


With three or more GPS antennas, velocity measurement 
is not required and only position measurement is enough 
for the system to have instantaneous observability. There 


is at least one unobservable mode with one GPS antenna. 


If the location of the GPS antenna coincides with inertial 
sensor position, there exist at least three unobservable 
modes. The error model is instantaneously observable 
with measurement from two GPS antennas unless the 
vehicle angular velocity is zero or parallel with the cross 
product of the two position vectors of the GPS antennas. 
A numerical simulation shows that low-cost IMU with 
carrier phase differential GPS system that has cm-level 
accuracy can be an accurate and reliable navigation 
sensor system. 
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APPENDIX 
Proof of Theorem 1: Introduce a new state x by a 
linear transformation such that 


x=T)x. 


Then, we have 


Let 

Gy = CA =[(7), (G2), (Ce), (Ee), (E),] 
with k=0,L ,n—1 where n (=15) is the size of x. Then, 
the proof of i) can be obtained from the relation 
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The proof of ii) and iii) follows from the relation 
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- LQ? -F)}, k=3,L ,n-1 
where 


Cy, = GRL, - RF + R(Q,L, - LQ)Q+ Q,RLQ 
Cys = GRL, - RF + R(Q,L, -— L,Q)Q+ O,RLQ. 


iv) can be proved with C* which can be obtained in the 
same way as before for k=0,L ,n—-1. The proof of v) 
can be obtained using 7, in stead of 7, in the above 


state transformation equation. W 
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Abstract: In this paper, the error compensation method of the low-cost IMU is proposed. In general, 


the position and attitude error calculated by accelerometers and gyros grows with time. Therefore 


the additional information is required to compensate the drift. The attitude angles can be bound by 


accelerometer mixing algorithm and the heading angle can be aided by single antenna GPS velocity 


information. The Kalman filter is used for error compensation. The result is verified by comparing 
with the attitude calculated by medium-grade IMU, LP-81. 
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Abstract: This paper designs an integrated system using a GPS attitude determination system and a low cost IMU. The GPS attitude 
determination system calculates the precise attitude of vehicle using carricr phase measurements. The integration filter should 
compensate gyro error and attitude errors more precisely. In addition. the misalignment error of the GPS antennas should also be 
compensated. This paper also presents an algorithm for compensating the alignment error. 


Keywords: INS, GPS, attitude. 
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Abstract: This paper describes a carrier tracking loop design for high maneuver applications. A 


third order PLL(Phased Lock Loop) is designed for accurate carrier phase measurement, and 
IMU(Inertial Measurement Unit) information is used as external aiding for robustness against high 
dynamic stress. From various test results, we confirmed that adequate loop filter bandwidth and 
PIT(Predetection Integration Time) design is needed according to application environment, and 
external aiding could be useful information for carrier tracking loop design, especially for high 


maneuver application 
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Abstract: In this paper, a software correlator GPS algorithm under asynchronous network channel is proposed. The signal processor takes 
snapshot of the sampled IF signal from the GPS RF section. The processing for code and carrier acquisition and calculating the navigation solution 
are implemented using the digital signal processing techniques. In order to achieve fast code acquisition, correlation of the incoming GPS signal is 
performed using the DFT method. After code acquisition, to reduce the Doppler shift effect on the code and increase the measurement accuracy of 
the code correlation, the interpolation method or EML discriminator compensation is used. The performance of the proposed processing algorithm 
is first evaluated using Matlab/Simulink. In order to analyze the performance of the designed algorithm, using the signal acquisition board that 
samples and saves GPS IF signal form the Mitel GPS RF chip set, some experiments are performed and the results are analyzed. 
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Abstract: This paper designs a GPS attitude finding system. In order to determine the integer ambiguity, ARCE 
(Ambiguity Resolution with Constraint Equation) is adopted. The test results show that more reliable navigation 
system can be possible since the system provides all navigational information such as position, velocity, time and 


attitude. 
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Fig. 10. Digital processing board. 
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Fig. 11. /Q accumulation in noisy environments. 
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Abstract 


: It is well known that ambiguity resolution is essential for precise carrier phase differential GPS. However, in practice, 


there are still many problems in resolving the ambiguity in kinematic mode, especially in the urban area. The multipath in reccived 
signal, the frequent change in visible satellites, and the cyclic slips make the ambiguity resolution very difficult task in real-time 
operation. In this paper, we concentrate on the differential positioning with the float ambiguity that is free from the integer 
constraint. The float ambiguity estimation is carried out by the Kalman filter. The float and fixed ambiguities are combined together 
to determine the position in real-time kinematic mode. To analyze the performance of proposed algorithm, static and driving test 


were performed. 
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This paper implements representative FDI algorithms for the GPS/INS integrated 


navigation system. The algorithms have been qualitatively analyzed and the performance of the 


algorithms was investigated through the post-processing of the experimental data. 


The FDI 


algorithms was evaluated by the number of false alarm occurrence and failure detection delay time. 
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Development of a map matching method using the multiple hypothesis techniques 


oF SA* Asa Aq Ej] 

*5)] 8) 3 (Tel:042-483-4072; Fax:042-483-4074; E-mail:toroso@navicom.co.kr), 
ar Hb} 6} 015" 4(Tel:042-823-9495; Fax:042-821-7705; E-mail:snoopy@sunam.kreonet.re.kr), 
west Qa else} Zul--(Tel:042-821-5660; Fax042-821-7705; E-mail:tksaint@ivab.chungnam.ac.kr) 
Abstract: When a vehicle is driving in a downtown area having high road density, it is difficult to determine the road where the vehicle 
is located on. In order to overcome this ambiguity, a map matching method is often used to resolve it. For the implementation of high 
performance map matching method minimizing miss-matching, additional information such as heading error, road connection, road 
direction, visible GPS satellite numbers and so on should be utilized as well as position errors between a car and adjacent roads. This 
paper seioaes a map matching method using MHT(Multiple Hypothesis Techniques) to determine a road driving in probabilistic 
approach. 

MHT is a method to track multiple targets under the clutter environment using a likelihood function. To realize a map 
matching method using the MHT, pseudo-measurements are generated using adjacent roads of vehicle position and the 
MHT is reformulated as a single target problem. Since pseudo-measurements are generated using the topological data of 
digital map, topological properties such as road connection, road direction, and road facility information can be considered 
in calculating probabilities of hypotheses. Moreover, In order to reduce the degradation of map matching performance by 
sensor’s measurement bias, Kalman filter is employed to estimate and eliminate the bias. Field experimental results show 
that the proposed map matching method shows the consistent performance in complex downtown areas, 


overpass/underpass areas, and in the areas that roads are in parallel. 


ES Craig Scott [3]°] 1994do] A}SS HSA AE awe 
1A, = GPS] A] Ye ANS SOB EBS FO ES UE A 
77 FARE $232 YU FqA-ESS(Maximum A 
ae SY AAA 74 FAG SSI BYE Veet Posteriori, MAP)°] AW7} FE ESAS] MAIS FASE YH 
AAS VE A]G. SB B7] AS Fey 74 Hse] Aq lt ay, F Ag Az WAR Peal AE Wage 
Fo] Fon], Bea] 1996\4e] 93S) ADVANCE[I](Advanced QEZZEOffroads APsE AP SE Wa ala ajo] 
Driver and Vehicle Advisory Navigation ConcEpt)=£= 4 Eo] 2jE ja] DRAA 2] 2A SS YB Os FAS] Balao] gc. 
WARS Bae 4) ABS B77) $a GPS/DGPS(Differential 2 ¥ Zadeh[4]9] USA SOS] AS WA tay Fo] Aes 
GPS), DR(Dead Reckoning), CD-ROM #2} 2)= 33+ ja] $23} Qc. USFAS lV] AE Wey gtaySe AVSae ES 
A S32. Ques Haste Aqgte VS 34, 4] aaa Ages 
AS BE AAW We #712 197d HAE DRAA] 2} et 3} yrs Sa qaAsaA don. dae aAIeze anes 
ASS] OS AA Oln] VAME 27] MAS 7ISOB AE (Member Function)\= A334, SE] Bs} Bao] Bal Aq 
3} w]wspo] vase FAA (Semi-Deterministic) V2]S>+ DR QAe FYPG. 1A, Vat JS AIYS ays ages 
AA) 2AS BAGS 9a] OAS Arez FE MANS AE 7] BEM PRA, Beal Paxyo] sae, Musso] 


ms) 


Jo MASE SSAI Ro] AASARAL,. Jt FACS (SIE H CHSES AAaIqAY ESQ AA 
a= Sol, MatsushitaA}[2]4] ASA EAQ AE yawWe AS DAG SS A445 FS DAG + ME EAs 27 

1988'4o] 4) =1 RACS(Road/Automobile Communication System) sich. 

ZEAEMA 744) Apo] ES} AbSAS 144 DRS SkmA 44 wepAl, Oe] et EA] HAS z$o] SaAs}7] FE ERS 

10km™}t} BAS] 4] Z(Beacon)\S Fst] BAW. Ag az BF Aat4 Het BE 2hB DAS + Qolok sri, 

50km HAE HAMA OH 44m BES] BMAQA BINA. Hyd Cesc qe Qay7eS vaso: QRS yy az. 


SF ULB EEINE Alga HAS ae] 7g 
Beso AzgseE GF 7} 7)¥] (Multiple Hypothesis 
Techniques)[5]& AZsto] abo) AMsHE BE EHS ABE 
agae AS Ways area. 

B ERA aletel GS 7a WHS SR AS AAS 
Quees Et VE AM] Des, WES BY wae 
YAo] AH YAS Ha VSAS Au WAs HAs BH 
PPE PANU EA AS yao] Aga ae cpsmraa 
2} 9)3] 8 ApAZ} vpolo]A QaKs 4s] Ae A Wel 7} 
ABSA. 


24S 7a AAS BSG AS WAY 
OS yea 7S Bee (Chutter)7} SAE SAE HA 
& #AS7] HH Reidel] aq ate Ach[5,6). 4S 7 7 
& AY Be] Aol +e S94SS A BY AS 
AS APRS-O] GLE] F-HS BYA(Validation Region)\S 24st 
| 2$dde FE ZAqAH AAAVS AA. Bal wes 
42 s4q2uqe 44 78S Bs] HH, HARMS 
a as) 2 HE SY AVS AAA ach. Hel 3g 
EAS Belsl HA(Confirmed Target)Z} A) HA (Tentative 


Target) °.@ Belo] MSS FAYS Qxajoj Baaqos w@ 


Bago) falda. 
as 
hea “ih 
: 


oft 


2919 AI (confirmed target) 


AB OS 7a ae) Ta 
Fig 1. The configuration of MHT 


4S 7AM AWS AM AH AbSst7] Hele GPS/DR 
AA, 2s} 2J= GoleyMolzS Sate. GPSMRINS 2 
HA Aa) SAAS aAlSsh= Gest ASAI ARS A 
SE DRI Ase}, clo] He TIS a2 ec. ans 
BS AA pars), FEC pers), AA) 2A Wako, ), Aral7t 


AB(w), SE AAP dg), AAA Bat VWaAl(o,), BA, 7h 
AMMA (n, ol. S31, GPSTAlZlo] Ya ASHE 14) & 
at WahS A) Fol Alloa ALAMES| BACH )oll 
a Os SHA VA P)Z HAS 249 MS TY BF 
GAs BIA4[7,9,10]. 
& =(H™H)'H" dp 
P=(HTHY' Oveee 


$0, 


GPS Receiver 
Channel Navigation 
Processor Processor 
Faeroe oar = | PERE i 
; | 4-12 | ete ee | - Pars: Yers 
| channel | | 
Least 2 Sure | 
9 a Sle bese nual | Navigation 
: _ Processor 
t Fal 
> N-State PO¥,5,,09 
| Kalman Fitter | tae 


Dead Reckoning | 
acetal 


| Gyo | Yin Poe 
_ Odometer ' 


9% GPS/DRAAL 
Fig 2. GPS/DR structure for map matching 


292. AS Was 


ake goleMozAS ae, FS, aay AY So] Sy AI 
AYS Yet Ko] JE Quolr, $3] SY BE Fey 
ALG2}9] 917] HA, BS chy PAM AS Fol SL als 
SE SRM Folejo|c. 4S AoA] APBsba] A AE glo} 
HE BME] BAS apHapr] 1yY39} Go] HHI AP LH, &% 
a8 qAaoe Fas, AXE Fe, WEpfaysZe, 2S 


AAs Sl Addy GRE AEH]. 


ANg3sae 


DIGeU 
NSLS KI 


2493. 7= moles] Au 


Fig 3. Digital map information 
21 4# 7279 By BaF 


AS aS Ms SAAS GPSOR AA SYR WFZLE, 
J A(link)tol] ~=Q AES (Projection)S1 BS At 24) Al(Pseudo- 
measurement)® 7} etch. AA, B24 z} SYA} SABalee 4] (2), (3) 
ol) 2}2} UE} SI 

Zrcr(t) = {pl0,60),0,(0,09(0.n,(0} (2) 
Z,(0) = {p,10.,(0- L(G nmnji(O.(0} 3) 

AMA Zesing» PO. Tp + Fo oMy = GPS/DRAIA 2, 
A AA, zz} a] Qab WA, zpaz} Qab Wa, ae} 
ANAAASE A EHTEL He AGMA BA Lee 8 
A(link St 7A ZZ UEP, Azbn Fo me Badr 4B 
Wao ZY VIA PSS Molec. YA} SOae Acie 2 
St GPS/DRE WAS AM ra BWelL]) Qab Serre] |e ef 


8 Zag & AF] Bal HY HS AS Sad Be 
BULT WABIS AAlsy] Ma AIA mpc} SAS og 
acy, 


AM 1H FES SSNS mAS 1 HAL AE SAA 
= AH Zo] PYZHSZ AIA, 44 GS AAs ¥ 
Asc, 47] A, YA, SAaS BASE 7A Age wg 
(mapping) BAS 7}21ch. HS AZ cA) BAA SAS] 
QEeeo}] st 7 sy7} #7}MC. Reid[6]o] 213! MHTE 
aso] AYD +] YE ESE spose de BAZAMS 7 
Asta] Sct. 

Z(t) = {Z,,(0),i=1,2,...,m,} (4) 

Reid= W°O]A)Sot 7A (Bayesian Rule) S AS] SAA] ql 
ot 7S 42Y SS4S FEW). o] AS BMAt SAA 
He BACs cp) gaa ges ge. : 

PIL.G 1Z,O.Z pera] 


: (5) 
= LPO Zaye e(O1 27,0 Wo | 2}- PLO) 


AN(5)2] SA} RPP az sj) 7A] BHA z, 
+ Zoe A AVE AZZ YEH, JQ) AGS deta 
OS3} do] Haat. 
PLZ,(t),Zpssae(€)! Pr gps ae +59; eps rae Or Hi()sF nm Ost (Dh oMey] 
= PLAZ/ gps) ar(®) 09) gps ae(®> X(OnSin ms Eat (D) Py] 
Mo|Alet FAA] ata] Oe Ges go] FA F Veh. 

PAZ | O79) 


= PLE goes ar(t)s OF onside ©.% (2),4¢n.myt (7) 
PQ | 


= PLEO! OP gos1d Or, psa, 110,29 Mel (7) 
PL ()! Br gps1ar()> 5, gos1ar(t) Xi() QP Moe] 
PLB, gpst ar) 5, gpssa(t)! L(O.2759 Mos] 
x PLZ (0/ 27,9 sn] 
Be4o28 fad ANS YASS, Aes WY wpe F 
AA Tes SSgys Aas7] Ha Wezs AE HA A 
J 


(6) 


ASE FE BHGTE A479] PA RAYS] BA (Condition) 


ol] Ja} Eis} Bo. AAS Pa} YPyY, aaa, aazpa, 


YF MN n, A ABS FEE nny 7b ME AGS Aalsto) BW 
AZ AR) Ge Bsr Teach ¥ lo Man, 7 


Qe ado wae, Ged So] VAS SEE, \E 7a 7 


Af AqNACZ SS +EA(a)S ZA AF, URSA RA 


WEE FL 42 EAB )S BA AG. 


(8) 


#1. ARS SES OB VHS 
Member Function On Thor 
a Y 
Te Sin mt =U EE), £0 
29 Vu}. AALS Mala}, 2ztah, AAs] Ala 7}aso] a 
+E YYY Hho] BS USS YS GoUS ZEB Bo 
(9) 
= (58, epsiar(t),m(1)) 


PLeinmys(O)>ey 1), goss ae):59, eps an), X(0)E™ @) 
Table 1. Member funciton using the facility code 
If (S,my¢ =—Underpass ), f,(¢) 
[ee ae 
(Et, a2B2y ) 
Ugo] AG. RZ= HY BH Se7t 7S] Bee ap apap 
ac. 
#2. FEY QRS ol Ss Bass 


= Fe Ao) sSin ms) 
GPS/DR #4 Aw FEAT 73 NEC, ) 
If (Em ==Underpass 2! 4 4} 4), (9 
wea da PES AMMA $A Saye] Zao} a} z 
ao] Ya AE WHEAS AoA. SAE AbAlzbab} S 
PLM () 8 gps ar O58 gps 1ae Oe LiL My] 
Table 2. Member function using direction information 


(4, a=B2yze, [(a,b):a,6& xyaHoe 


re aaa 

ANA FBS] AAPL ABA] 1S A] Apo} Aha za}7} 
AAS MPEGS vepac. 4, ALA 7a, aa 
Aro] Arp. Of 7}oAy Apel a B+ N(x, B) 2 
Ants ELS Bet 4103 Fo] Boyer. 

PLP gps tart) 5) eps 1dr)! Xi(0), 2 Nyy | 

= Sp,50P1,epsi del), 99) gps tart) = IN (OZ) gos sae — EX, BY] 


ayy 
rlo e 
r2 

lu tt 


(10) 


AN] PAS Wal Be Aba 70) Bal 7-AdAFo} 


of) dad Bao] AP] ual HES BITS BaAaala. 


Pla (10°",9 }=—exp(-t/1y) (11) 
ada 1, BF aa ga Fe udue, Ela a 
} Aa Haze] 2 


GZle 5) $M] SARBS Aa 7a QT jl] es 


Vole] As 7] PHs|=S HSS vec}. o] Bee aa 


(12) 


A744, Vy = F089 #4 Gas GRY, VN, = AM 


geld HA Yar AY WS Afo MA HAS vepad. 

AS YM} aAgnly aAa2gAas Azlarea oA, 44 ahol 

Jet A 7A At JS FSAe] AIHA Gol Bel adct. 
Ponta) =f Ona) 58, gure HO) 
X[N(&Z, nei de ~ HER, BYP)! (1— Py)\7*P! FBV) 


(13) 


St, AZ HAL QEBEMARMS ALM ¢-12] A 7B & 
Bee WANE A(x = He Pog - Road () = x 


P(t-l)& AAtesc, 


22 7 Ay BH] 


HS 7A WS 1a AE AY 
Wa BRS quday 
FAS olFoAdh HE AISS UF UeogyIIOe Asay 
AAA 2 Ea) BHO 
probability) 22 B$-S} ool ASML Bel 7) st 
ch. 


9] a] Ae Azo] 


(maximum threshold probability) ©] 72] 


4] 24-8 3(minimum threshold 


23 YE YAY Ye upolyA HA Ae BA 


OS 7A AAS OBA AJS oYAYe Ao] co] Fol ys 
49- GPS/DRAIATS] HA) BR zbdZzp wpolola Fas Ae 2 
a MEYS oo] RSA] Bc]. wold A FAS Met Au wD 
ES GPS AS7 B71 TASHA] SAU DRE] LAP FAL e 
Qe GPS/DREA2] 2ABIS BAS] HA ASA, Beh 
waaAe Age Baad. 

10 Of dp, ) |v. 

&(t+1)=|0 1 Of} dp, [+], 

0 0 1] 66 w(t) 
=@-dx(f)+ w() 


(14) 


7A, dp, Hp, = AA BolA, AAD, do= as vt 


jz, oF THA] PA, UAV OF Bet 0% FRAO 
2 AAS At ag B+ Sse] BATS (White Noise) 2 & 7+4 
arp. SQ wAaas Ass} Zo] HA srt. 


OZ(t) = 2 ie (= Zum, it ) 


PLO ~ Prim il) 
= Py (1) Pim, {t) 

A(t) - yy, AE) 

1 0 Of ep,] [v. 
=10 FT Of} dp l+)¥, 

0 0 1] 6a Vy(t) 
=H -dx(t)+ v(t) 


oe 


(15) 


He 348 


Ox 


2] 0] pj, Lomi 2 = ZS al eo} o]Fo}zl 
AP) AVI WL iE AvMeh HH wine VE oF SR 


ARS We Bt ag Wess] MAPSoR 7AVeItct. 


3. 2a 2 aH TA 


eS 7A WBS BS AE vayo] VSS Brtst7] 4 
3}o] GPS/DRAIAS ALE 
& TAso HI PO PUSS Asta] Fase} Hal, F 
3 


yout Age mae PY DS Yaqsp go] MAMI M+ 


nao} AME Mol7} Baye abe 


3 
7S) Me 2]4o]o], 7AsEe WN syuk so] AVE 
HE Aso] Sole AAS 

29 4 28 AS ull GPsMRaAst 4S 7A MNS 
leet AS HAR] AS wwe Holt. Ad Ade 


uw 


AE SEACH. 


GPS/DR AIA] & 4 & 
JG AAS Sec Adar} AAG Te, way gu a 
AQ QBS SS M+ BAAS] A Deol BssyaA 


ch ES o] APMVoAS] ae Awe! Ye reals) m+ 


njspoj, et AAS aH oy gel a] et mf 


je 


Was Alo] SE esd Abel, ley eb 7) 
HHtACS Sr VSS BF Vd. 
#3. Fa Uy Ae 


Table 3. Field test area 


Agee 


ne Bday gr 
24km 


rea 


ae 


Zs} RE 5km 
Ale ae 6km 


= YA 5] 52km 


meme VN Tatehe fe 


neve CE PSIDR 


134, GPS/DRAIA) SF} MUT AS aay) ae) 2a} 
Fig 4. Test result between GPS/DR and MHT Map matching 


225 QAQs mAs Aas 7A 


a os 


Fig 5. Created hypothesis and confirmed hypothesis number 


A962 ARES SE BH] Ye AE Axold. o] J 
Fe VES UWWEZ FYB FEE AA], 
GPS/DRAIM OH YA AAARAE7E 00] OFC. UpepAl, Ala 
= 3E7; WES AAS BHtrd Aa BS ASA 
ASD, 442) WIS WEA BA ASSO} Azo] WA 

= AQe + AV. 

19g 72 D9 PH] 9a Ag Ado. yey g 
HHEo] Yel Heo] o|Fo]z]= 2] 9A GPS/DRI} ELSI 
AMZ At7t ASL, ESS] SAH] UPSQYel EE} Bdge 
BRS AEA ASA} HAS] QP Vast yao] a 
BAA. | 


2146 A8=S PHd+r7} 43a agdas 
Fig 6. Test result applied facility member function 


meer Map Matching 
wom GPS/DR 


4 8ntr7; 43a aeas 


247, We 


Fig 7. Test result applied direction member function 


442 


He EEAME OF 78 78S OSS AS WARS ale 
Act. AS AS ABS AS WA Asst7] Aso ag 
EAS EZ Bs] ASaact. S|Ab |WAzle GPS/DRS] 


AAS NSLS QA SHA BAS} AS WA Bey Se 


E29) YAY, BEY, AE SE SS Bra pol Hse} 
BAZ, VISLAY So AS WAI ASS Baa + 
VES Basa 

Jaas =a yg aq gs qdaee 7 


Ja 


AL QB Alo] 7}FsVO, QHSZE 7}Aoll 2a GPS/DR 
ANZ4°] ESI AA VAs] QE SHAE Apa] ag 

g@ Eso Beso] oMees Age AME GS FA 
SATh. 


UE a 
[1] Yilin Zhao, “Vehicle Location and Navigation Systems, 
House, Inc., Boston, 1997. 
[2] FUMINORI MORISUE, Evaluation of Map-matching Techniques, 
VNIS, September, 1989. 
[3] C. A. Scott, “Improved GPS positioning for motor vehicle through 


” Artech 


map matching,” Proceedings of ION, September, 1994 

[4] Zadeh, L.A, “Fuzzy Sets,” Jnformation and Control, Vol. 8, No. 3, 
1965, pp. 338-353. 

[5] S.S. Blackman, Multiple-Target Tracking with Radar Application, 
Artech House, Norwood, MA, 1986. 

[6] D.B. Reid, “An algorithm for tracking multiple targets,” JEEE Trans. 
Automat. Cont., vol. AC-24, pp843-854, Dec, 1979. 

[7] D.H. Shin, S.B. Son, T.K.Sung, “DOP Relationship between the TOA 
and the TDOA positioning”, ION, June, 2000. 

[8] Hirofumi Ohinishi, “Map Database Generation System for In-Vehicle 
Navigation System”, VNIS, August, 1994. 

[9] B.W. Parkinson, , J.J. Spilker, Global Positioning System : Theory 
and Applications, Vol. I, If, AIAA, 1996. 

[10] ELLIOTT D. KAPLAN, UNDERSTANDING GPS : PRINCIPLES 
AND APPLICATIONS, 1996. 

[11] G. Minkler and J. Minkler, Theory and Application of Kalman 
Filtering, Magellan. 


The7th GNSS Workshop- international Symposium on GPS/GNSS 


Session 3 — DGPS 


ZAASYS : 
SaulkelS] AASAAALH OSIBSHY AS ------------------------.--.-----------------------+---------------------_ 7] 
=) 


oO 
USA HA (HBFBLEATA HBAS 


ZAaYG : 
sS DGPS SZ2} NDGPS 2) -——————-———------—-—----------___------____________-____--____--_ 8g 
SHS SeANSaSet (WSEAS COA BMeIA!) 


OSSAAOMAS 18-801Z 01S IDGPS PH —————————____________119 
ASE BS pelt, OSU CARH AUSSI, “ABW #SS*BSH) 


aNa*, 2aet 


wend 


ae 13 Haa*, £33", HHS", Wae* 
CASH AASSSSS, “SES 


o*« 
wa, wy 
Ch AFHSS yh) 


= 


2HAl QHYSAAM(LAAS)- AAMSAABS ANAS HA 127 
Ags, SHAM" OAS, SAM YES BAS, vss 


(HBPBOTAUDA, “GMSOIA BSA) 


GPSS OSs MYNSSSAMONAG Ol CHB SA —————---_____ 135 
4 


Ben, ASS (SAH WMSHD, “HRAGAPA Avra pBsreesy 


“ 


1. APHQ 
Oates 
on § 

SAW 


[=] 


3. HWE AlUe| 2 
OAlLtel2 =LI2| Wey Abst 
OBaBAlltel& 


4.NUci ool OE #2012 Dl 
Pelltbcro} We AS 


OS PZer Al 
O Tae 


H 
ial 
i 


ASSAAAG SA 
(GPS/GLONASS) 


=a 


- 1995 BZOlS 
DH Bt- BAe! 

A KI (30me Kt) 

- AQAl2t+2A(140ns) Far ageial 


- Ha zor 
(1997. IMO= 200918 


DEA) Future GNSS 
FSS AH Bl) 


BMA AlaAS 
BASIL -120%S, 
BBSy! - 4028, 
GSES! - 30°, 
DSB Bet - 20%, 


AISI] - 204$ ) 


74 
i 


aval 
oS 


2S OM DHE 
SaAS OCASSi = 
u 


Pret Tress 
Of 


SCH OH wt Ml oI 


eS 
= 
otLhe AH ES) 


(S32 20018 FE 
EQIIAHIA ALA) 


2} O} I BEA et 


2} HI 


@ O22 BYa/e+4 SSALSSEMAAA BSS 
@ AMSHAAo APSA R/S B/S A/01S SH 


> 


AS si(IGEB) 73 
3 26 


AHS DEM SS 


YS 10 AS"HS FHA ASSAANAS USNS Bass 
sal pa 
e UMNSPISUNS AS SHABNS SEAS 
@ SEHYAAIAY AASAINSEDS OA APH He SBSH 
SH BIAS) WAYAS MAS AS 
fo ers” | aa a 
2 
———— a i ts Specie ee | 
mB >U AASAANAE OSIAISA SIA! WEAS Bk | 
@ ASSAIS WHALE AAW SAS Soto 
@ S2ili2lo| SAE o AIA SASS 
e YEE WI, Sa 


J 


ele AASANAB BIse WHeas 


ASS Be WAIIS ZA 
ASSA SHISSASOS 
27} SBI] AASMIIS He 


OAaPAZyee2 
© GNSSH EA SHA I/E JHB 
e2 He SESH WEA! 
AB+A° WANES BS 


OlS (<tus), S92! Al2tulw BF! SI|(<50ns), 
BKatHel, AB AALS +e 


S SQUAY, +24, 


- M-commerce : CNS, AVL, 2A 

- L-commerce : E-911,q1H,A 
O dae 

- HASSE (AA, 4A, O24, al Boat, 2 eS! ZA (ZAA, hot 

AA, BA Bela, AR BAA), SB St IS, GMDSS, ECDIS) 

Oo As 

-GPSS 0182 We +357/2 AS, GPSS O1Se As GQUAZ 
O GPS + Al 

- GPS + Wireless 71, Low Power & Small Size7|= 


HOuas, SxS WOH 
%,mobile yellow pages, personal tracking 


2 
Cc 
q 


Salers) ASSAAAB SIs WAS 


fall 


mee | She | WSR 
CoM) | Cs) | CON) 


suuad AASSAAAE Sle WeaAs 


as) 


e HHH OS AAW AMS aa tes 
- Ble ISN, AMA HWLAS Be 


- B89 BIA SP ASMA MSS, SI ABS Se 


@ GPS} ASSHAABS SYSS A 


e ASASANAB OSAAD SAS 
- S/ADt SHAM SISLALSH OFA CEA HAPS SX 


- GPSASAUANSS AS (AUW DSBS 


S2lad ASSAAAH OSIANe WeaAS 


es 3. AME AILHEI S - ALfel 2 2} 2) G14 Ab St 


e AASAINSI OBS 
O| 52] GPS A018 SSIS as, GPSS! HC SH(L5, L2+Code), S/A Off 
#82] GALILEO B28 #4 
2A] AASASA OSSAS BHISA SL 
GPSB URS SASS ISS YOR 


Ass = a (AMA Al&) 


© ASSAAABS S384 
- DR SSI: GPS PEMS 8 S/A Off SHI AA AF 
- R82] GALILEO7} GPSS SCH MIA Al OL 


Bl AS=z Io Set SA 
HeoOls AS AS 
OAL AlLtel2 aa 


Pelee] ASSAAAB OBIS WAS 


#8 OSA ASHNA AMS Ce 


e =a ee 
- SAAS USst SOF: SS 2 HALOS) DGPS PAIIASA 
- ASSlOS OUActSE I] ZL] iSO HIS GPSS ois 
M2A, SMA Sv 
gSet $0 Say, 
A2tS |, SA, SS, 
AIALA S| BS) A 
WEMNSHPS D Al 


@ 2S 1 2UM 
GPS#S2 72 
ArAICH GPS = 4 
AAS SAM SMS 
GPS# 419] JHE 


‘ealuatel HASAN AS olgvie Hes wet 


| |S. HEEAILHeAI@ - alae? | 


O AUA? 2 
e@ AAS 
- GNSSOISA YAS SAGAA AA 
(HAS) PS 


e £ue 
- GPSP OSaALI SS2SO TE 


@ NSH 
- XIAMIKSI, SUE) 1S 
AUS SAIS 


- LORANOIS DGPS S44 


oT 


Sas ta 


- AlLte]|2 12] GNSS(GPS, GLONASS) @ =! 013 at 
- RIA ASAZOILE RAMA SAA SS MA St (GPSS| MAIS) 


FOolS SAAC! SY AA 


IZ Xl 


Petters ASSAAAB OBris Weas 


| | 3. HEAILF2AI@-au2e3 | 
ee = 


O ALt2123 


eAASHE 


e =um2 


- AlLkel2Q 2+ 24 GNSS ZAIS 


@ SHLD 
- MSAS2 AIS, EGNOSZAI=, 
GALILEOS Al S/S ABS 
FA, 2S 
aes 2Aap 
ASALAA BAIS 


- 23 GNSSO JEGE VSAOReS SO 


- CHEO| BI} lato GNSS YS = 
a | 


ASAOl AAA OS 


ion en ieee nd 


SBI RAS ANSIBS 44 


C= FEMS Js RA 


, 
» 
a 


eee 


Selle] ASSANAB OBI zs 


—76- 


| 3. WHALE - AUIS 48 


- SUS 


AAO SAISSUM 


. rls ee 
-atuaigas seg AAAS 
+ BE IGSO Me 


IGSOW S718 SHBSA> Ue eae ee 


; AlUeI2 5: 
-AUGLI+ SA AHRSSOUA SSM BH 


- Galileo », Program 


iP 
i 
I 


ee HHetAlLbel@ Hla 


4 ChSHAS Deal 
OLS At Atal At 
a2 ct 


{ 
i 
| 
ot 
{ 
‘ 
! 
i 


© null 
saagy 


O AlLtel2 5 
osasHe 
-AM AA SNASS 
ozaAet 
- GALILEO 2278 820 
o FHS Dt Al 
- GALILEO =AlSS 


A#ene8 


“ealuyetel a4 


SANAS ogee x 


we 


ea 
= 


+ GNSS BMS SS 7e wes 
a3, BAS ot ne 
. s01SZ : 
epi] 8S FCS SII ame 


‘Et ENSSOISAL VAS 
RAIA Soe 
21 SS aA PS 


SAE 
Smee Ot 


QIRGNSSH NEStHets 
BO ASAC! AME AB 


AIPBASAAS 
we 


AAS 


Ads Jed 


~ SBS FS Ala] We 
“| = SxBpGPss! SIRS Atul 
“KAM SZA2S IDs (S/SA) 
cist] WSK OLS Al2+S 21/LOMAN OLS 
*OGNSS See 
- BAAS SSA2tS 71 Ae 
= GPS/GALILEO/AI2} SBt+ dale 
“NSASLATS, EimsaNe: GALILEOZ AT =/ 
ASZAZ RAVES 
“Mas 2a ap 
-ASBVANBAINS 
* 1690 AAs BMSSAP 
- WRAP NEA? Awaz sea 
SAAS 


AW ASS zene 
0 


* GALILEO RMSSA 

- Integrity SatAuIA Bol 
- BMRZS 4G @ +41 S 
- 4248 Sane 


SA ASHES 


Sellers] ASSAAAS 


To a ALelool ge Qe 


Q 
o 


Nae 27a 2 


QCA FAB 
VSAAINS apa onwen eaee ONE ay 


DPS HAASE IAS WH 
(20004) 


e AL at OR / AM 


ROSA RAR 


WAZRS 


GPS APAPREAR ATION 271(1999%3) 


WAR H 
AI74 Dt 


GPS BAL B® 329K 
Hi C1g9gtt) 
OGPS BHAA BUM AAR WH 
¢ 20004 > 
wre A FAPAARee SBR VD Et 
ABaAy AMY ASS Sse We 


Aes 


ABR FS 


ARM2AP GPS QAIZS> BH CISD 


APR AP. CPS BARR BIHCISET) 
BIS HII Sas, Ew. wi 
CHV RREAl FAL 
SS A MAE SAL 

28d ON DIOL. 


MAPA NS bine Het 
PSSALSaCay) 


Bey z7eves 


ALA Sol GE F QW IAS - BI 


© D122) Aral 


~ GPSHSA ARO CHSt AAAI A(IGEB, CSGIC)I} USA Sr, 
2 ADIAWAM HEROS SAI 

— 19944 DOTS DA (NDGPS FS a/s 

- 199644 GPSLIES SRSA Het A+ (NRC) JH4 | 

© HL2SONECEUUS ASA, HAIISS sad #0 

* YEA SS A588 1s 89 sro 1 | 


r9 
rQ 
a 


30 min Trimble 
RON Fates 


NOAA Profiler 
Network Site 
with GPS & PSOS 


Constellation 


USCG DGPS Site 
with GSOS 


Pellitetol 


a 4. HEAL CO OE FRIIS - HAC GPs* ell Jet a 


ZEMICH GPSAAID| JH Yt . 


O Wa 
e@ O|=2! GPS&ist 
- 200344 Block-lIF GPSOl Ai L2(1227.6Mhz)0l CASS AS (L2 Ol 
- 20054 Block-IIF GPSOll Ai L5 (1176.45Mhz) As a S(s Qt 
e@ S282] GALILEO 22728 
- 200044 5,62 WRC-20000 Al GALILEOS 4it= 2S 
- 20004 122) 0+ X) IH At! ZAI, 200744 Full Service 
e SU 247/82 
- SASIA, ea A, DTS 


- ASB AA 


O AAI GPS | 


@ 11/12, L5, Galileo = AI} 
@ MSAS, WAAS +419] 
e Ad, SAA SMS +A) 


: ney Darts requesht 


Salers AXSSAAAE OBIS WHeAS 


| | 6g AUIS SQW IS-namagazse |  «/\a | 


ASwWols AVA ASI WS(AAeS) 


Ow a 
@ GPSO| S37} 
- GPS O18 2 GPSO|S 7/3 We 
- GPS2| Ol= Q! St} 2] Back-up JIS 7S 
@ 0/2] GPS Back-Up 
- ‘8041-971 GPS IN Fe (sole means) C= 0/8424 
OLA = &t 8 (primary means) )1'8 CS AS 
- 20008 2H HA) A101 LORANS! AIA AO! SAS(GPS/DGPS-beacon/ 
LORAN 2&9] Jt 
@ cYAlOt 
- 2000 102 Oll Planning for Global Radio- navigation (Integration of Satellite and 
Terrestrial Systems, including Eurofix) MOSCOW 20003 4 = 
~ EuroFixS NSotO ALO ASO st AMIAS Da 
e #82] Naws ose GPSeStay 
~ 19901 CH 24461 LORANA! SO 2S 2! DGPS, DGLONASS4! 
QU= EuroFix ttt =A 


fosaualorancesee 


Korea LORAN-C Chain GRI-9930 


Pohara Contra: Center 


Pyon: 
Moniter © 


Tovane-3 Staten 


| : 

Pas! Pohang ‘station 
i 

| 


[itjima Statice 


Sates ASSANAHL OSIAS WHAS 


t 


4NUACH OE FQHLIS-AADAaIgntsra | 


‘oO BuPcse 2 Zies0e 
| PCS MHIATI 
| 


Fl 


J2 
jo 
Qo 


SDGPS ABIA AS 


ee 
OHS ou 
ae ers 


Saud ASSAAAH OBIS WeEAS 


|, ONATWOWS 

|) @ LORAN-C 

a - Saat SY AC BANS, 49) SAIS BQ 100Khz HAT AS 
- GPS back-up ISLS HA! (OlS 80,000 ACH #SI| SA, SA! EuroFix) 
- HB] (9930 East Asia), BE, WAKA), SAAAADAAOH SB 
- ELSE APSe! OF Ot OILED MMIA IHS (2,200Km) 

@ 31 DGPS 
- =HAl°]2s}(IMO) #2! O HE AMAR 
- 1191 dGPS D/2== (Reference Station) 2 Sol 2AI= PSS 
- DGPS SAYA : MSKBIA!, RTCM 104, 285~325KHz, 200 bps,300~500W 


dl ts 


&3(RTCM-104)0l SAMA ASS Ib, GPS2I 
S0lE th NAY OS, GPSAUSS SSA 


3/+4l Is THR, GPS OS 
- GPS SI SSANE AS AYA JS 


Selltes] ASSANAL OSs WHeAS 


|| WALCO Ge FQ Wwe am aes | 


(0 Gate Aso 24 
| 


@ MSK: DGPS 2S2aA2 (RTCM-104) 
@ ASK: A4SS (Hae) 


<SAe> 


Aztae AW [fartitude 
—> 
Gos Modulator Outsut 


<A HSS OSE B> 


MSK 
/ Soectrum 


Sellicts! ASSANAA O BISWA 


a 
S 


~ 8i- 


“~ 
2 


|] a, HEEAILE2| 2 Of CHE = 2 EET 


SHSAAS SASAGASS 


-#HSAAS SNS Le 


O Ha 

e@ GPS SU0IBS dela Ba |’ 
~ DISS S/A Of 01S, GPS2] SAA (Integrity) BAS Fle! OGPSA| SLS 

@ 2°) MSAS, S82! GALILEOS SS SOMO AS SAIS SB 


O IASAIS FAS SH ps8 Seostat sant y Satelite GPS Satellite, atc ij 
“AS TT A : ee Signa, = | ; 
2 Dt ; Hic pana hier ener gee / i\ ‘ 
ecpseza2zNazes ff mst / | | \ \ | 
@ ASHEN SSI Ae Bi” gee ely eee a 
Olescrct fete | \ | beter] | 
eras i Station | eK] Station} |; 
ee 22 ee 4 | 


Ground Statran ARI 
NN 


Inkegerty and 
Defferartial 
taifuemat- or 


GIT Network mg t8! fore 
4 oy 
Concept of GNSS ‘other nets ork 


Sauety AXYSANAE OBIS WeAs 


To] a We AUelool GE FQUBIS-caumeozmes | 


GALILEO =H 8S Aa 


O wa 
e ops WSee2 S8el GALILEO H2NS FA 
® 22o] 2S OSA VACA SSNAS FAI HVS Jel |: 
Oo 2288 


i 
| 
| 


@ GALILEO 4 Ut 


~ WRC-20008} 2|(2000.5-6)O1Al PIS SPS 1iMhz FUE AS A ItWWS B7Mhz toe 
(GALIELO),24MHz(GPS L5) -> GALIELO I+ D)MCHS > 102Mnz(B7Mhz + 15Mh2) | 
2 


@ SHee 22 
- AAW BHA BETH (ZI SH FO, WBS SHH) 
| - SH BOAB: AWS/AAO}, APIA/K LAO, LICH, OAcHS , 
- FHLECHE CSAOIA 3008! FHLICHS HO BA AA(20004 3B) 
- 9} BOARS = JHLECH AMA, OAH, LSMOl, SOS, OFVAS OS SF, | 
WH SASHDIS YO, YUAOWAGSS, MAD, SS, OMI! SIS 
O Selltetel Az 
@ 871 BOS A 
| | - GNSS2 2 2/NATO 2H EGNOS #O2/ESA HAUB/RHSSAP BH 
G 


@ Z2ilick= EUSAAA SO SAAMI IA SHAE SST GALILEO 
DEH S2zaAISSH FHSVDAGAS ASS Laat ers 


S2UAay ASSAAAH WSIS WHeAS 


ASB We Wt ols 2121 
fe] AAO VE Allteies 
Ses AS Ne CH ALDI 
BAIA 1S) ZI ee! 
AIAAOlD AEs Ses BAe! HVS AMYeSe Sa 


aac 2a 


O AAS SHS vss os Tig We 8! O18 BASS SoH 
BO PS MAS VSHSD cael Aeae: Sa om 
~Divt Dig Deo CHS! BB SAS Ss Be B® as og r7+ 
DtSst S83 Jee wa 
-PStet @ BSD Bee ASB Avlstol wsiwe weer 


eal gl Ale 


CI RE WHS All HIS Vso BSBas Vsrsr71 Heh 
Ss BSS! Tel SICH 


OA Sia Ber wrey, 
AZ SABA VSst 


AO] SB! SAB SHS Ata S 


HRS cersi@ 71x) e 


Sauer ASSANAB SIs Wtas 


age 


AGHA APMS/SSAIBA SLAB 
214ol RAMP HB, SA, OBA ART 
RB AHMRA Swe APIS AS ws 

BS 1), MASA) MSAPH AAA) S 


So 


28 OSs S/S WY, HID. 
2 iMeI 92 OlSst 219) WY WY ASAA 
2 QAl G BES) Wo Setz zl 
BEY MAS, GStarvs 
S0lS BPeAlelSs Sass DSi rH! 
2018 ASHMSHWWPAOSEE eae 
= ASO MME! S) Deony 
SOHZSrHe B PARMA a 
ETD) 
SrDD| SNH ty ALSO! wey 
CAO BHI HS. SHEA DW HSA AS 
CABMS[ HAAS O1f AMA So CHE 

LAI A&11S8 OlBst NARS, NARS Pel 
@ WSS ANSH APS 


SY do Bato 


4COm E WTY 


J 2+ 4o 40 


ce +40 


Suu ASSANAE WSBIs WeaAsS 


ofo| 72 
1x 
a 


oe 
> itt 
“” 
oy 
FQNZJ 
ig 
ong 


1O¢2M0 5 HL 


0% oF 0202] 00 


fOSt4y 
S 


a 
> 
HOI2 GE 

on SIA 
vie 
sos 
Hur Srato: 
HATE 0 EJ 
Pal 
(G3 


nf OF 2 


Nm 


~ 020 
Es} 


| 24 
nO 
fQ) UJ 


0 Aeeel 
Dee 


ane 
Ass aa 


ol> 33 2410. 
“O28. fier: 
ZA 


he 


AC Oe HO 


FASASBAUIA BAS Ai Y W218 


= 2 
0 ESF SEUG SS 


Ariat “G5 SS /TM NE We oA 
seats | ENE HEQESM TNS 


Sales AASAAAH WSIS WAS 


— 84- 


ETSA MSH 
; p SAREE sist Na TIeNe 
0 gels SIUIGM2|e. BUAtOWE. AMADIS 

0 S24] BS p AMES SMHS IIS MHYCAIEESII 
WSsAesl| DH ¢ SADSsplse eS 
X SBS} Galileo, VSS 1SGOH] BAe 

- BES WS ARS 
jo GNSS SSIS (MIA ait) wet 

- HASHISS OSs! Viol Y at 

PASUNS OLB SEH 


u} 
at 
3 
3 
=) 
» 


palate ASSANA WZ Is Wer ae 


Integrity Monitoring er 
MSAS, EGNOS) 


Paes Pes 


gaHor) 


-S aa a@Aizt ve (SAS BAZ Shh) 


-DEPSS AN] HWS 


-ASBSS USAF - 8818 AE 


(QUA® A) 


~-LORAN/AWSKS BA] 2357 
OSS SAE 

-SRANAS SBAASANSHS 
DPS 2 ASH BBStN718/-88 


2h ee WE 


KOS ASA S407) 718 


nH. GPS Meta 
AS Aare 


SN MS 


(LI/L2/L5} 


@. ASSA 44 
88 ASHE 


-PSBA LEV AST 
8 4471 4 


MSHA OY a $s 
(SUAS x SAH) 
S518 ACMBH SH 
SRUVSHA BA 


SENSAS 


weet 


Pejltctsl ASSIA AR ol2 71S 


The 7th GNSS Workshop 
Dec 1, 12 


Park Jae-Hyeon 
Director, Aids to Navigation Division 


Ministry of Maritime Affairs and Fisheries 


oweFS DGPS AazAl 2 Ole 
SA AHlHA = DGPS # 
NDGPS #4] 2 o 
Loran-C 22 

Eurofix 344 

A het 8 Al A 


-ow#ASDGPS S 


Differential GPS(DGPS) | 


eos 


DGPS HAA SS 
(BSS, NAAZIE SY Bs Way 


OGPS JF OGPS MSA 

1, BB HAS IM INA WERS Sil . GPS ASB eal 

2. GPS SB +A! . OGPS HAINIB +A! 

3. GPS Satz Aldt . GPS APS YA 

4. DGPS HAIN YS . SS ARAL AS B GPSAWIAC AHS 
met 


B= He 

-ITU DGPS JIFaD 
-IALA RDFOIS OGPS2!D 
* IMO Resolution A.815(19) 
- IMO Resolution A.860(20) 


* ECDIS, AIS, * 


-Hetg, ager = 
* MTS (Marine Transportation System) 


BASAS 
> 98 - 20008 7742 SX 
* Controll Center (CHA, 1942) 
* Reference Station (89H 2) 
&ZOSC, ASS, ets, HES 
4c, #28, 22,355 
* Monitoring Station (89124) 
% AFL SOW OLAOIAT Al 


Tiandashan 


x 


Taiwan 
a 


Ml Al DGPS BS 


mw Dual Coverage 74 


Bsa 


L—J 
peas 


= 
0 
oD 
Ki 
Mw 
=] 
Ss 
KO 
iD 
ro! 
oS 
Mi! 


st 


WES 2H AX 
> BOSe AAS IA OS At 


- Sol-sal AA A S20) 


— 
=) 


2 Ie 


EX} Ol 


:1,71481/ 984,000& 


DGPS Korea Coverage 


—90- 


H2FS DGPS 


Mi Radio Beacon(283.5 —325kz)2 01S 
* ASS 200bps / 4 300W 

- USHAIES BA 100dH21/100M/m 
@ Format : RTCM Ver 2.2 

- Type 3, 5, 7,9, 16H BS 

@ Accuracy 


° +H : 25cm 86%, =4l 70cm 84%, 3D 1.4m 99.5% 


@ Availilability 

* Broadcast Availability 99.7%= 
* Signal Availability 99.9% = 

* User Availability 99.8% 2 

* J}== Dual Equipment 


 Intergrity 
Hoe stHAl : 12.6m 
* Time to Alarm : 200bpsOIlIM 2% 


SAMI H = DGNSS #2 a(1) 


M Ol=ONA 2000. 5. 24°48 SA BAI 
© 100m SHS 20m OWS APS WA 
SH AUCH, AA ATAY OK 


-adcls- dea 


magc, sas 
* 20034 2% Gl2tdlS(L2) JHA 


*2005H Bea Je 
3A} A! S(L5) JHA 


al l= 
L2 


GPS dle Ji 
Ine 


SJ ALSO) PA 
(NAJIIR) 


L201] C/AADE 3} 
(IIR-M) 


L5 


L5o] BI2kS Ale 
(IF) gee 
1176 MHz 


(las 


P(Y) 


1575 MHz 


Block IIR- Modified Satellites g 


L222 
*P(Y) BE S49) B43 
*HT2] CARDS BH 
awe se hee 
CIA AS #7} 
S28 Pp(-yDeS 
SHH} Se 29) 
N22 M,Z2ES 27 


- 22 dla? ZEB A/S (Mon L1 and L2) 
-1H1 2 Aly? Bl2t& AIS (C/A code on L2) 
-HWEY GPEAFYS NEY S454 (HAMAIY HFF4Y AA 2O/) 


- Basic IIR Launches i wh 
" (JPO a hata Profile) : ioc 


earth 


T A} ist HR SV RetrofivDeliveryaunch 


ae 


A = 


Ast! nF ‘Delivery/Launch 


T "42th SV Retrofit/Delivery/Launch 


"Legacy Update 


-—-93-— 


: ‘MM, 2E2 A 7e FF 
L122 er *C/A BE EF 
M; 252 W7e =F A? S240) HEF 


AF SAH Wat AWE Seiad 
SLi Ol At Wat al SHAE ABAD 


ZAI) EF a7 


L5 Ale 
“Ar? AAS Ae 


- 2H2 dlr? LE A/S (M,on L1 and L2) 
- 2H Alf? G2 Als (C/A on L2 starting with 
Block IIR - Modified and new civil code on L5) 


-SY ASS YIPE BISA 


I Transition Ops ] 
ntti VE SW Del 
H OCS Version é | sig 
: : Development Transition I 
: SV 1-6 Modification 


| SV Deliveries he 
pi aa perianal 


Block IIF | SV Launches _ if | JPO Projected 


(eps oe Xe - ea 
_ Last IR Mod Vehicle : 


| Launch Profiles : 


SV7-9 
tL “AP ; 
| SV 10-12 
LL /\ Prod 


OV RY ALA EO Ch Bt M_2 Spot Beam 
22° : 
we fo | + 20db (-158 to -138) : 
eOXeo more power in 
ope ye ps new M-code Ls Als 
eHAFe . + WF WAS MS 
“ALB Al 


Cite f 
CANAVERAL *) > —., 
cw. “<Se 
“. ASCENSION , DIEGO. * 
> GARCIA _- 


ea 
Ne 

7 ~ 
ont 


i i i 
petitive Trade Studies 
(Two contractors) 


{ 


Li | | 
Sn ie ter ae | 
| (On-Orbit Test 


Systems Arch 


Rats Definition | 


Decision for Production 20 dB | 
i based on Qual Vehicte Testing | K 


| 


i ALL/\Prod' 


GPsuiProd] sd 
/ : / KicDere 
| po sviews 
iA 


NOAA National Geodetic Survey 


DGPS /GPS WAI tH 


Al: 2000. 11. 7. 


1 Sigma 4099m 
Positions 7952 


1 Sigma: = 0.134m 
Positions: 7760 


2 Sigma: 0.268m 
Ouration: 23:49:29 
“Updated 4428/00 12:41:17 PM 


3 Sic 


2Sigma 8.199m 
Duration: 23.4651 


Updated: 11/8400 12:38:37 PM 


3 Sigma 


DGPS/ GPS SIAl aH 


Oe} 
2000. 11. 7. 


5 1 Sigrma: = 0.149m 3 Si 1 Sigma: = 3.274m 2Sigma: 6.548m 3 Sigme 
! Positions: 370 a 368 Duration: 01:01:32 ; 
ere «Updated: 1177100 £3345 PM 


GPS Scatter Plot 


~-L10 


1 Sigma: 0.137m 2Sigma: 0.274 3S. = 1Sigma: 1.314m 2 Sigma: 2.628m 3 Sic 
Duration: 01:03:14 ft itions: 38 Duration: 01:03:27 


detest T7700 T24S17 OM: 


4 Sigma: 0.189 2Sigma: 0.378m 3 Sig 1 Sigma: = -1.135m 2Sigma: 2.270m 3% 
Postions: 389 Duration: 01:04:46 ; Positions: 352 Duration: 01:05:12 
ala ala al : Updated: 1177400 1247:52 PM. ome ipdated: 14/700 12:48:22 PM, 


mw 2000. 5.4 GNSS Flight Recorder Approval 
Commitehl ASAS OlSotO] SAet By 

100% 32m O|LW, Hat 13.2m 

99% 28m OIL, Be 13.0m 

90% 20m O|LW, Be 12.1m 

80% 18m O|L, Ba 11.3m 

70% 17m O|Lt, Be 10.4m 

60% 16m 0], Be 9.5m 

50% 15m O|LU, Bae 8.3m 


SAI + = DGNSS 4 =e = Al a a) 


@ 2000. 5. 12 IALAQ] AD 
“SAMA AA LAA 20m OLULS WAARSLE IMO 
BSlA815(9) 01 Al 95%2] SSAA 10m O12] ATASS 
FAot7| Aso GSia 20] OGPS AS SA Sel 2 Ha 
SE SO HSS WEAAS 
Ret At ase AS 
\(Integrity) JHA 


SE 199.7% 
5 99.7-100% 


99.977% - 
99.999% 


HES SO} 

-82 We 95-99% 
-ZUPA Pach es aAIDdNssS 
BBN Se +2 AIT SS 
-32 tAlz2t aIcgnt ses 
-A2ta Ss, Dat AAI ss AAA es 


NTIA Special Publication 94-30 


EJ 
OHM 


Cre GPSO|S AIAB 


S 


Integrity | Availiblility 


Accuracy | 


GPS GPS 


GPS 


GPS GPS 


GPS 


HU) St} St) OHM 


NDGPS NDGPS 


| NDGPS 


TO) 4>-| 098 004) Fe) 


2a! 2) & NDGPS NDGPS 


NDGPS 


2] AN £2] 02 | Of 


ee ze GPS GPS 


GPS 


(M1 eS) 008) £2 (2) 2 


la WISE GPS GPS 
2 He Bell | 


ox} TU 


GPS 


_NDGPS © 
_NDGPS 


GPS 


WASSILASS | WASS/LASS 


GPS 


WASSILASS | WASS/LASS 


GPS 


Om! Jat) HJ) £2) 2. J 
JeH | Hl |S] 02 | D4 


| WASSI/LASS | WASSI/LASS 


WASSILASS 


oH2t DGPS AI] PS 2! ANIA AMS (2000) 
oH2t OGPS2Y] Double Coverage +A 2 (2001) 


DGPS 0IS#A0l Set 
-oh2S DGPS SSZ#ES 
*NOGNSS 74 
*DGPS > DGNSS #C 
*DGNSS S# OSS +419] WY 
*DGNSS $87|3 WH All| 


=i100= 


= He (2001) 


Ea SAO ae Os As 


¢ = 71 OGPS2] Aas 
* NOGPS S2 |S Jt @! 
* DGPS2I Infrastructre: =Olw 


DGPS2] NDGPS?3Al 2S 
+401 AA O18 Its 
ISNS0! ABI SsO!| AS AS LF/MF OS Its 


D ITS, SABES(PTS), SS2tcl, Emergency 
Response 
st: JA, A2es 
> NAASAAR, HAScl, AGHAS HS 
>: ABs 


=O 


= 102% 


BxXId SA 
* oHAFSDGPSS| A 
XESS AAAS A 


SOGPSIESD AAS WSOGPS 
JO] AALZ2 SAM Sat AAPZANSE 
BOA QAAl AS AMS 


PAO ASHE Beet wane 


mM NOGPS 73a 
et] NDOGPS JIS 
NOGPS da SAS 
Loran-CS$3 


Aale peo a 2 (2002 — 2004) 


i 
Se S HSSA HE 


Ol 3-494 WE 


-103- 


* DGPS #419] #B2] Sash 


Seg sce 4a] Wet 
* DGPS 38 J17) WS 


Chayka, Eurofix, GPS 


¢ Loran-c, 


Galileo, GLONASS, J] €t 


-104- 


= RAAL SE AAO 


HW OBSLLeSA SAA Hol 


SF 
* GLONASSS! 3a 2801 
* S38 °22 WADGPS > 2! 


DGNSS Format  $} 


¢ & Ai RTCM v 2.2 
*RTCM v 2.3 


—Eurofix, GLONASS, ephemeris 
new Type 18 - 21 


* RTCM v 8 integrity (CRC) 
*RIK 
¢ interleave, embed, new station 


05> 


LORAN-C 


mw LORAN-C 212 
* 90-110KZ2] BUS OSs Saad 
= BAcl PASASAl 


- O12 A: 1,000Km, AAILA : 90m 
- a AM sBsHyoes gs 


cy 
3 78°22 FA @ OA o 


S045 ASC Ae 


* Scilketol 3 


Siem ain? -~ 
meme ae 
te Se re 
eo =~ ’ aed 


= 106 = 


ae 
So, 


a 


=107>= 


LORAN-C COVERAGE 


free eee! 
} 


fosraee 


“KOREA FP ORAN-C CHAIN 
heal CUVEE RAGE RHE A 
GRE USO 


wee Eee a he 


i 

H 

i 

; 

i 

H 

i 

! 

! 

i 

ere € ! 
ee i 
_- ULM wage be te : . | 
; 

H 

i 

i 

: 

H 

| 

; 

| 

i 

i 


+. 


es Soe ae ‘ 
POCA RY CONT CHE RR ot DPE > 
TOM. Atte es FS RNR: 

TRAMSHETTER ETAT TOES 


% Wexurt dake 
COT RVL/MOMETOR STAT TOK 
*Pohandw Gantrok tnt tors; 


ePyongtaek & Kan jo laah 
dt 7 oa ifur 


=108= 


Eurofix Ji & | 


mo 2 


* Loran-C2} DGNSSS 28 
SPO A2HAS 8 tod 


HOe Loran-CéAls 


mal Wa 
-GPSol He AAA AAaAo Wet 
SH#42 S)2, SHAAN he 
Sot Loran-CS A221 Atl SABA 

* SARA Loran-CAIAS(NELS)S SCR = 
-feEAO0!, AOR, $2, BAA, O 
-FA AAS | AA OMA SHS Al 


ne 
> 
in 
a 
ng 
rw 
In 


Eurofix 5 @4 


gS a 
- 22 OSS : 191 SAIS 99.8 ~ BI SAlzS 99.9999992 
* Accuracy : ## 1.5 — 2m (95%), #4! 3 — 4m (95%) 
-AWIA SSBHA: SASOCSHE 1,000km 
* Integrity (time to alarm) : LAAS mode 4.5 - 6% 
RAAS mode < LAAS 
* Continuty risk : 1x 10-4 per 150seconds 
¢ GNSS 2! DGNSS& Sel GNSS} SSolll RE SSA 


=I}M STA MBIA AB 


SIPAOL Al 2H 


B Additional Time Modulation 


Late 


Early 
Prompt 


148 PPM modulation 


110 


Eurofix Coverage 
(NELS,SELS,Chayka) 


Eurofix Coverage 
(USA) 


a 2 


Largtude 6 [neteret 


4 
a 
4 
4 
Sj 
4 

aL 


ee | 
cir 


o Vig ne 


Loran-C2} GNSS Ul wW | 


* Loran-C _—- 
* Terrestrially based <> 
* Low frequency +> 
* High power +> 
* Repeatable accuracy «* Absolute accuracy 
*High penetration <> 
¢ International <> 
* Robust +- 


GNSS 

Space based 
High frequency 
Law power 


Line — of — sight 
National 
Vulnerable 


¢ The perfect combination!! 


=1125 


NELS2! Eurofix 3 24 (1) 


A, LEAMOBHAS SASAA 2S 
B Ast OA: O.7M HS 

* Eurofix S24 YSS Ao 20004 JS4 Ala 

- '99 Het BrestOll M MIME 42 I 

*HASse SAS © oran-C AISHA ASA SS IH 
(+ AlD] MAA SPAS PA AS MS) 

° SAID]: GRAS) Be es We Aw 
~ Disk Il — MAN/Vicon—Locus—Reelektronika—Diginext> 

Gloria~Novatel>Marconi 


* NELSS} AS 2] SAS WO Al 
- 44, 2H, © cClOKE 3} 


g AIA Sa AS 
* Loran-C/EurofixS Al=!S2! Galileo? 2st 
~ GALA 
— Helios study 
* Loran-C/Eurofix2} EGONAS &t 2 ét &1=2(ESA) 


=Lle= 


* AIS Loran-C APISS SBA 

- OSHA FSA —- OFAO - SO 
AAgce HA 

- ZOH2o SASHAABGO ASS + WE system GNSS 
oO Ost QP RA AR 


- HoS2°2 VSB system 7A Js 


“RAM ASA! HelraS 
* NELS #22 OFAOS PSHE 
- FERNSS AIS chain slo] G4 
* OS0NS SAH! Eurofix ASI 
Al HAS Sst Loran-C AA Hast 
—- Remote Automated Intergrated Loran Ji# 


— Aji, StH, Software2l JH B@ Az 


23808 de & AI Global #HAIABHOS OA 
- 2k YW AABS HSI AS SM ala 


=114= 


MOMAF 2] LORAN-C 444 


* AASAANAHS! Back- up AAGWoCe é 
AAs WA 2 RA 
~ ARAFS BW Is WAS ASO AA Ss 
* Eurofixet AB Alo CHUISt CH2t Brot Jp 
- 3X4 Loran/DGNSS S# AAW APS SOA OA! 
M&S Ga 7S A4e 
- 7} FOIE] SAI BW! SAAS 
Eurofix SAAl 2a Jie’ A 
AQJls] 2 gaye 
- FERNS 3] 2 Sa} go 


FERNS 


(Far East Radio Navigation Service) 


2=, clAlOt 49401 Loran-c / Chayka SASS 
4 (22 BA > BH2 Hye wa sz) 
° OH (OFC Ss} S]O} WHA 
* DGNSS / Loran-C 5 @ 
-Ssqa soaitsol a 
-SaS AHA SAS AA AY AMS 
-OS4A #an AAAS 
- 3A AMSA AAS 


aaa GS Thea 


Ol Mists AAS! on a 
HYASS OlSolI| AoA H&E 


—) 


- 8 We 
-575=, 174(S Ake 74, AA10H), 624 A('98 - 2000) 


: Sellte: AHS B So( GHA SA) 
: GPS, GLONASS, DGPS, Loran-C 


-AS, WAY ao SE, Fe BA, tNess, aH 


Om 
gt 


i 
J SSD 
Ou Ox 
& 
0 oF J2 


AS, BAT} AMA, AUS) wl S}(Qi et WBA) 


SE IWBFNUAM APS SS BONS LSA ASRCe BA 


ANG 


2! DGPS JI 


int 
N 
nN 
i 


ax 


Al dil 


— 
— 


OA Ol SB 7} 
= software 


Go 
Ga 


St At 
MOMAF 


a 


Ol 
= 


2A 


& RSA0N=E 


-AAl2t Ae BRANES FRISBadaA SAaerg, BCS 


(2001H FA AS As 


tO OS IS 


fe) 


[=] 
2s 


Internet 


homes 
[aes 


Ace] AS 


Ss 
+ 


ae OF 8 fo 


GPS 2 DGPSOll 2tokOo] SSI AFStO] USA 


CS Af SSurseudc. 


OH SF ATS? Ste HAI at 
02-3148-6340 
b6215h@momaf.go.kr 


—118- 


0] Al Yo] A] 2] Is-801& °]2% IDGPS 74 


IDGPS using IS-801 Standard in Mobile Network 


uo * = Paced 
. eee", 24e",— 


2 


Zl FFO1"** 0] a3 a 


+A] Stal A Az}S SF} (Tel:02-452-7407; Fax:3437-5235; E-mail: mphyun@kkucc.konkuk.ac.kr ) 
+421 Stal 21 2}-S 3} 5} (Tel:02-452-7407; Fax:3437-5235; E-mail: drhskim@kkucc.konkuk.ac.kr ) 
#442) SS} a A Zb-S 5} F}(Tel:02-452-7407; Fax:3437-5235; E-mail: gijee@kkucc.konkuk.ac.kr ) 


4+ 


ASU] Spal BY} S-9-F SS F}(Tel:02-450-3358; Fax:444-6670; E-mail: younglee@kkucc.konkuk.ac.kr ) 


Abstract : One of the methods to meet the positioning requirement for E-911 is to include GPS receiver in the mobile cellular 
phone. Lately, several position location solutions for E-911, so called network aided or server aided GPS are being developed and 
proposed. The standard for this is TIA/EIA/IS-801 by Telecommunications Industry Association (TIA) announced in October 1999. 
This standard defines the formats to transmit and receive the GPS and communication data between the mobile station and the base 
station. This doesn't define which positioning method should be used and how they should be implemented on the mobile 
communication network. There would be many different positioning methods and implementations satisfying the IS-801. The 
subject of this paper is to investigate the positioning performance of the IDGPS using IS-801 Standard in CDMA networks 


Keywords: E911, IDGPS, 1S-801, CDMA network 
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Abstract: 
DGPS correction message must be transmitted from the reference station to DGPS users. By now, Radio 
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this paper, we evaluated the performance of the Internet based DGPS by various experiments. 
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Using DGPS we can compute our position with meter-level accuracy. For DGPS positioning, 


beacon, FM Radio have been used for DGPS data-link system. But the prior data-link systems have some 
disadvantages. There is limit on the available distance and data transmission rate. Therefore, we suggest 
another DGPS data-link, the Internet Based DGPS data-link. Nowadays the GPS users can use Internet 
easily, so if the reference station is connected to the internet, they can be connected to the reference 
station without special equipment. And with the Internet based DGPS correction, multiple users can receive 
DGPS correction from single reference station regardless of their distance from the reference station. In 
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ABSTRACT 


Since 1994, KARI has been conducting the research on the Local Area Augmentation 
System for precision approach and landing of aircraft based upon the concepts of 
CNS/ATM (Communications, Navigation and Surveillance and Air Traffic Management), 
which was previously called the Future Air Navigation System (FANS) at the International 


Civil Aviation Organization (ICAO). 


The LAAS system includes GNSS Satellite Subsystem, 


LAAS Ground Subsystem (GS), and Airborne Subsystem (AS), and makes it possible to 
satisfy the precision requirements for approach and landing of aircraft, CAT-I to CAT-IIL 
In this year, KARI has tried the first flight test at Ulsan airport with the flight inspect 


aircraft of MOCT (Ministry of Construction and Transportation). 


We would present here 


test procedures, test results, and pictures taken during the test, and also our lastest version 
of LAAS subsystems-GS and AS would be shown. 
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Abstract: This research is concerned with the design of LQR/LQG controller for automatic ship 
maneuvering, which is in fact the automatic path-tracking control. We _ utilized the 
non-dimensionalized equations of motion to generalize the derivation of control law. Also, the 
equations are linearized to apply the control technique. The optimal path-tracking controller is then 
obtained by applying the LQR theory and the estimator is designed by the LQG control technique. 
The problems arising in applying the GPS measurement to the automatic ship maneuvering are 
discussed in detail. Numerical results show the effectiveness of the proposed control algorithm. 
Experiments will be carried out to verify the theoretical results. 


Keywords: Automatic Ship Maneuvering, LQR/LQG, GPS, Gyro 
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Abstract: 


The rapid growth of technology as applied to such areas as precise navigation/position, 


digital communications and space applications has resulted in needs for time and frequency 
synchronization over large areas at nanosecond or better level. Satellite technique may offer best 
chance for improved synchronization of time and frequency. 
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Synchronization of Remotely Located Clock 
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Abstract: 


Baa e)] B7) SH EHF (Tel:042-868-5146; E-mail:kimy @kriss.re.kr) 
Ay si} St F(T el:042-821-5665: E-mailidcpak @cuvic.cnu.ac.kr) 


In this paper, the clock synchronization between reference clock and remotely located 


clock has been carried out by using commercial GPS time transfer receiver. The remotely located 
clock was locked to the reference standard clock by GPS common-view technique. The results of 
clock synchronization experiment have shown that the accuracy of the synchronized rubidium clock 
can be kept within a few parts in 10’, which may be used in the network synchronization. 


Keywords: synchronization, common-view, rubidium, accuracy, GPS time transfer receiver 
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GPS Round-robin trip Calibration for APMP 


SAS, “OAR, “ol Sy 
Skt ESS =F 1 (Tel:(042)868-5147; Fax:868-5287; E-mail: shyang@kriss.rekr ) 
“shat Fv} Sh 4 = 91 (Tel:(042)868-5140; Fax:868-5287; E-mail: cblee@kriss.re.kr ) 
"Sh at SE vy} St =F 21 (Tel:(042)868-5142; Fax!868-5287; E-mail: hslee@kriss.re.kr ) 


Abstract: The purpose of this intercomparison is to compare the electronic delays between (and 
including) the antenna and the reference input of the GPS common view time transfer receivers used 
for the realization of UTC at several laboratories in the Asia-Pacific region. 
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Fig. 1. Clock comparison system using 10 kHz and 1 pps signal 
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Fig. 2 plot of the experiment setup. 
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Fig. 3 Results of comparisons between the GPSCV(NML) and 
those of the participating laboratories. 
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Table 2. Assumption that the internal 
delays of all participating GPSCV rx are 
identical (50 ns). 
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Fig. 4 Results of comparisons between the 
GPSCV(NML) and those of the participating 
laboratories.(internal delays : 50 ns) 
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[1] Report on the APMP Round-Robin GPS Common-View 


Time Transfer Receiver Intercomparison Experiment: 


Round 1 - October 1999 to May 2000. NRL 


Table 1. Results of comparisons between the NML travelling GPSCV receiver and those of the participating laboratories. 


aren Slope of line} RNS Paso 
(annentdey Start End Difference of Bi tnmee fat Beriat pen # of Tracks |Error in the 
MJD MJD REFGPS: N Mean 
host-travel ling (ps/d) o (ns) o /VN (ns) 
NML(Sept 1999) 12.3 -34+ 14 3.4 0.11 
TL 16.3 -13 +432 


NAO -32 £70 : 0.09 
CRL 30 £30 3.9 0.19 
NRLM 31430 0.19 
KRISS 30+ 29 0.15 
SCL 4.0 ~22+60 0.2! 
NML(June 2000) -27 + 16 0.18 
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Development of Position Location System using INMARSAT-C 
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Abstract : This paper deals with the implementation of position location system using INMARSAT-C. We propose a system design 
and software program for a real-time based position location system. The overall system designed in this paper consists of 
INMARSAT(International Maritime Satellite), GPS(Global Positioning System), LES(Land Earth Station), and Packet Data 
Communication System. The developed software system is composed of LES calling procedure, Position Data Decoding, and Data 
Communication program. The experimental result of the proposed system implemented in this paper is evaluated via real-time 


experiments, which show that the test results are satisfactory. 


Keywords : INMARSAT, GPS, Position Location, LES, GMDSS 
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2] Frame Xo] °]°] Packet FAH WE Dole] = Hinet-Pys 
Arey FAzq] AAMOE Ago BAA] aAlzdo] 
X.25S 2)23}2] YOU Packet2o] PAS Ia] AES] PADS A} 
Sarl, FAlel X.25 Packet S Gal ws] doles Basta 
HS7] AZM FS AAA ASA. o] MUHE wH7} 4 
BH EYEE] SARS WolY=S AA BES Dole] we, 
= WES TAF ASF YES Also Dole wes #e 


aq. 


HE HES #S3}e AVA ASAE SAlB Feo] 
E]S} C++, JAVA, HIMLS o] 28 HB AAP] Sule] z}a] Ze 
Agee S29 A) WE YES A B2hoq aAsaAyg = 
= Data Basel] 4A ZAC, Daw B aAAdo] Ags za] = 
ZBAGs ASS] DB] AAs Aolys A AWAoe Asa 
| ol SEA2] AAS Fa | Vey. 


3-2 HAH) Ze2Ig 


3A BYAly (Telnet Server) 
Fst] A] ols} PS Bgol 


INMARSAT Mobile Number 


Password : 


SEND 
T= AUF tc3 S) 


24 2. AA] ZBUG SEE 
Fig. 2 Flowchart of position location program 
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SANZAFOLE BED HF “Poll & Data Report” AH!2~ YF 
Zoe fl $s}, HAS yrs SE OFA] DWH Ez 
AE VAM. 471A OSA DH SAAD INMARSAT-C + 
ALT} Babee Muyo! BD, INMARSAT Al7]PHC} Holshe 9h 
2] 224-8 Ber}, Yao] ASA FAl(Send) FHS AAspo} 
INMARSAT AVS Fs] HGS FAs ac. olm Sal 
2 ABso AY 7] (Enter Key)S BAS WM, “lS ao} = 
WA7} & Ab] 2] So Atte INMARSAT FAOe Sad + MS 
Yop.” ye QD WAZ 7 ULE APs A 2s 87)aa E 
= A 3SA7AF oe ASA PHS aAlSafors etc. 


SA YSo] ASAOR FANS Bole Telnet Server} 
t<3H o]jo] cha] BUA SE Bes Fala anoles Fel 
acy, ol, -Alad mole} BAS 162%) (Hexa) HEH FY 5/9] 
Q7) We B aAdo) # AAG oA] Wat THUS A} 
R35 BADLS WHA Wc}. WS} A Raw Datac] PHE HS 
3} 329 35 Poul, AHA FUE FolejMol2zd Ags 7 
Uy AES hol A) AS ACH. 
Bit No. 
8 7 6 5 43 2 1 


PiCys Type 


. LESID 
Member ID 
ca pL | 


1 


Lat. 


Freee | Mines | 
[Fraction Pat] W] Dew | 
ae | 
a orcs 
ap oew 


-- Attribute -- 


Data Report Packet 


24 3. 4a) dole] a4 
Fig. 3 Position Data Format 


Byte 
oo nyt an pwn 


Position 


Long. 


4. J4)81 244s} 


ls Sale] 14S Helse AMS Sells) Aste] $4 
LES(Land Earth Station)olA] Faz) dlolej2] Aas wa ag 
3x} Ach, of7|Aj ES] OINMARSAT-CS ARB SHE Ae 
BES] SAS 71! Datal AZ SVC. 


B Ao] AVS AY] TrimbleAtS}[5] INMARSAT-C 44) & 
aa] as 7eusdaw Aga] 2415)¢o Ve AS WAZe 


SO} Mobile ID7} ¥12} Zo] 444076712, 424838210, 444050612 
S| 6a2] AYOCe APS HAct. 


21933} Zo] INMARSATS BLE] 4-418) Datacll4] 6 - 10 Byte 
o] ey Qe HA) datas 2hal MAS Programs o]-B sto 4] 
zZ] aa Sale Ho] HA) Datah AN, BE Sl] FA HLF 
Zo] seo] HA) eh. 


162? 


444076712 


#1. Ag sole 
Table. 1 Experimenta! Data 


135,12 


34,39 


flat__[UpdateTime 


2000-11-03 2S 5:01:43 


1628 424838210 121,09 26.12 2000-11-03 2S 5:01:43 
1629 444050612 118,03 38.55 2000-11-03 2S 6:01:43 
1630 435716510 = 131,39 32,26 §=2000-11-03 2S 5:01:43 
1631 451445810 120,00 31,58 2000-11-03 2S 5:01:43 
1632 444076612 118.53 24,23 2000-11-03 2S 5:01:43 


12] WS FS DE SEAS AY 3714S User IDO] a2 
Address INMARSAT-C 24]2] Mobile ID o]tk. © lon® QE 
(Longitude), lat 9 &(Latitude)S UE} Ac}. Update Time® 
Aay2l azo] zAQM ol SL ale] 4) Data Also] Deol 
AAR Be ALZES HALSEA] Sc}, H1z} Zo] INMARSAT-C 441717} 
Ala] S]o] QE S72] AB Alyto) Mae Sel SF AZ Mobile 
ID7} 44407671291 AVS] Hae BR 13512", AR 34.39%] 
31, Mobile ID7} 42483821091 M¥Bbo] MAE Be 121.09%, F 
© 26.12% $US Hel S FA UWRCH. 


os} BE Hal Moles ajolElujo|Ac] 4 Aste] ley y 
HBAS Sato MYM 7] aA z4EAoA| Bole 
WSFA AIS Aste HAE F Art. 


5, a= 


7\ 4 Zo] B EES GMDSS Heke] wet 300 EH ol ss] 
BE Apo] YAO Bays] AE INMARSAT-C +217] Sh oll 
EPH VRE PS FANE lsaste] FALe AAA dal 
UY AALS VSSU Sa] AIS Helse AGS TH 


HB Mate AADS BSA ANVUALS Falls) AAl 
B08 + YVOes, BA NS SAMA MA sy LS O} 
Soy te aa ast qesty GAAS Sa zy] AAD tas 
Ee ASAUN BHAA cola=S MAE MAA) zt 24 oO] 
SEMAMES SE, 1 UAE Use AZ Seq 
yy] Qset Mae SF Qa sale. 

res HB Uo] Ha) Sel AZGS 


1%) 


Ya Ea 
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Bugg) qa asaya aa 
A Design of Car Navigation System for Incidents and Events 


4) 3] zt 
Ax qty =EAASS4# (Tel: 031-299-0873; Fax: 031- 298- 2737; E-mail: hkchoi@hyupsung.ac.kr) 


Abstract: Car Navigation System(CNS) of ATIS(Advanced Traveler Information Systems) which is a branch of 
Intelligent Transport Systems(ITS) may provides drivers with real-time traffic information. But traffic congestion is 
no longer limited to a peak period and particular roads. Recurring congestion is predictable and occurs in locations 
where the traffic volumes routinely exceed capacity. Nonrecurring congestions is very unpredictable and is caused 
by incidents and events. The object of this study is to analyze incidents and events as CNS(Car Navigation 
System) provides the shortest path information. This approach could be applied to other CNS technologies using 
the Global Positioning System(GPS). 


Keywords: Car Navigation System, Advanced Traveler Information Systems, Real-time information, Nonrecurring 


congestion, Incidents and Events. 


A= : é 
LA = AZS tse BAe] CNS] SY24Aincidents and events) 
ACSAS SHs7] Wet ALAS AAseED BETS 34 

FSG DEA S79 SHGE zS4} BH STS] SAS 2 o] gtk, © 7)=5] CNSS Zils Ase + ME BVIYS 

so] 2S4e1 ERS] BSYE BFAD quer wea a S Page BRS Aetsar. SY Iso] Be Barz 

Ae 22 AST Wd. ASIE 2340 SAIME AEA Incident Duration? JS - Ua H22 CNS WY Ye 

3) S249 Bey} S$QATD YOO, oat EAS Waet7) = BAS = YWHeE AS AAlswz Hr. 

Ae eee 2s% mee (ITS: Intelligent Transport 


g a 
Systems)°] HAISS SALE BPA ASNT Nth. ITSE 2. S43 Bz] AA 
ATE 7S EEG AABYS AAS ELMSSES 4 
Ase7] Aa Atsst2] BH] sh 
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ash SMa eae Se Sas Sel, ae 20 % S38] Segal PRA Be Ra 
Ee ERM Ao] BA7IESS 48121 ITSE Sarda HB AT)S TEAS AVS Yo SQM Qelo] Fa gy 
BoA F AA Be AS UERRS Asta EAI y ole eye jd fe Sst AAGqA a wHEso} 
BIS SOD SUA ay is Yass 24 F wagye dao TEEA QdgHe Be] ySo] oq} B27} 
+ 27] 88 SS WAS +A SA ITSY SHE vad yao sees ASE] SHE AypaEe wv 
oFe) sean e dBA) WATS Aavenc ed Traveler information de 23) vee Ryo Seyseo wees qa aa 
Systems)2 WEARS FAt4 AEA F SAAN MYA 2a oeryy EL SIS Bas qaEzs. Pa Soe 
Se Hofolch. oS APst7) HAA LaAAcda PSgQus Al ls SYR Hao] das $2 age gq, ase evry 
Se - Qe aes ya A~Aa(CNS: Car Navigation Aen 7) HAqS Msqy Soiady ws) Faq(SoIa qapa, Fy 
ce 7A dest Beet. CONSE He 77 eel RY oo} Baa eaah, EE Qalgal H)} DSRU Saas 
a Ael7] Wel aa7lee Fad Aaa soe SUS 2] 4}(traffic incident operator/manager) 22 %o] UAT 2 
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Bee Wes AUD eh. G Bo) Strode EEA wasta ya, seg 
Al z}4 2] 2] =(DRM: Digital Road Map), 973 & 4A) 2 S(GPS: ol AS Yaa Sa Be TEMA saAo wwsjojok a 
Global Positioning System) $2 E¥sto4 Bz] Bal $4 th, Baas ByAbashi} ACF SS ol gst Asawa 4 
HEA} BAZ AM (real-time information S Ass} #7} 2) 24 HsA SVYS Yoo Weve eqzq adaa, Sw 
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Normal Traffic Flow 


Incident Verification/ 
Identification 


Incident Response 
Incident Clearance 
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Fig. 1. The Concept of Incident Management 
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3.2. HAZ] 27 Aq al 

MEA Boke) uAlz] 2B Bale BE-YEMA, Ve 
Na-Aq, YEMA-oME, WE-BA, ae-Aal-ag, ae- 
alo], 4, UEMa-ME Az) 2H) & TRAE MLA 
ate GS. ITS ek BS FURS MAA MAA 
1s, MA) Bet, Malzte) a 7e4) He each. APY ASS 
UA}z] ol at, HSS ASe wMalze aAgslcp. cela uA) 
Zl¥- ove] BB (Data element)#} 4422) Betahwa, Ala 
QYee Testa) cp. 4%, Zhzb2] eAlazl7} Be] Srp] wr) 
HY QBRss], SHAS) Babee) 7)SE& Gepults Blo]rt. 
ret BMP Set vaya 57tzye] vAlal 4 ets} 24 
FEAL S] Al z)e AAVSt +e Qh. 

YEMa-oME Wz] 2HS BBY Wes LH 
wAlSH(roadway incident), 5. Aplanned roadway event) 
Sy} es AMS AS}. AASB Re] Sal ayes 
au, Abal HE o]MeEs} Wales cele SHS} Se] ula 


AG Meshal eich BWI val 2] Bara ais —24) 2st B 


Oo 


BD, YEMNa-ale Wz) 22 5749 WA APLE TF 
Add, yesa- oe WAZ] 2F2 aq-YeEda-Sad 
3, AS-ES-WE, ME-AI-8E, YEMS-Seda- 


2@ewovs 


a, BE-OME-WA, S SMS WAZ] HOR Tyact. 


19 2. Sw AS} BA] 2] BF a 


Fig. 2. Classification Scheme of Incident Message 


Message Set? BE 7]8t JAMS] AA 2 Ls Ara 
rod 3} 0] Tn, A) 22 Zt EAA) AS TIRyE 7|2 421 AWS 
E(VRAOe WAZe}y Bs Per}. S45], s+] Message 
Sette SYA Fas WHA SH, o] GF AAA ABS AG 
a7] Bato Aw] Aes](Semantics)2+ BAl(Syntax)S Ales} 
A 2th. 4, HZ) Ashe] Bo] sol] BAe] Alle Message 
Set AISA, WA Weezy AAgts] FS LIVE oF 
Ye DEAqNAags Ge ITS +4 29a WS BU Se 
Ate] AR TBS Ma] BRA Wf Saet LE F 3B 
ITSS Ht HZHA Message Set2] ASL AAG 4 
22 FE S457 Asa Ass ANAS A 
Rs Ao]rt. 


3.3 Abstract Syntax Notation One(ASN.1) 

Suggeadag 2 agg Base SYgs WAS 
23} Go] ASN.1 BAZ BAA; Aq-vEMa-Sus 
3 3A). YEMS-ME DE} 57] WAZ) AH, 24719) a 
AA= GSH So] AMAZE S45 BAEZ FE Ags 
FU. HF SLUSH 4] WAZ) Petol= Global message 
= AAAA WANA} SHS] QT, ITS SBF 
# OS 2Gs7] At Fe] MAPA7} BRsee 
1S 48 gee) WAZ) SASZ FYS LRMS7} VSS oO} oF 
Seq, °c] ZZaUS LRMS Dole] HSLLY ZBEBS B 
22 AS 7AQHhD AB. SUIsaqay Ze ayy w 
Bese SH9S ANAS 423} 2°] ASN] GAS ZAG 
(2). SSIS 92) WA) Bol = Global messagel*Al As} 
= 9A)H3A BAALWLRMS: Location Referencing Message 
Standard)7} ZMS]o] Beth. ITS SRHoks] eYWMs os 
AMs7]) At Fe] Azlpe7} BAges ofS He 4g] 
WAZA] BASH FY LRMS7} BAAS BEG, eo] ZZR} 
V2 LRMS dole] FSSBY ZHESS FRA AS 79S 
wz ASH. 


MESSAGE GROUP : NETWORK-EVENT 
MESSAGE SET : CURRENT~-NETWORK-INCIDENT 


MESSAGE : INCIDENT-IDENTITY 
INCIDENTIDENTITY DEFINITIONS AUTOMATIC TAGS 
= BEGIN 
IncidentIdentity 
{ 
orgn-ContactOrganizationIdNumber 

Integer (SIZE(0..4294967295)), 


n= SEQUENCE 


trnt-IdNumberNumber 
UTF8String (SIZE(1..256)), 
orgn-ContactOrganizationNameText 
UTF8String 
evnt-DescriptionT ypelncidentCode 


(SIZE(1..128)), 


Enumerated 
{ 
Other No Additional Information Required (0), 
Other Additional Information Required (1), 


StalledVehicle : (2), 
Vehicle Fire (3), 
Roadway Debris (4), 
HAZMAT Spill (5), 
Vehicular Accident (6), 
Transit Related Incident (7), 
Weather Related Incident (8), 
Earthquake (9), 
Landslide (10), 
Flood (11), 
Tornado (12), 
Hurricane/Typhoon (13), 
Unplanned demonstration (14), 


}, 
evnt-DescriptionT ypeIncidentOther 
UTF8String (SIZE(1..256)) OPTIONAL, 
Integer (0..4294967295), 
ENUMERATED 


evnt-IdentificationNumber 
evnt-IncidentStatusCode 
{ 
Other No Additional Information Required (0), 
Other Additional Information Required (1), 


Incident Detected (2), 
Confirmed and Responding (3), 
Cleared and Recovering (4), 
Over and Done (5), 
Other (6) 


}, 
evnt-IncidentStatusOther 


UTFS8String (SIZE(1..256)) OPTIONAL, 
orgn-ContactPhoneLandlineNumberText 
UTF8String (SIZE(1..32)), 
evnt-UpdateUtc 
UTF8String (HHMMSS:ssss), 


evnt-UpdateTy peCode ENUMERATED 


{ 
Other No Additional Information Required (0), 


se 9 does 


Other Additional Information Required (1), 


= 


“ 


2 4c 


New (2), 
Update (3), 
Delete (4), 
Clear or Closed (5) 
}, 
evnt-UpdateTypeOther 
UTF8String (SIZE(1..256)) OPTIONAL 
} 
END 
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Network Time Server 2.4 22 #4) 


An Analysis of Error Factors in Network Time server 


AT SI at ie ay Bas Tel 640-3179; Fax: 640- 3038; E- mail: ere sites ac.kr > 
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Abstract: Network time server performs the synchronization to universal standard time and provides 
the time distribution service through the Internet. This paper introduces the essential technologies 
and each inherent error factors and presents a mechanism of time generation in a network time 
server. From that point, we obtain the periodic timer interrupt and interrupt from the external 
reference clock source have some important influence on the time anomalies and then, aperiodic 
network interrupt does also. Following the illustration of the test environment, the measurement 
analysis of interrupt processing delay, interrelation among interrupts, and influence of interrupt 
service routine with long execution time in a network time server are presented. The analysis deals 
with interrupt characteristics of embedded systems based on aspect of time and can provide the 
practical consideration points when designers who develop the embedded systems. 


Keywords: time synchronization, network time server, interrupts, RTOS, GPS 
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Abstract: In network protocols and applications, like distributed system 
synchronization or network performance test, network delay time is one of the most 
important parameters for correct operation. But TCP/UDP/IP does not provide these kind 
of requirements whereas almost and more and more network applications are going to be 
based on Internet. So common methods like PING or NTP(network time protocol) use 
“roundtrip delay time" to get network delay time. But "Network State" is not constant 
and unpredictable under TCP/UDP/IP. If we don't have additional aid OTT(one-way 
transfer time) has no meaning. OTT includes clock offset, clock skew, and network delay 
time. We will show that to get the network jitter and.clock skew gives more stable and 
accurate clocks. In experiments, we could get network delay time roughly by using NTP 
and the range of network delay time was more wider than general NIP accuracy. 

In following section II, we will analyze time error factors in NIP especially network 
delay errors, and section III, we show experiment results of network delay time, 
section IV, effects of network delay error correction will be showed, and follow 
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Abstract: Airborne laser mapping system, developed in 1980's, is considered as one of the major 
state-of-the-art technology for generating DEM with high density and accuracy. This has been done 
by technical researches and innovations in hardware performance, positioning accuracy, data 
processing and application. In this paper we have studied overall status of airborne laser mapping 
technology, application and feasibility of domestic use. 


Keywords : Airborne Laser Mapping System, DEM, LIDAR 


1. A= jo me, BE ATINE SSqAsso] Bet aAw7yel 
7S BR BSaAS AE, NE FIBBAS AA7I1E9 WHE 
1A7] JAABASAINE coke AAgqayS 4 27az 4, SsalASsss] SUEY SGI GES TVsAct. 
Be GaAty7 SBS, TEMA Het TESA BSA 
22 2 UIs Aq 2 FAyo] SZ SAAD Qe. 
APAAVES AAQBARS} 7]RFO] HE AVSttyws| aA 2. Lidar 7] 2 
Aone] FSAI TSS Ao] Ba CaA2d7] S71 g ual 
A TSAS) 2S] BARS Sse At. H7 ASAAQR 2.1 FSHMUAMGF 72 
f LH5o VE FIEDABE AVE WEA) BLsyoy 


SEAS 2 tz, ES, 24, Aa, S4, TA S 4S Lok BS qASeys AAWS Sa z-EA a7} Wests Aqsa 

yA) 380) FD VM. ATA) TARUARE APSee] 4 7} He SSA AAES BH FHS 32] Bar Sed 

a ue, 727) =o Be Hy, AMSsto] Bs HA SOB FM. Ht AA7ISASSs SH] YAS Bet, 4 S33} 

AFARQOA, ATE aso] de At Faewtse Ase SF aAgzbo yy Pawe zs, EA ASIANS] FABUVSAY A 

BARS Ho} REO} Mr}, Ao] fz] st ue Qo, a3 Aszo YAR +AgAge 
gJ2 H2e Saaseeze veo] ad, FUSES SA,  loem, FBYSE B0emS BAHE FAFA S%o]LII[71 

ASH, VA AWS ASolFe Fe OBS ASD Aq. 

244 FAADSBee oS FARDARS AEASL AY 211 FSF HUAFS PI] Ae] 

ESS St Sao] PRR ASE} B27) Ss] low () tsaoAsese] az 

texture2]4, HAlz}9} Zo] 7HUA7 AD 287A Yes BS qol Ase alzde2 Ye4°# LIDAR(Light Detection 

Bol we zASaqyjye ASE ArAsA GAM. and Ranging) A}2"o]z}a2 Bz}74Ub ALTM(Airborne Laser 


At 200 Set TF ARDS BP] BAAS Bolsa ur yA Terrain Mapping) *)24o]#}31 #2]>), Laser scanner, mae 
3e +TARUASS zASosR AAAS AES Wes7) INS 2 478 HW, S/We 745 VHF 1. 


ya Be esAo] Boon, a WHA Ysxt7} LIDAR(LIght AAAAS 772 Aae AolA 7)7/9} & ae co 
Detection And Ranging) 7 2@2] 7eolc}, LIDAR A2AWS °]t}, GPSS} INSE FH BYRSS Bep7s= vlay7)o) ypskz} 24) 
Airborne Laser Scanning J28oc}TZE Hey, 2 Ato SASH. 42] waAASS GS ABYSS] BAAN 422) 
AE NAGS “8S a]AS B77)", JAS A882 ASASH Se ZAstA) Bch. 

AS “HSao|ASeW olehe BZ ASa7Z Br. Beal 

ASSES GPSS 7S] EMS 1980dq FUE BAW! (2) Sze] 2] so} 2 A} 

Ao] ALAS|O] 199001] Sto] 22S} SB aso) SAS Art. YSAqoAzeo] Fao] VIS FE QABE Aah, oly 
YRS BP as SA] OF Buss SAQHES] St, zB FH] 9} GPS B INSE Bis = qoaqrAzd, Sa, det 
AQ PRS SS At YEA SALTZ Ms AW TWae oq} 2254, 2m AtoZ BA Ao] Bee 2aSRe 
/DABEY FARWAS YS Beal4Ase 7)Zo] BE + U1]. 

Ss] tt. Lae Alas $AOAZEZ FHS F YOd, ds Bu 4 


eayepoy ae Calida S7}2le]gwaal teat, S 
ae F415, FAS S GE 1G eau TH sI94 
4] Rofo BSS. Mech. er} Eo} ase FAIS 
$ aed #9 EUAES Asses Saw AM, 8] Bt 


QA Wola B2)7)2] BW (bias) el S(gain), Ae] 4 B22 
ala, INS?) £4) °] G(misalignment)*} 2°] = W#(drift), GPS 
A924, d7]aare] GPS Aad, za BAS, AQ, Malate 
} Sho 2 Byatd. -Aaozs= BA BAL, pointing 


5 34a FURNES AINE Mo Baan gee jitter, INS, GPS, AV] AIS, Wise FZ GIS WE 


,es 


aaa BarsAy, 44 cis #87 Sage 


= 
re 
a 
a 
Ali a 
a FA)EDABY FRE YH AAD Ath 


=160> 


YEYo|ASBo] Mel 


29 1. 
Fig. 1 Airborne Laser Mapping System 
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A study on the crustal deformation between the southern Korean Peninsula 
and Kyushu Island by using GPS . 
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Abstract: 
mantle diapir, 
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To confirm the crustal movement of the marginal sea between the Korean Peninsula and Kyushu Island by 
Data for the period January 1999 to May 1999 were analyzed, with data from International GPS Service for 


Geodynamics(IGS) global sites, to estimate daily coordinates. To fix the estimated coordinates to the terrestrial reference 


frame, the reference site is the Suwon IGS site in Korea. we find that site velocities are (Vew, Vns) = 
(0.50, 0.30) em/yr at Hagi2, (Vew, Vns) = 
~0.22) cm/yr at Meabaru,(Vew, Vns) = (-0.92, -1.10) cm/yr at Nagasaki and (Vew, Vns) = 


cheju, (Vew, Vns) = 


(3.03, -1.44)cm/yr at 
(-0.42, -0.38) cm/yr at Kamitsushima, (Vew, Vns) = (-.28, 


(-0.24, 0.01) cm/yr at Tamanoura. 


Key word : diapir, volcanic activities, MT(magneto-telluric), precise, orbits, displacement velocity vector. 
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Fig. 1. The map of tectonic and GPS stations 
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Fig. 2. Flow diagram for GIPSY-OASIS (after J.S. 


Hwang, 1999) 
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Fig. 3. Time series of horizontal coordinates at selected GPS 
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Fig. 4. Displacement velocity vectors of GPS sites with 
respect to Suwon IGS station estimated from linear 
regression of time series for the period Jan. 1999 to 
May 1999. 

1a 4 #9 IGS BFAq We 4 ASsa9 Yq wAse wy 
e}] (1999. 1.1 - 5.31) 


Site Lat. Long. Vew Vns 
(Periods) (degN) | (degE) | (cm/yr) (cm/yr) 
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+ = + 
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(1999.1.1- 5.31) | 7” oyrare cae 
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Nagasaki ¥ x : fa 
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Tamanoura 
-0.24+ + 
aaetason 128.619 | -0.24+0.04| 0.01+0.02 


#1. 24 GPS #429] BY SEME, BE YL BEBA 
Table. 1 Site velocities at the GPS sites, Values in the 
parentheses are S.D. 


22 43} ¥ 1A ME HS} Go] zABaA/B4 ag, 
HA2A*Y (Vew, Vns) = (3.03, -1.44) cm/yr, Hagi2 4S 
(Vew, Vns) = (0.50, 0.30) cm/yr, Tsushima #42%(Vew, 
Vns) = (-0.42, -0.38) cm/yr, Meabaru #4A"(Vew, Vns) = 
(-0.28 ,-0.22) cm/yr, Nagasaki #22(Vew, Vns) = (-0.92, 
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en/yr2}] 42}2] HE SE GS BaF D Wrh. 

ayy 9} AMSe dA AEBAeal Ye Axols 
2, AF2449q He 724 ASS GPSAASPE ESs}7] 
AIAV=E YOR Be ASH B79 ARS Aq se 4 
¢ $2 NaHS AST US ALE Bead. 


fy & 


ba Te 


S-l. AZAA/SAAS, ATBs+E(Vew, Vns) = (3.03, 
-1.44) em/yr, Hagi2 42% (Vew, Vns) = (0.50 ,030) cm/yr, 
Tsushima W#42AE (Vew, Vns) = (-0.42, -0.38) cm/yr, 
Maebaru 242% (Vew, Vns) = (-0.28, ~-0.22) cm/yr, 
Nagasaki #422 (Vew, Vns) = (-0.92 -1.10) cm/yr, 
Tamanoura #22 (Vew, Vns) = (-0.24, 0.01)29] 229] 4 
GEIS HAST Ve. oe 47/4el B44, BAl elo] Set 

Aoe A7)Ad ABABS S Bl FH 
ye BO Swe] zjz2] 4d qe zes Yt $e Ha 
x 


[1] Karasawa M. et al,. Japan-Korea GPS continuous 
measurements for studies of plate motions around Japan, 
S4th meeting program and abtracts, The geodetic society 
of Japan, 1996, pp.11-12 

{2] Mogi, T., Nakada, M., Ikoma, Y., Min, K. D., Chung, S. 
H., Shimoizumi, M. and Yumoto, K., Mt survey in Jeju 
Island, Korea, Autumn Meeting Soc. Gemag. and Earth, 
Planetary and Space Sci, Sapporo, 1997. 

(3] Teruyuki Kato et al. Initial results from WING, the 
continuous GPS network in the western pacific area, 
Geophysical research letters, Volume 25, No.3, 1998, 
pp.369~372. 

[4] Thierry Gregorius, GIPSY-OASIS I, How it works, 
NASA + JPL, 1996. ; 

[5] Zhu, W. Cheng, Z., Wang, X, and Jiang, G., 1998, GPS 
measurement of crustal motion in Chinese mainland and 
its surrounding areas, unpublished manuscript, 1998. 

BR, Ea, $S5, URE, Teg Atas gale 

@ GPS network? (1), IES PSP HERS 2B AG 

Sth oW TT), 1997, pp.365-368. ‘ 

[7] aS, 8S, 2+, GES, GPS Yet AAAs Ba 


[6] 


-173- 


[8] 


[9] 


AA, 9 EFYSPHYWSS 2 PITH), 1997, 


pp.369-372. 

22%, Teruyuki Kato, GPS= °]-2% 
Awad sz apuadza fe, 
pp.43-48. 

Salat, TRF BAM BAA GPS4Ao] 
Of pail AAP SP 9] =, 1999, p65. 


Sea z YEAS 2] 
A121 ALS, 1997, 


ge at, gate 


~174- 


GPS4 TS443E9 STAad Be AT 


Combined Adjustment of GPS and TS data 


“Add sal 22) =F S3}(Tel:053-850-7313; E-mail:yjlee@bear.kyungil.ac.kr) 
“Add sa d]et 4 (Tel:053-850-7319; E-mail:geo5055@hanmail.net) 
“cH 3} 2+] 3] (Tel:02-679-4042; E-mail:ehjung @kasm.or.kr) 


Abstract: 
GPS vector data. The initial results of the adjustments show that the three-dimensional adjustment 


The purpose of this paper is to adjust the combined networks of TS(Total Station) data and 


techniques are powerful tools in computation of national network based on existing Bessel(Mean Sea Level) 
ellipsoid and/or new geocentric coordinates. Also they are very useful in accuracy and reliability analysis. 


Keyword: GPS Vector, TS data, Combined Adjustment, Observation Equations, Coordinate Systems 


TEAS Ca PAW FS, SS, AVVSBs} Bose 
né, u Se 4S 4 (Deal Va ATSAasAdsgease Za 
gt(24 1 #4). 
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GPS(Global Positioning System)o] 2]3% 34:99 24]}aaqe 24 
AZ HBAS GYM AH] esa Garo] GFYoU 
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(Geocentric Coordinate)?] WGS84E} 21 a} 2] 324-21 2} AJo] 2} 3} 
oo X,Y,Z7F SAL ole} Ary. 

He, GPSS TSH S VB SH4Ce 24S AA 
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yA, FH#S3RR] AP By det Agee ux Pi a 
2} 4A (ellipsoidal normal)#} 2 4] x}, 
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ag 4= a TSYS SGU AOS FS 3272] 44 
2 FYE] ar 7JAME] 2071 SA, Aa 6144, 847 
SAO Z ol] Fol Ar. 
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A(TG336)S DAGL=S sAc4(6,8]. 
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Abstract: The Old-Primary-Triangulation Networks are established by National Land Survey in 


year 1910's and they are nowadays used as geodetic reference system in Korea. In this paper, the - 
estimation of the triangulation Network, GPS observations data of 60 points are re-processed by 
commercial software and the transformation parameters are computed by without height model. The 


results show that their accuracies are suitable for transforming a middle scale map to geocentric 


datum. 


Keywords:_ Old Triangulation, GPS Surveying, ITRF, RINEX, 7-parameter transformation. 
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A Study on the Pass Adjustment of LIDAR Data 
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Abstract: The vertical accuracy of LIDAR is quite good, but the horizontal accuracy is somewhat 
poor compared with ground survey or photogrammetry. Especially, the drift among observations of 
each flight path is one of the major sources of the positional errors. In this study we have proposed 


the pass adjustment method for adjusting the drift, and have applied it to the pilot project. 


Keywords : LIDAR, Pass Adjustment, Positional Accuracy 
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Abstract 


Resolving the GPS carrier-phase ambiguities has been a continuing challenge for sub-centimeter-level high- 
precision GPS positioning. Once the integer ambiguities are fixed correctly, the carrier-phase observations are 
conceptually turned into sub-centimeter-level high-precision range measurements making it possible to attain 
high-precision positioning solutions. Consequently, this topic has been a rich source of GPS-research over the 
last decade. 

A brief review of the previous work on ambiguity resolution and validation which has been carried out by 
many research groups from all over the world is presented in this paper. For a general understanding of these 
contributions, we have classified the ambiguity resolution techniques in terms of their characteristics. Current 
research trends and issues in ambiguity resolution and validation are described and a bibliography of 
representative papers is provided. ; 


Keywords: ambiguity resolution and validation, classification of ambiguity resolution techniques, research 
trends and issues in ambiguity resolution and validation. 


1. Introduction 


For the last decade or so, it has been a 
continuing challenge to resolve the GPS carrier- 
phase integer ambiguities. From ‘the first 
demonstration of the use of the GPS carrier-phase 
observations for sub-centimeter-level precision 
positioning and = surveying, to the latest 
development of the instantaneous ambiguity 
resolution on-the-fly (OTF) techniques, a lot of 
ambiguity resolution techniques have been 
proposed by many groups from all over the world. 
According to a recent report of the Special Study 
Group (SSG) 1.157 of the International Association 
of Geodesy (IAG), 300 of the 325 catalogued 
papers on ambiguity resolution and validation have 
been published during the last decade. Furthermore, 
research in this area tended to increase recently 
(Table 1). Why is interest in ambiguity resolution 
and validation growing? The answer may be found 
in the interesting scientific aspects of the problem 
such as integer least-squares estimation theory, 
ambiguity search algorithms, and the extensive 


number of scientific applications which can benefit 
by the highest-possible GPS positioning accuracies. 
Also, there are large commercial interests in the 
problem as well. 


Table 1. Number of papers on ambiguity resolution 
and validation 


Years | _No. of papers 


325 


The GPS carrier-phase ambiguity represents the 
arbitrary counter setting (an integer value) of the 
carrier-phase cycle tracking register at the start of 
observations of a satellite (phase lock), which 
biases all measurements in an unbroken sequence 
of that satellite’s carrier-phase observations. Once 
the integer ambiguities are fixed correctly, the 
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carrier-phase observations are conceptually tumed 
into millimeter-level high-precision _ range 
measurements and hence it is possible to attain sub- 
centimeter-level positioning solutions. However, 
fixing the integer ambiguities is a non-trivial 
problem, especially if we aim at computational 
efficiency and high performance (or success rate). 
Therefore, this topic has been a rich source of GPS- 
research over the last decade. 


In general, algorithms for solving the integer 
ambiguities have been developed for two different 
applications. The first group of the algorithms has 
been developed for applications using multi- 
reference stations in static mode. Multi-reference 
stations are occupied for several hours or even 
several days. Inter-station distance can reach 
thousands of kilometres. Whereas the second group 
of the algorithms has been developed for rapid- 
static, kinematic and navigation applications. Only 
two stations are usually involved with at least one 
station moving. The maximum distance between 
the stations is a few tens of kilometres. Time of 
occupation is the order of seconds to minutes or 
even instantaneous. However, conceptually there 
are no differences between the two applications, 
and the research directed to one application can 
benefit from research conducted for the other. 
Moreover, particular interests such as real-time 
long-baseline kinematic applications integrate the 
two approaches. 


For a general understanding of this research 
area, it may be useful to take a look at previous 
work in the field and various research projects in 
progress. One good source of information where 
we can get an overview of research trends and 
issues on ambiguity resolution is the IAG Web site, 
particularly the reports of the SSGs which belong 
to Section I. IAG Section I is concerned with the 
scientific aspects of measurements and the analysis 
of regional and global geodetic networks as well as 
satellite, inertial, kinematic and marine positioning. 
One of the main activities of an SSG is the 
international coordination of the ongoing research 
in its field. The following SSGs are of particular 
interest: 


1) SSG 1.154 on “Quality issues in real time 
GPS positioning”, 

2) SSG 1.157 on “GPS ambiguity resolution 
and validation”, 

3) SSG 1.179 on “Wide area modeling for 
precise satellite positioning”. 


This paper attempts to review briefly the main 
issues and research trends in ambiguity resolution 


and validation. First, we will summarize the 
previous work which has been carried out by many 
research groups from all over the world. Then, we 
will take a look at recent research activity in the 
area. 


2. Ambiguity Resolution Techniques 


Comparison of the ambiguity resolution 
techniques is not easy and not always feasible. The 
terms of reference of SSG 1.157 clarifies this point: 


“Despite the large effort spent by many groups 

from all over the world in devising various 
schemes, knowledge about their theoretical 
foundation, and how the schemes are related to 
each other, is still lacking. Different 
terminology is used and comparisons between 
methods are rare. Due to a lack of knowledge 
about the various methods, the 
implementations used in the comparisons (if 
made at all) are not always complete, thereby 
making the test results unreliable. Moreover, 
results reported of one particular method, are 
often difficult to relate to the results of another 
method, due to lacking knowledge of the 
characteristics of the data and the type of 
computer that was used.” 


As a result, comparison of the ambiguity 
resolution techniques in terms of computational 
efficiency and performance is not always reliable. 
To avoid this pitfall, we will restrict our interest to 
categorizing the techniques in terms of individual 
characteristics. 


Classification based on measurement information 
utilization 


So far, there is only one paper which has 
presented a comprehensive study of the 
classification of the ambiguity — resolution 
techniques [Hatch and Euler, 1994]. In the paper, 
many ambiguity resolution techniques are 
classified according to how they attempt to make 
use of the information contained within the receiver 
measurements. Although some papers briefly 
describe (usually in an introduction) their own 
classification, there is no significant difference 
from that of Hatch and Euler. The classification is 
three-fold: 


1) Ambiguity resolution in the measurement 
domain, 
2) Search technique in the coordinate domain, 
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3) Search technique in the ambiguity domain. 


The first class is the simplest ambiguity 
resolution technique which uses C/A or P-code 
pseudoranges directly to determine the ambiguities 
of the corresponding carrier-phase observations. 
Since the precision of raw C/A or even P-code 
pseudoranges is, however, generally not good 
enough to determine the integer ambiguities, inter- 
frequency linear combinations of the L1 and L2 
observations are used with a smoothing process for 
the estimated ambiguities. There are a few papers 
which give comprehensive studies of the inter- 
frequency linear combinations (see, for example, 
Cocard and Geiger [1992], Collins [1999)). 


The second class of algorithms includes the very 
first ambiguity resolution technique developed, 
namely the Ambiguity Function Method (AFM) 
[Counselman and Gourevitch, 1981; Remondi, 
1984]. This technique uses only the fractional value 
of the instantaneous carrier-phase measurement and 
hence the ambiguity function values are not 
affected by the whole-cycle change of the carrier 
phase or by cycle slips. Despite significant 
improvement of the original algorithm by Han and 
Rizos [1996], the technique provides relatively poor 


LAMBDA 


Null Space Martin-Neira 
Chen and 
BASE Lachapelle 


Kim and Langley 


Independent 


computational efficiency and consequently it is of 
little import other than of historical interest. 


The third class comprises the most abundant 
group of techniques which are based on the theory 
of integer least-squares [Teunissen, 1993]. 
Parameter estimation under the theory is carried out 
in three steps — the float solution, the integer 
ambiguity estimation, and the fixed solution. Each 
technique makes use of the variance-covariance 
matrix obtained at the float solution step and 
employs different ambiguity search processes at the 
integer ambiguity estimation step. The following 
are some representative techniques in the class: the 
Least-Squares Ambiguity Search Technique 
(LSAST) [Hatch, 1990]; the Fast Ambiguity 
Resolution Approach (FARA) [Frei and Beutler, 
1990]; the modified Cholesky decomposition 
method [Euler and Landau, 1992]; the Least- 
Squares AMBiguity Decorrelation Adjustment 
(LAMBDA) [Teunissen, 1994]; the null space 
method [Martin-Neira et al., 1995]; the Fast 
Ambiguity Search Filter (FASF) [Chen and 
Lachapelle, 1995]; and the Optimal Method for 
Estimating GPS Ambiguities (OMEGA) [Kim and 
Langley, 1999a]. 


Table 2. Characteristics of ambiguity resolution techniques. 


Conditional 


Transformation 
/Conditional 


All 


Single/Multi- 
Epoch 


Transformation 
/Conditional 


* Refer to the section “Classification based on methodology” for details. 


Classification based on methodology 


The ambiguity resolution techniques which 
belong to the first two classes delineated in the 
previous section are so straightforward that further 
discussion is not needed. As for the third class, 
however, an overarching point of view may be 
helpful in understanding the various techniques 
because each technique has more or less different 


characteristics. Therefore, in the remainder of this 
paper we will restrict our discussion to the 
techniques in the third class. 


Before we investigate the ambiguity resolution 
techniques further, we will consider the general 
goals of the techniques and discuss both the 
common and disparate aspects of the techniques. 
As can be seen in many papers, the general goals 
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that the ambiguity resolution techniques try to 
achieve are to increase the performance of the 
integer ambiguity estimation and to improve the 
computational efficiency of the ambiguity search 
process. The performance of the integer ambiguity 
estimation, which means the capability to 
discriminate a correct ambiguity set from all 
candidate sets, can be generally increased by means 
of more realistic stochastic models for receiver 
system noise (or observation noise) and well- 
defined validation procedures for the solutions. 
Unfortunately, it is not easy to distinguish one 
technique from the others based on_ their 
performance because almost all the techniques have 
common aspects as far as performance is 
concerned. In fact, it is the computational 
efficiency of the ambiguity search process rather 
than the performance that distinguishes the 
different ambiguity resolution techniques. 


Another common feature of the ambiguity 
resolution techniques can be found in the objective 
functions which are derived from the integer least- 
squares problem. Many ambiguity resolution 
techniques which perform a search in the ambiguity 
domain are basically based on the minimization of 
the quadratic form of the residuals. Principally, this 
minimization problem is referred to as an integer 
least-squares problem due to the integer-constraint 
for the ambiguity parameters. The same objective 
function, which is related to the minimization of 
the quadratic form of the residuals under the integer 
least-squares problem, can be derived using 
different approaches [Euler and Landau, 1992; 
Teunissen, 1993; Kim and Langley, 2000a]. 


With respect to the computation of the objective 
function, there are two basic approaches: a “single- 
epoch” (or instantaneous) approach and a “mutti- 
epoch” approach. A single-epoch approach uses 
one epoch’s observations while a multi-epoch 
approach uses many sequential observations in 
computing the quadratic form of the residuals. As 
far as the performance of ambiguity resolution is 
concerned, the single-epoch approach may not 
attain a higher success rate than the multi-epoch 
approach because the single-epoch approach often 
fails to find correct ambiguity parameters. This is 
due to the difference between a local and a global 
minimum for the quadratic form of the residuals. 
When the observations are _ significantly 
contaminated by biases such as multipath, residual 
atmospheric effects, satellite orbit error and so on, a 
local minimum which is determined using one 
epoch’s observations is apt to be biased. On the 
other hand, a global minimum which is determined 
using many sequential observations over a 


relatively long time span is unbiased as long as the 
behaviour of biases can be assumed to be a random 
process over the time span. Therefore, the single- 
epoch approach should be used with a procedure 
which improves ambiguity resolution performance. 
It has been reported that the performance of the 
single-epoch approach can be improved when a 
linear filter for the residuals is employed [Borge 
and Forssell, 1994] or when the time average of the 
objective function is used [Martin-Neira et al., 
1995]. 


As turns out to be clear (to some degree) from 
the above discussion, the ambiguity resolution 
techniques have different characteristics in terms of 
the computational efficiency of the ambiguity 
search process and hence we can classify the 
techniques according to their approaches for the 
search process. In classifying the techniques, we 
kept in mind two questions: 1) How do the 
techniques describe or limit the ambiguity search 
space? and 2) How do the techniques deal with the 
ambiguity parameters? 


The computational efficiency of the ambiguity 
search process can be improved in several ways. A 
general approach is to reduce the ambiguity search 
space which comprises the ambiguity candidate 
sets. In this case, the correct ambiguity set should 
be retained in the reduced search space once the 
original search space includes it. We can find two 
methods for this approach. One is the search 
domain transformation method which transforms 
the original ambiguity sets into the corresponding 
ones in a transformed space. The reduction effect is 
usually gained through a  “many-to-one” 
relationship between the original and transformed 
sets, and/or through redefining a more efficient 
search space than original one [Abidin, 1993: 
Martin-Neria et al., 1995]. The other method for 
reducing the search space is to define the 
conditional search ranges in multi-level searches 
(e.g., FARA and FASF). This method is based on 
the fact that the ambiguity parameters of lower 
search levels are conditioned on those of upper 
search levels. Some techniques use two methods 
simultaneously (e.g.. LAMBDA and OMEGA). 
Another approach for improving the computational 
efficiency is to find a more efficient computational 
algorithm for the quadratic form of the residuals. 
There have not been many investigations into this 
approach except for studies on the modified 
Cholesky decomposition method [Euler and 
Landau, 1992] because the search space reduction 
methods usually obtain higher computational 
efficiency anyway. Moreover, we have already 
well-known computational algorithms for the 
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quadratic form of the residuals, such as the singular 
value decomposition and the Cholesky 
decomposition. Note that the modified Cholesky 
decomposition is frequently used in the context of 
the ambiguity resolution techniques due to a 
symmetric nonnegative definite matrix in the 
quadratic form of the residuals. Recent studies by 
Kim and Langley [1999b, 2000a] show that 
computational efficiency can be significantly 
improved using alternative algorithms for the 
quadratic form of the residuals. 


As mentioned previously, parameter estimation 
under the integer least-squares problem is carried 
out in three steps — the float solution, the integer 
ambiguity estimation, and the fixed solution. 
Regardless of which ambiguity resolution 
technique is used, the same computational 
algorithms can be employed for the float and fixed 
solutions. On the other hand, the procedures of the 
integer ambiguity estimation depend on the specific 
ambiguity search technique adopted. We can 
classify the ambiguity search techniques into two 
classes according to how they handle the ambiguity 
parameters in the ambiguity search process — an 
“all-ambiguity-search” method (e.g., FARA, 
LAMBDA, FASF and the modified Cholesky 
decomposition method) and an “independent- 
ambiguity-search” method (e.g., LSAST and 
OMEGA). In generating the ambiguity search 
space, the first method uses all the ambiguity 
parameters while the second one uses only 
independent ambiguity parameters which provide 
an unique solution for the system (or observation 
equations) and hence dependent ambiguity 
parameters are determined once the independent 
ones are given [fatch, 1990]. Generally, the 
independent-ambiguity-search method has a great 
efficiency in generating the ambiguity search space 
but computational efficiency of the method is not 
improved as much as expected because of the 
computational burden of the dependent ambiguity 
parameters. Kim and Langley {1999b, 2000a] 
proposes alternative algorithms to overcome the 
problem. 


3. Research Trends and Issues 


Over the past decade, the ambiguity resolution 
techniques have been improved with the ambiguity 
search process made more efficient in order to 
make the ambiguity resolution techniques practical. 
As a matter fact, it is “computational efficiency” 
that the ambiguity resolution techniques have 
aimed at. Unfortunately, the race for more efficient 
ambiguity search algorithms seems to be already 


saturated because current achievements in 
computational efficiency seem sufficient for many 
practical applications. For example, some 
techniques such as FASF, LAMBDA and OMEGA 
provide ambiguity and positioning solutions within 
a few tens of milliseconds according to their 
developers. Such speed is more than satisfactory 
for most applications. Although many of the 
currently available techniques are quite efficient, 
this does not necessary mean that further 
investigations for better ambiguity search 
algorithms are pointless. 


On the contrary, ambiguity _ resolution 
techniques have not usually been described in 
terms of a generalized (or standard) procedure 
which includes: a functional (or deterministic) 
model which describes the relationship between 
observations and unknown parameters; a stochastic 
model which represents the noise characteristics of 
the observations; a quality control scheme which 
handles cycle slips (or outliers); and a parameter- 
estimation scheme which determines ambiguity 
parameters as well as navigation solutions. As a 
result, it is not always clear if the results reported in 
the literature are isolated (“best”) cases or can be 
reproduced at any place on the earth at any time 
during the day. In fact, we need to consider the 
“reliability” of the techniques. IAG SSG 1,154 on 
“Ambiguity Resolution and Validation” had 
investigated the problem during the years 1995- 
1999. Unfortunately, the lack of a generalized 
procedure is still present although progress has 
been made. Furthermore, there is no general 
agreement on the form of the generalized procedure 
up to now. 


As was mentioned previously, the general goal 
of the ambiguity resolution techniques includes the 
“performance” (or success rate) of ambiguity 
resolution. Although no comprehensive study has 
yet been performed, it is evident from many papers 
that concern about this aspect of ambiguity 
resolution has been growing. The issues related to 
performance are two-fold: 1) How do we increase 
performance? and 2) How do we evaluate 
performance? The first issue is involved with the 
qualitative realization of performance, whereas the 
second issue is related to the quantitative 
evaluation of performance. 


To obtain optimal solutions in the least-squares 
estimation, both a functional and a stochastic model 
should be specified appropriately. As long as the 
models are correct, in principle, the optimal 
solutions are not biased and hence we can obtain 
correct solutions. The same is true for the 
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ambiguity parameter estimation because it is based 
on the integer least-squares estimation theory. 
Therefore, the realization of high performance (i.e., 
the determination of correct ambiguity parameter 
values) depends on how correctly we can establish 
the functional and the stochastic models. With 
respect to the function model, specific error 
modeling and parametric modeling for the error 
sources — ionospheric delay, tropospheric delay, 
satellite orbit error, multipath and so on — are 
essential concerns. Although intensive research has 
been conducted in this area, interest is still growing, 
particularly for the applications requiring real-time, 
long-baseline or kinematic solutions. On the other 
hand, stochastic modeling has received less 
attention than functional modeling. However, the 
stochastic modeling tums out to be a crucial 
research topic as the interest in the performance of 
ambiguity resolution increases. The race for more 
realistic stochastic models is in full swing at this 
moment with many recent developments: the 
elevation-angle dependent function approach [Jin, 
1996]; the signal-to-noise ratio (SNR) or 
alternatively the carrier-to-noise-power-density 
ratio (C/No) approach [Hartinger and Brunner, 
1998; Barnes et al., 1998; Collins and Langley, 
1999]; the (adaptive) least-squares estimation 
approach [Wang et al., 1998a; Wang, 1999; 
Tiberius and Kenselaar, 2000] and the real-time 
estimation approach [Han, 1997; Kim and Langley, 
2000c]. Fundamental discussions about the 
observation noise were given by Langley [1997] 
and Tiberius et al. [1999]. IAG SSG 1.179 on 
“Wide Area Modeling for Precise Satellite 
Positioning” will be closely related to the issues of 
functional and stochastic modeling during the years 
2000-2003. 


There are two approaches for the evaluation of 
performance: 1) the performance evaluation 
function approach and 2) the discrimination 
function approach. The first approach tries to 
evaluate the performance of the integer ambiguity 
parameters using the probabilistic properties of the 
integer ambiguity estimators. Teunissen [1998, 
1999] has proposed a performance evaluation 
function, particularly for the integer bootstrapping 
(i.e., sequential conditional integer rounding) 
technique which was adopted in the LAMBDA 
method. On the other hand, the second approach 
tries to measure discrimination power between the 
best ambiguity candidate and the second-best one. 
Although some discrimination test procedures have 
been proposed during the past decade, more 
teliable test procedures which can quantify 
discrimination power have been proposed only 
recently [Han, 1997; Wang et al. 1998bj. 


Compared with the activities for the qualitative 
realization mentioned previously, the research on 
the quantitative evaluation is less intense so far. 
Nevertheless this topic is ripe enough to start a race. 


With an increased level of research on (real 
time) long-baseline (kinematic) applications 
recently, some of this research has shown up in the 
ambiguity resolution arena. Although we have 
already noted the extensive ambiguity resolution 
research over the past decade, this research was 
usually directed towards post-processing and static 
applications. But all the issues mentioned above 
(i.e, computational efficiency, reliability and 
performance) are inherent in the real-time, long- 
baseline, or kinematic applications. For example, 
functional modeling tums out to be difficult 
because of the decorrelation of such biases as 
ionospheric delay, tropospheric delay and satellite 
orbit error in long-baseline situations, and the 
quasi-random behavior of multipath in kinematic 
situations. As a result, stochastic modeling also 
turns out to be more difficult. On the other hand, 
those problems can be handled relatively easily in 
static and short-baseline applications. 
Conventionally, the multi-reference _ station 
approaches have been proposed to overcome the 
problems in (post-processing) long-baseline (static) 
applications (see Fotopoulos [2000]). The 
capability of real-time and kinematic processing 
has been investigated using the virtual reference 
station approach [van der Marel, 1998; Odijk et al., 
2000] and the parametric estimation approach [Kim 
and Langley, 2000b]. 


4. Concluding Remarks 


A brief review of the research trends and issues 
on ambiguity resolution and validation has been 
presented in this paper. At first, we took a look at 
the work carried out by many research groups from 
all over the world. Since it is not always clear if the 
results reported in the literature are isolated (or 
ideal) cases or can be reproduced at any place at 
any time, comparison of the different ambiguity 
resolution techniques based on the literature is not 
always reliable. To avoid this pitfall, we tried to 
give our viewpoints in classifying the techniques 
solely in terms of their characteristics. This may be 
helpful for a general understanding of this research 
area. 


Considering the general goals of the ambiguity 
resolution techniques, it is the computational 
efficiency of the ambiguity search process rather 
than the performance that characterizes the 
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individual ambiguity resolution techniques. For at 
least a decade, there has been an_ intensive 
competition to develop more efficient ambiguity 
search algorithms because computational efficiency 
was essential to make the ambiguity resolution 
techniques practical. As a result, current 
achievements in computational efficiency seem 
sufficient for many practical applications. This 
trend has switched the race to the performance of 
ambiguity resolution. Currently, work on two 
issues related to the performance is in full swing: 
i.e., the qualitative realization of high performance 
and the quantitative evaluation of performance. 
Both issues are equally important in terms of 
accuracy and availability in designing a system. 


We are facing a new challenge in the demands 
of real-time long-baseline kinematic applications. 
There is an attractive driving force in the challenge 
because all the issues related to the ambiguity 
resolution (i.e., computational efficiency, reliability 
and performance) are inherent to these applications. 
These applications require the integration of the 
knowledge and experiences that we have obtained 
from the conventional approaches for short- 
baseline and long-baseline analyses, for static and 
kinematic applications, and for post-processed and 
real-time evaluations. 
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Abstract: To determine attitude in real time, it is very important to reduce computational efforts by 


making searching space efficiently. There are many algorithm developed for real time RTK or 
Attitude determination. In those one, We used LSAST(Least Square Ambiguity Search Technique) 
which is introduced by Hatch, and determine attitudes using baseline constraint to reduce searching 


space effectively using constraint equation of ARCE(Ambiguity Resolution with Constraint Equation). 


Using Ll carrier phase, we could determin attitude with 0.4 deg. accuracy in elevation and 0.2 deg. 


in azimuth 


Keywords: LSAST, ARCE, baseline constraint, integer ambiguity, attitude determination 
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Abstract: With GPS LI carrier whose wavelength is about 19 cm , very high accuracy can be achieved. If integer ambiguities of carrier 


phase measurements are removed, it is possible to realize precise positioning within a few millimeters. In order to determine ambiguity 


set, however, it needs large computation. For the reduction of computation, there is way to reduce search range of integer ambiguity. 


This paper presents a method that is able to improve a Ambiguity Resolution with Constraints Equation(ARCE) algorithm. By selecting 


the allocation of satellites observed and the operator of satellites differencing, it is proved to be able to reduce a search range. 


Keywords: GPS, Carrier Phase, Integer Ambiguity, Search Range 
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Abstract: Nowadays, the range of GPS navigation is being extended indoors. Pseudolite is thought to 
play a key role in this area. The indoor navigation system is one of pseudolite application.After rigorous 
efforts, we have successfully implemented the indoor navigation system using CDGPS. Error analysis shows 
that the RMS value of horizontal kinematic positioning error is below lcm. In this system, a miniature 
truck is equipped with all the devices for real-time operation. Our final research goal is automatic control 
of vehicle using this system. To achieve this, we have mounted 2-GPS antenna on the truck to obtain 
the yaw angle. From experimental results, we can concluded the indoor navigation system can be sufficient 
position and attitude sensor for miniature vehicle control in factory, hanger, indoor amusement park, and 


etc. 
Keywords: pseudolite, 
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Abstract: This paper presents a practical real-time Carrier Differential GPS(CDGPS) using low-cost C/A 


code GPS receivers and a communication link based on mobile phone and internet. Such a communication 


link is ideal for a high data rate. Since the mobile phone and internet systems are universally available 


nowadays, use of such a system extends the relative positioning application to virtually unlimited areas. A 


Kalman filtering technique is used to estimate integer ambiguity states in real time. Field test results are 


presented to evaluate the performance of the proposed system. 


Keywords: RTK, OTF, CDGPS, Internet, Kalman 
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GPS Control Segment 


“Fakon AFRO wee 
Colorado Springs’: Pies Xwajatein 
Master Contra 6 soepeah ents 
Monitor Station. ¥ fo ; aes eee 


ee oo arma 
po Hawaii oy 


u ik scension Istand 
Monitor Station ©... Mdéniter Station « 


ae 
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GPS Control Segment (Al 4) 


e Replaces 50 Hz Navigation Message 
Data Stored in Memory 

e Broadcasts the SIS PRN Codes, 
L-Band Camriers, and Uploaded 
50 Hz Navigation Message 


oh SATELLITE 


Ui &ts 
DOWNLINK 


ie Gasren CONTROL STAION 


MONITOR STATION 


GROUND ANTENNA 


e Checks for Anomalies 
e Sends Raw Observations » Computes SIS Portionof UERE ° Sends New Upload 
to MCS © Generates New Predictions to SN 
¢ Builds New Upload 
© Sends Upload to GA 


Control Application Lab. 


GPS User Segment (Receiver) 


Time Keeping) += 4! | 


* Low DynamicS 


* Medium Dynamic® 


* High DynamicS 


Ajou University 


Ld 
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“@. Ajou University 


cl & 


1,4NAd PS 
(Aaa, o2SNA) 
1. Down Conversion 2, GP3 AAA We 
= 3, GPS AA NAQ 
2, OS na 
4. 40A AN AX 
: (28 fads) 
Signal 
pipcessing Data Processing 


2.20 44 Iola 
1. SY GIOLE 


Control Application Lab. a s Ajou University 


2) A} cl(Pseudo range)o] SA Z! cl 


© WAHel 


PR=c:-tT=c:-(T,,.t+OT)= R+B 


true 
c: AAS 
: GPSA AD FA] AAS Pest ATA 
: & ATAl2t 
: AAD FADO ALAA] AO! 
: & Ael 
: AAAS) SA Hel 
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=—200> 


AA cl 2] SA al el — Code Tracking 


ASA SUE Code} SAE Code “I = all 


LH OW Ad a8 2 St Ch. 
Al] OIA Ale Code AIS : p(t) 
AAIIONA, HAGE Code A'S : P(t-7) 


- SASAIA2AA SAI AGE rSRS 
— Delay Lock Loop(DLL)2|'% : PN Code Correlator 


Control Application Lab. « Ajou University 


PN Code2| & 2(Correlation) 


F(t) = vale) 
ptt) vat) LOWPASS va (ty 
Neen . FILTER . bg Say’ (fame as 
eae -T T =" 


RECEIVED 
CODE 


REFERENCE CODE 
pit -T) 


Vath) = pct) p(t -t) 


Val(t) = Ave [ pipp(t- t 
Fit. 


A jou University 
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GPS Vel 


1OMAASBYLSSYLSSAAYBIA, PPS) 
— 3.3M(11 NSEC)2} 2] AFA 2] 2 XKGDOP=3 J/Z) 
— 10.23 MBPS code Al A! (P-Code) iH 3 


30M 911 SL S(HEEA AYIA, SPS) 
— 1.023 MBPS Code AlA& (C/A Code) HS 


Ba AWIAS| ASA et 
- JSSOPY)AS 


- HESAMBIAY BVS RAIS 


— Selective Availability(SA) HS :2000H 52 122 aAll 
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GPS {!S 4A GIS) 


© OAC! ASBOHL} DEO Chet Bola 
- 2A AE & (Spread Spectrum) 4! S 01S 
— Pseudo Random Noise (PN) P-Code, C/A Code. 


Aa @Aaes ASA Hae SAO Saeet yt 
— L-Band Radio Frequency(RF) 4!S 4S 
- ASA OSS AUAA Sts Fld 2712] L-Band A 
— LI AlS: 1575.42 MHz (=154 X 10.23 MHz) 
L2 AlS: 1227.6 MHz (=120 X 10.23 MHz) 


lA ASS MsS 


lSSAOl Slot ts 
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ee AT fs 


GPS ASI FPA:Lids 


S,(t) = ApXP,(t )D,(t )cos(27f,t + $)+ AcXG,(t)D,(t) sin 27f,t + $) 


Sf, + Li Carrier Frequency (1575.52 MHz) 
Dt) + l& (binary) 84 G0! E1( 50 BPS) 


XP,(t) : P-Code (10.23 MBPS, 49/72 
XG,(t) : C/A-Code (1.023 MBPS, 2£| 1 MSEC) 
Ap ,Ac : P-CodeSt C/A Code SUS SS ASS 
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2 
SIN X é 


x ' 
SPECTRUM 3 


PRLIIP 
ee dy 
LLXY Ly Lez ‘ 
| 1227.6 MHz ! 1575.42 MHz 
| | 
I t 


[pene jamz | 


L2 ul 
SIGNAL SIGNAL 
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GPS AS2 SS 


(SPACE VEHICLE) 
" EPHEMERIS (Xs, Vs, Zs, ts} 
CLOCK BEHAVIOR (dts) 
HEALTH & STATUS (ALL SATELLITES} 


SYNCHRONIZATION: 
CODE 


ALMANACS 
(ALL SATELLITES) 
MESSAGES 50 BPS (1500 BITS/FRAME, 5 SUBFRAMES (300 BITS)} 


(Xs, Ys, 2s, ts) 
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at Gl 0] A (Navigation Message) 4! S 2] 2 


50bps 0/ 2! GIO! Ef 
& 1500 BITS, 322!) 2O| 


Li = (Navigation Message) 
GPS AIA2| SAIS! Al Zt Wheto El 
HANDOVER WORD(HOW) : C/A-CODE& ¥ Ej P-Code 3 S01] 01S 
ALMANAC: @ GPS #i4 2) GAME! SAI BAYS AEH Weta 
GPS A 2OS4E UTC Al2to] He tetolal 
CS! Frequency ASAMS Sle OSS ATAA SA Met olel 
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at 4 G|0|/G % (Navigation Message Format) 


SUBFRAME }—_—— 10, 30-BIT WORDS ——_- 
# 6S 


[ TLM BLOCK | | 


CLOCK CORRECTION 


BLOCK II 


EPHEMERIS 


BLOCK Il 


1 FRAME 


30s, 1500 BITS —___-__--__4 


EPHEMERIS 
MESSAGE + IONOSPHERIC MODEL 
HOW 


+ BLOCK Il (25-32) [PAGE 18 OF 25] 


GENERAL 


HOW | BLOCK Ill 


ALMANAC & SV HEALTH [PAGE 25 OF 25) 
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22els 


MSA SH GIGS OSH AAS! Al2tD AAl WlAt 
— Clock Correction Parameter: SA! Sl AALS] Al 2k BA at Al At 


— Ephemeris parameter: SA! Fl 42! Sl Al AA 


4H OlAFS] GPS FAO Het A AtAecl SAAS SE 
AZ(XY,Z)2 AlAlL ABS 7S 

— Deterministic Solution Method 

— Kalman Filter Solution Method 
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= AA! (Pseudo range Measurement) 


PR, 


(X,,¥,,Z;) 


(XYZ): 


Control Application Lab. 


¢ Instantaneous Solution 
— Iterative Method 


* initial state required 
— Direct Method 


* no initial state required 
* Koch(1980), Bancroft(1985), Krause(1987). Leva(1996). Biton(1998) 


¢ Sequential Filtering Solution 
— Kalman Filter 


Control Application Lab. - eB A & Ajou University 
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Instantaneous Solution (Koch Algorithm) 


INAS 4HS OH AAAe] SBAEZ 
(x; -x)? +(Y; -y)y +(Z; -z)= (Pp; -1)? for i=1,2,3,4 (1) 
(4-unknowns, 4-equations) 
OMA (Hs Aso, 


x? —2x;x+x7 +y? —2yiy ty? +27 - 22,242" =p? —2p,t4+7? (2) 


AAS Hop Asso GCSiw Ac. 
ax+byy+ez-djt=f; for j= 1,2,3 


AQ3)s tO] Aso Sa, 
x=k,t+l, 
y=k,t+l, (4) 
z=k,t+1, 
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Koch Algorithm 


(4)3 4 Oo CQ otOl H10 CHst AAO & LEU 
at? +2Bt+y=0 


21S MHS. ( [r[<150km_ ) 
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‘Kalman Filter 


¢ Dynamic Model 
x(k +1) = x(k) + w(k) 


* Measurement Model 


z(k) = Hx(k)+v(k) 
Where, 


[x y z SOs B A in ECEF Coordinates 
: State Transition Matrix 


: Observation Matrix 
: Process Noise 

: Measurement Noise 
: Time Sequence 
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Kalman Filter Process 


Initial states %. 


medion CES es Be - Measurement Update. 
_ Time Update a 


Ree ronnie rat 8 “(correction)” 
ae (Prediction) = ——— 


ee ee | RAHA AAR 


| a= @.RO +9, Pe en 


PaU-K HI 
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tO} 
—_ 


1. GDOP : ASA St GPSFl 42'S] Dl ot el 
(Geomatric Dilution of Precision) 


UH Al af Ed 


2. UERE : JA He] SALI 3B) 
(User Equivalent Range Error) 


et 8 © At (1-sigma) = GDOP x UERE 


Oo 


K Ajou University 


Control Application Lab. 


C/A Code 2| At Hel SA & At (1-Sigma, meters) 


SA Sa SP 


eis SAS He ae 
ource 
| Bias {Random | Total Random 


Ehpemeris Data 
A Al Al 
Jes AMAA 
ea anaes 
Multipath 
ANI BS 
UERE 
Filter = UERE 
+2 2 Xt(VDOP=2.5) 


+= © ii(HDOP=2.0) Al(HDOP=2.0 


“@ — Ajou University 
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GPS HA WAS] ASSO Het S 


RANGE ERROR) _ig¢ +lo 


GDOP 
lo. RANGE ERRO 
la + POSITION ERROR 
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GDOP2| JS 


GDOP2Y] 32 
. _ _4—-Det otto 
GDOP = rite Bes a(10) 


¢ GDOP Hck 2] =a 
cov(x) =E(4x-4x")=0:[G'G]"! 
i gt Oxy Oxz Ont 
oy Dyz Oy 


ae 
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GDOPS| SF 


_ ABA R OAL) 
HDOP = -orip5jzj BO at(1o) 


_ ARIA MEIRIOA(1O) 
VDOP = -ofA3]2]/ SS At(10) 


_ 3 —D9 a SRI 
PDOP = ~sixpe]Ba ert) 


Al Z+2 AC o) 


TDOP = “oppyaaw eat o) 


_ _ 4— D8 24t(10) 
GDOP AAAS 2 A(1o) 


Control Application Lab. 28 Ajou Universi 
BA ty 


412] GPSAA B2tOll GE PDOPS| 4 st 


ZENITH 


MINIMUM PDOP 1,633 


MINIMUM WITH 10° 
MASK ANGLE & 2.0 


© ONE SATELLITE OVERHEAD 
@ THREE SATELLITES 120° 
APART IN AZIMUTH, 
)° ABOVE HORIZON ¢ 48 20 


ANGLE ABOVE HORIZON ¢ 
(DECLINATION) 
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2 # © Xt (1-sigma) = GDOP x UERE 


So 


UERE = 5.1m , VDOP = 2.5 , HDOP = 2.0 


(16) = VDOP x UERE = 12.8m 
TH © ZH 1c) = HDOP x UERE = 10.2m 


&. Ajou University 
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GPS2| & CH St+(Modernization) 


es Ajou University 
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GPS2| & CH st(Modernization) 


History 
- 1996S OSS AAS YE (PDD on VS. GPS Policy) 
OA Ol ASA! GPS SSA, 1044 OL V2] SztSA MH 


- 19974 Interagency GPS Executive Board (IGEB) 24: QO] DoD 2} DoT2] SSIF 
SIZIS AHIA(SPS)S AS MES AS BIB US SS AIS 


- 199714 Ol] National Defense Authorization Act for FY1998 
inZs 4!S (M-code) 


- ae 19994 Ol 42 Gore & B(GPS Modemization Initiative) 
2 Tt44(1227.6 MHz)Oll 2! 2! C/A-code A! 77} 
cd DEAL Mi 32) QIZS AlS FI Safety-of-Life S) 


- 19994 IGEB Working Group 
A392] O12 GPS 41S & 1176.45 MHz (L5)ol St 


- 200014 0] HES, SAGMH SH - 2000 58 1 SAHA AAI 
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instantaneous Error (meters) 


Selective Dvalability(S):0 oH Al (“Magic Moment”) 
20004 52 12(0:00 EDT, 4:00 UTC) 


i 
SPS CEP AFTER TRANSITION: 78 meters 
7 SPS SEP AFTER TRANSITION: 4.6 meters 


-— Horzontal Error (rreters} 


ae oe ia Vertic at Error imetars} 


L ee Vite bea woettaie oe 


Time of Day (Hours UTO) 
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Selective Availability(SA) of XH] Sa} = 


SA 4S Al (2000. 2.15) SA oll Ail = (2000. 5. 9) 


Nedirection (meter) 
Nedirection (meter) 


0 
E-direction (meter) 
1a SA AWA SA 
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GPS & C1 St(Modernization) 01 C= 
Kt ACH GPS Al S 2} BHAI 


L5 L2 


, P(Y) P(Y) 
emorsNs ame aw 
C/A 


L20l 
C/A lS FD} 
(20034 ZALIHAI) 


1 
C/A 
C/A 


Al3 Gl2tAls(L5) 
Ainz2sas(m) 
(2005 EE APTA 
1176.45 MHz 1227.6 MHz 1575.42 MHz 


* C/A: BAS, Ply): HHL2BAS, Midnese ds 
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GPS Block IIR Satellites with New Signals 


f tet yn weet 
* Possible increases 
PY) and C’A-coda 
povrert OVE? Current 
specifications 
{-163 and -160 dBW 


* New Earth coverage 
military code 
(Mi, code) 


Up ta teetve HA satelhtes 


if 


© Possible increase in 
PY)- code power 
aver current 
spacificarans 
{166 dBW; 


© Addition of C/A-code 


© New Farth coverage 
nultary code 
iM, code) 
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GPS Block IIF Satellites with New Signals 


oh se gah gis easy 
* Now M, code added 
* Complete Picture 
Undetnes 
Boemy des 


Control Application Lab. 


Twelve IF satelites 


L5 


Lt 
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By Seas 
* New robust Civitan 
Navigation Signatin 
ARNS band 
¢ Centered at 
VI7G.45 Miiz 


eo code 
; Code added 
@ Complete Picture 
fined wire 


Ajou University 


GPS #0] GE l2tS AUlA(SPs)2I 
3S Ba (derived from 1995 SPS Signal Spec) 


Typical Range Error Magnitude (meters, 1-sigma) 


Error Source with sa | Without | Multiple C/A-code signals 
SA (L1, L2 and/or L5) 


Selective Availability 


Atmospheric Delay 
Tonospheric 


Tropospheric 


Clock and Ephemeris 
Receiver Noise 
Multipath 


Total User Equivalent Range 
Error(UERE) 


Typical Horizontal 
DOP(HDOP) 


Total Stand-Alone Horizontal 
Accuracy, 95%(2-sigma) 
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GPS & C1 &t(Modernization) 2 3 


Activity Implementation Date 


SA set to zero May 2000 
GPS IIR enhancements 2003 — 2006 
-M-code on L1 & L2 

-C/A Code on L2 

GPS IIF enhancements 2005 — 2010 
-C/A Code on L2 

-M-code on LI & L2 

-L5 


GPS II] enhancements 2010 -— TBD 
-C/A Code on L2 

-M-code on L1 & L2 with greater power 

-L5 

-Future Capability 
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a 


GPS 3 &(Application) 


- AA Ss (Static Applications) 
_ STI, ISH, GIS 


= anematie Applications) 


ast 7) 
Ket, Aut SB Oe 


|2} & S (Time Transfer) 
2t S 9|(Time Synchronization) 


Control Application Lab. 


t KSASAABGIS) PS 
Al: LoOlS, AAS, ASN 


A|2t S| 
- NSA S| 
- SSA ASI, 2! 


« Ajou University 


Control Application Lab. 


=272— 


oe Ajou University 


40 


ru 
oT on Sot 4H my 


S| KHAIBZA P| 24 Et AE 


2o, e468 


2 4&5, 2OSAl (Battle Management), # At34, AAAI AS 


2 


40 
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AAICH GPS SS =O0t 


¢ GSP + Internet + OIS SA! 
— GPS AL SSA SA (SnapTrack) 
— MASSA Ut 


¢ GPS +E119 
— 119 PA AS Al AS Sal QA! 
— O}]= E911 Mandate 


125m(2dRMS)OILIS! IRIS 4A 679 ABBE W 
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GPS Internet Web Sites 


F pas 
ee by, Reter H- DanayiThe University of http://wwwhost.cc.utexas.edu/ftp/pub/grg/gcraft/notes/gps/gps.himl 


Richard Langley's Navstar GPS Internet Connections | http://gauss.gge.unb.ca/gps.internet.service.html 
ttp:// -inmet.com/~pwt/gps-gen.ht 
Paul Tarr's GPS WWW Resource list a ae a he 


http://galaxy.einet.net/editors/john-beadles/intorgps.htm 


John Beadles's Introduction to GPS Applications 


Wolfgang R. Dick's GPS-Information Centor http://www. potsdam.ifag.de/server/geodesy/tech-gps.html 


Ww http://www. gpsworld.com/ 
GPS World on-line gpsworld.com. 
N h Si http://www.navtech, 3 if 
freveensenines | jeminars p avtech.gps.com, 


http://www.navcen.uscg.mil/ 


USGS-Navigation Center(NAVCEN) 


tp://tycho. ; mil html 
US Naval Observatory hitp://tycho.usno.navy.mil/gps. htm! 


The Institude of Navigation http://www.ion.org/ 
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GPS- SADA NS FSBe 
Set MHA MS 
Coverage, Accuracy, Availability, Affordability 
HAAO|D ISAC So 
“GPS BUish"Ss set a 
GPS S22] OSs & 


WSsAN, Ceo Cl 
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GPS FAl7| ZjStt set 
(GPS Receiver Technology & Trend) 


Cc =. 
Sr 


dima _ chung@samsung.com 


St 


my a a 


GPS +4al7|2| +a 2 HW 


- GPS FA17/ 9] FA 
— GPS FAl7/2I 7S 
— YEE ol GPS FAI7/9] FAI 


° GPS = 4l7| Hardware 
— RF/IF 
— Digital Correlator 
— Tracking Loops 


a 


GPS =4II/2| BIS 


* Pseudorangeo| Sa 

> SM FAO] GAst GPS AA ASS 
Transmitted Time ? 

< GPS PRN 4! 2&9! replicas SAGA Ol 
replicaJ| 4! | GPS PRN SS ASASCEe 
Track?! = 2S = et (Tracking Loop). 

< HL St #210 Ol replicaSl| Phase 2 PRN 

1ms epoch counts SAO SA otO 


Transmitted Time3 7 &. 


=216= 


Transmitted Time2| BA 


GPS Channel #1 
GPS PRN XAIS 


X\ PRN X Replica 
Tracking 
ere 


wed COE 
maa NCO 


GPS Channel #N 
GPS PRN WA! 
PRN W Replica 


Tracking 
cop Control 


2t 2{ Aol] CH St 
Transmitted Time 


am 


ELECTRONICS 


Sides Navigation 

; \—7| Processing Processin 

RF/IF i i | - Acquisition 2 
| | - Tracking 
- Synch- 
| ronization 


& 
SAMPLING 


NGHANNEL | 
ORRELATORS | 
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RF Input 
from 
GPS 


Antenna 


RF 
Filter 


Power 
Regulation 


Al 9| > 28(Mitel) 


Dw9255 
Band Shaping 
Filter 


Logic Analyzer 
Access Ports 


GP2021 
12-Channel 
GP2010 Correlator 


RF Front End & 
Support ARM60-B 
Functions 32-bit RISC 


PC Link serial ports 


RF/IF 


¢ Translate spectrum to convenient IF frequency 


¢ Establish thermal noise floor for system 


¢ Reject out of band interference prior to correlation 


« Reject internally generated spurious responses 


=-218- 


| RF/IF Down Conversion. 


Antenna 


Digital 


Reference 
Oscillator 


Digital 
Digital 
converter 


integrate t 
&) & dump 
©) 


& dump 
integrate 
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Code Tracking Loop(Example) 


Integrate & 
dump 


Bk shift register 


Signal 


Processing 


controled 
oscitator 


Code Loop Discriminators 


DX Ces ~ lis os + (es — 1s )Oes 


— Dot product power. 


» (2s + Q2,)- dD U2, - Ors) 


— Early minus late power 


Vo.2 2 als 2 
fee AeA Pom nen's 


— Early minus late envelope 


v2 


Often 


— Early minus fate envelope normalized 
by the early plus late envelope 


~ 


2 
1) 


705 


A 
: 
; 
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Code Correlation 


Incoming signal 

Replica signal 
Early 
Prompt 


Late 


% Chip Early 


Normalized correlator output power 


Code Correlation(Cont.) 


Incoming signal 


Replica signal 


Early 


Prompt 


Late 


Aligned 


Normalized correlator output power 


1-14 1143/4 
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Code Tracking Loop Design 


¢ Design Consideration 


— Discriminator design: Computation Burden vs 
Performance 


— Code Loop Filter Design 
— Predetection Integrator Design 


— Correlator Design : Combined Early-Late 
Correlator, Tau Dither Technique, Variable 
Correlator Design, Narrow Correlator 


— Summation and Update Rate : Doppler Aiding 
from Carrier Tracking Loop 


am” 


ELECTRONICS: 


Carrier Tracking Loops(Example) 


Digital Signal 


Processing 


Carrier 
Loop Filter | 


Prompt Replica Code | 


Numerical 
Controlled 
Oscillator 
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Carrier Tracking Loop 


¢ PLL(Phase Lock Loop) 
— Minimize Carrier Phase Difference 
— Accurate Tracking 
— Sensitive to dynamics 

¢ FLL(Frequency Lock Loop) 


— Minimize the Difference of Carrier Phase 
Difference in Time domain 


— Less Accurate Tracking 
~— Less Sensitive to dynamics 


Sign (Ips) “Ors 
Tos" Ops 
Ops lps 


ATAN(Qps Ips) 


—> sind 
> sin2p 
—> tand 


76 


Costas PLL discriminator output (deg) 
BRSee,.8288 


Abd Tt PAROS 


_ 
ny 
oa 


+150 


“1805 


' 


i i? 
i iF 
if 
ie 
f H 
“ | 
| “a } 
: i 
a : 
f i 
fa i 
H 
i 
i 


eee ae OO io” : § 


True irteut error (degrees) 


18 
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Frequency Lock Loop Discriminators 


Discriminator Algorithm Output Frequency Error 


sign(dot)cross sin{2(¢; - ¢,)] 
foal LP 
where 


dot = Jpsi * ps2 + Orsi * Opsz 
cross = Ip * Ips - Opsi * Qps2 


cross sin[(¢, - ¢,)] 
(t,- 4) Pare: 


¢,- ¢, 


ATAN 2(cross,dot) G, -1,)360 
oe = NEA by 


(1, - 1,)360 


Carrier Tracking Loop Design 


* To tolerate dynamic stress: 
Discriminator - FLL 
Predetection time - Short 
Carrier Tracking BW :Wide 
¢ Carrier doppler phase measurement to be accurate: 
Discriminator - PLL 
Predetection time - Long 
Carrier Tracking BW :Narrow 
* Design Consideration 
— Discriminator Design : FLL/PLL 
—» Carrier Loop Filter : Band width, velocity/acceleration 
sensitivity 
— Variable Tracking Plan 
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GPS =Al)| Software 


¢ Signal Processing 

— Acquisition 

— Tracking 

— Synchronization & Navigation Message Collection 
¢ Navigation Processing 

— Satellite Acquisition and Tracking Planning 

— Navigation Message Decoding and Processing 

— Satellite Position Calculation 


— Navigation Computation 
¢ Pseudorange Formulation | 
* Weighted Least Square Method, Kalman Filtering 


Penusuned Abel UKEP APL 24 


ELECTRONICS: 


GPS #4l7| SW 7A 


° 1kHz Tasks 
Correlator read, accumulate 
Correlator Control 


¢ 50Hz Tasks 


Tracking Loop Update 
GPS Navigation Message Collection 


°1Hz~ Tasks 
Navigation Processing 
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Signal Acquisition 


* 1/2 chipEt | 2] code binS coded +4 
* T=dwell time per cell (1.0ms~20.0ms) 
¢ Dwell timeS cto] Wo S AlAtS Threshold?! at2t 4| a2 


> 


| Cee) Weanparr 


a | 7 


<uig Aguanbel4 


Bit Synchronization 


* Message Bit 2} S7]| 


1 ms C/A code 


C/A Code 


17 18 19 20 


20ms { message 1 bit} 
Message Bit 


t 


Bit Lock 


=2o0 


Frame Synchronization 


¢ Message Subframe Synchronization 
Subframe2| Word 12} Word 201] CH o}0j C}S 2] ParameterS0| DS S| 
Frame LockS “121 
1. TLM preamble (10001011) 
2. HOW subframe ID 
3. HOW zero bits (bits 29 and 30) 
4, Parity check 


i Message Data Bit 
57 58 59 60 


| Word 1 Word 2 


Frame Lock 


an” ASAT SMAPS 


ELECTRONICS 


Message Subframe 


Navigation Processing 


¢ Navigation Message Decoding and 
Processing 

* Satellite Position Calculation 

¢ Navigation Computation 
— Pseudorange Formulation 
— Weighted Least Square Method, Kalman Filtering 


¢ Satellite Acquisition and Tracking Planning 
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Navigation Message Decoding and 
Processing 


¢ Navigation Database 


SV Clock & 
Ephemeris Database 


Decoding 
Routine Almanac Database 


Ephemeris 


Collection Buffer ; 
areas meas : UTC Parameter 
oe Data | Database 


Collection Buffer | 
ase ea lonosphere Database 


Satellite Position Calculation 


¢« GPS Satellite 
Position 
Calculation 
using Almanac 
and Ephemeris 


\ Satellite Position 


Reference Position 


Geo 


S2885 


Navigation Computation 


¢ Pseudorange Formulation 
Raw Pseudorange 
< SV Clock Correction Calculation 
<- lonospheric Delay Correction Calculation 


< Tropospheric Delay Correction 
Calculation 


Navigation Computation 


® Pseudorange to /s satellite 


P,=|,-7[+0-b, +4, 


@ Pseudorange Estimate to /’s satellite 
B=, -A|+e-b, +2, 


. p;=Pseudorange to /’s satellite 
@ Pseudorange residual r;= ECEF Position Vector of /’s satellite 
- r= ECEF Position Vector of GPS receiver 
=p-p=-(-H" 1 Ag c = Speed of Light 
TREE : ee ; 5, = Receiver Clock Bias 
€,= Pseudorange Error 
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Navigation Computation 


* If there are N Pseudorange measurements : 
Ap, = h,Ar, +h, Ar, +h Ar, +c: Ab + Ae, 
Ap, =h,,Ar, +h,,Ar, + h,,Ar, +c: Ab + Ae, 
M 
Ap, = h,,Ar, +h,,Ar, +h,Ar, +c-Ah + Ae, 


* Weighted Least Square Solution 


At =(G'R'G)'G'R'Ap 


«a 


Weighted Least Square Method 


eof ME vy SSR He 2, SSR ASH 


z= 
“Ty 


e/S a[rsste Zoe 1 S Aste, 
(Least Square) 


AS 
x=(H'R'H)'H'R"z : weighted Least Squares 


ELECTRONICS 
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Kalman Filter 


X,., = AX, +G,w,, w,~N(0,Q) — : System Model 
Z,=C,X, +V,, V, ~N(O,R) : Measurement Model 
¢ Time Update/Propagation : 


S,= 48, Ply = AB’A+G,0G] 


¢ Measurement Update : 


St = 5,4 ICTR '(z,-C,3;), PX = Po +CTR"C, 


Kalman Filter Implementation 


¢ System Model 


X44) = AX, +G,w,, w, ~ N(0,Q) 
x, =|, (k) r(k) nk) 8) vk) vy (kv, CE) d,(k)f 


1, > Recv. Position in ECEF Frame 
V,,, > Reev. Velocity in ECEF Frame 
b, — Recv. Clock Bias in meters 


d,, > Recv. Clock Drift in m/s 
QO — Process Noise 
Covariance (depending 
on Receiver dynamics) 


20s 


Kalman Filter Implementation 


¢ Measurement Model 


(6, @-*, 0) +, @=7, OF +6, O=7, OF +80 | fe, 
2, =| V1, OF +6, =, OF +6, 7,0) +20 |) & 
M 
(r,, (k)=7,, BY + (ry, Kr, OD) + (ry, 7, OY +8, (K) 
= f(r, (K).b,()) + Vy 
= f(n)+y, 


M 
E 


m 


| 4 =h.® 1 1) 6&0» vk) 4) 4,(of 
— Extended Kalman Filter using Measurement Model Linearization 


— H(x)= oe 


a 


ELECTRONICS. 


Kalman Filter Implementation 
¢ Time Update/Propagation : 


Xyu = AX, Py, = AP? AT +G,OG; 
¢ Measurement Update : 


Sf = 84 PSH (SR Ue (R)). PE) = PO + HT (RRA) 


Hye TZ 


sede" 


Kalman Filter Implementation 


¢ Implementation Consideration 


— Computational Capacity of Microcomputer 
— Numerical Stability 


—Covariance Matrix should always be 
Symmetric and Positive Definite. 
— UD Algorithm 


GPS F4l7/21 4 


- UHCI GPS EAI7/ 2] Als 
¢ Basic Functionality Test 

¢« Parametric Test 

¢ Environmental Test 

eJFEF ALY SE 


| os Oe ee | 
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QlHFAol GPS FAl7/2I A 


¢ C/NO : Carrier-to-Noise Density 


C/NO(dBHz)= Carrier Power(dBW) — Noise Power 
Density(dBW/Hz) 


— Carrier Power(dBW) : AntennaS 
GPS 4! 2! Power(Antenna Gain 
01 S73) 
— Noise Power Density : 10Log(kT)+Noise Figure 
Oil) -160dBW —(-200 dBW/Hz) = 40 dBHz 
—+ GPS #Al)| AlS Al C/NOES Monitorot 3 Al GPS 


RFS AS 8 OIA OFS BCH, 


a” 


ELECTRONICS, 


TTFF(Time To First Fix) 


* TTFF: GPS =AID| It SH ob Al ot 
O{ Validét Navigation Solutions 4aeZe 
CHORAL ASTE Al2t 
— Hot Start TTFF : Approximate Position, 

Time and Valid Ephemeris Present 
— Warm Start TTFF : Approximate Position, 
Time and Valid Almanac Present 


— Cold Start TTFF : No Position or Time 
Present 
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TTFF AS 


e Q8t4 Ol Hot Start Flow 


: GPS = 4!9| Power Up > Acquisition > 
Tracking/Message Decoding —> Frame 
Synchronization > Navigation Solution 
Computation 


¢ Hot Start TTFF 
: Power Up/ Acquisition/ Tracking Time 
+ Frame Synchronization Time 
+ Navigation Solution Converge Time 
H)y42+42412=9% 


Pensuncd AAMT PMA At 


ELECTRONICS 


TIFF AS 


e« 214+ 01 Warm Start Flow 


> GPS +419] Power Up > Acquisition > 
Tracking/Message Decoding — Frame 
Synchronization/Ephemeris Collection > 
Navigation Solution Computation 


¢ Warm Start TTFF 
: Power Up/ Acquisition/ Tracking Time 
+ Ephemeris Collection Time 
+ Navigation Solution Converge Time 
—6hAWDE B+ BOR +1 ZH37R 


am 
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Navigation Accuracy 


SPS(Standard Positioning Service) 
— Predictable Accuracy 

< 100 meters horizontal 95% of time 

< 156 meters vertical 95% of time 

< 300 meters horizontal 99.99% of time 
< 500 meters vertical 99.99% of time 


Basic Functionality Test 


Check the Basic Operation of GPS 
Receiver 

— BIT(Built-In Test) 

— Measure Test Points 


— Antenna Connector Matching and Return 
Loss 


— Predetermined Factory Test Scenario 


a 


ELECTRONICS. 
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Parametric Test 


¢ Evaluate the performance of GPS 
receiver for TTFF(Time To First Fix), 
C/NO, Navigation Accuracy. 

¢ Test Methodology and Test Conditions: 


— ION STD 101 : Recommended Test 
Procedures For GPS Receivers 


TTFF Test 


¢ Calculation of 
Minimum / ________Warm Start TTFF Results(Example 
Average / Sa Te eee 
Maximum / 
Standard 
Deviation 


Even Start Time 
among 1 GPS 
Frame 


a 


eA fae 


Ephemeris in Navigation Message 


eo Ten, 30-bit words, 6-s subframe 


1 TLM | HOW | SV Clock correction + Satellite Health 


2 HOW Ephemeris 
1-frame 
3 | TLM | HOW | Ephemeris continued | 30s 
1500 bits 
4 Almanac+lonospheric UTC correction 
5 


TLM | How | Almanac | | 


Pamsuned MAMIE PAA 47 


ELECTRONICS. 


Static Navigation Accuracy 


¢ Calculation of 
95%, 99.99%, 
68%, 50% 
percentile / 
Minimum, 
Maximum / 


Percent Valid 
¢ Position Plot 


a 


ELECTRONICS. 
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Environmental Test 


¢ Temperature/Humidity 

— RF components, especially "= 
RF crystal, are severely 
effected by changing 
Temperature 

— Monitor Position Accuracy 
and C/NO with respect to 
given Temperature Profile. 


a 


ELECTROMICS 


Environmental Test 


* Shock & Vibration 


— RF crystal is subjected to Vibration and 
Shock. 


— Monitor C/NO with respect to given 
Vibration Profile. 


— Monitor the mounting of each components. 


ae 1 a 


¢ AMPS/CDMA/PCS Interference Test 
— Layout, Shielding 

¢ EMC(Electro-Magnetic Compatibility) Test 
— Radiated Emissions Test 
— Specially important for Automotive Application 
— Layout, Shielding 

* ESD Test 


am” 


ELECTRONICS: 


412 GPS =ADD| JDlISo set 


© EAS/HAS GPS FAD 
— Motorola, SiRF, ST Electronics 


¢ GPS/DR Integration 


¢« Network Aided GPS 
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EAS/IN M2 GPS AAT 


¢ Motorola Oncore(Features) 


— GPS RF Front End and 12 channel correlator in 
ONE chip. 

— Dual Conversion Architecture with Fractional 
Synthesizer onboard : Crystal reference range 
from 10 MHz to 26 MHz 


— Low Power Consumption : 60mW for the 
combined RF and digital sections. 


— Small Size: 64 Pin MAP BGA(8mm X 8 mm) 
— Full production : Q1 2001 


«am 


ELECTRONICS 


Motorola Oncore 


¢ Block Diagram(Preliminary) 


Data Bus GPS Antenna 


Microcomputer Motorola Oncore 


Seria! 0 


Crystal 
Reference 
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AABN MS GPS FAI 


L_ 


¢ SiRFstar Il 


— Enhanced GPS Core based on high 
performance SiRFstar | : 1920 
time/frequency search channels 

— Integrated ARM7TDMI 

— On-chip 1Mb EDO DRAM for GPS navigation 

— Real-Time Clock and battery backed SRAM 

— TricklePower Mode for power saving 


PS AMsunced Atal MT PAARL? 


SiRF starll 


¢ Block Diagram 


=302> 


EAST MH GPS AAI 


¢ ST Microelectronics 


— GPS RF Front End(STB5600) and CPU with 12 
channel correlator(ST20-GP6). 


— very low cost oscillator (no TCXO) 
— On-chip 1Mb SRAM for GPS navigation 


— Real-Time Clock and Peripherals 


ST Microelectronics 


* Block Diagram(ST GiPSy) 


_ 
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One Chipst : RF 2 Digital°! One 


“SRAM &! CPUS! One Chip3} 
cA oO}: Kt A & OFA AG AAS 2A Al 


7| is Lt Power Saving Mode(Low Power 


Mode) 7 &. 
—Embeded GPS solution for Mobile 
Application (RF crystal & host processor | 


SF) 
& TT 


a 


GPS/DR Integration 


¢ DR(Dead Reckoning) 
— Current position is calculated from a knowledge of 
an initial position and measurements of speed and 


direction of the vehicle.(uses gyros 
accelerometers, odometer, wheel tick sensors) 


x 
North 
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GPS/DR Integration 


¢ GPS vs DR 


2| +4 Sto] 8} 
-Blockage 
-Multipath 


GPS/DR Integration 


* Switched/Filtered GPS-DR System 


— The navigation solution is switched/combined 
between/from GPS and DR. 


= When GPS solution quality is good, GPS velocity 
is used to calibrate bias and scale factor errors. 


GPS Receiver 
Position, Velocity Switch/ 
DR System ee Kalman Filter Filtered 


Position, Velocity 
DR Sensor Error Estimates 
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GPS/DR Integration 


¢ Integrated GPS-DR System 


— Uses DR sensors to propagate(time update) GPS 
receiver's navigation state. 


— Uses measurement update to compute navigation 
solution as well as DR sensor error estimates. 


— Provides more accurate state estimate and 
improves reacquisition and multipath rejection. 


GPS Receiver | Position, Velocity 


DR Sensor 


Conexant GPS/DR System 
¢ Block Diagram(Jupiter/DR) 


Digitat Processor(Scorpio) RF MCM 


GPS Antenna 
lo IF Down 
hecvetee| 


Crystal 
Oscillator 
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Network Aided GPS 


FCC E911 Phase II Rule(Background) 
— Automatic Location Identification(ALl) 


capability requirement for handsets sold in 
Unite States. 


— Position Accuracy (Handset based 
solution): 50 m(67%), 150 m (95%) 

— Deployment Schedule: start selling ALI- 
capable handsets beginning from October 


Network Aided GPS 


Restrictions for Conventional GPS Receiver 
— Does not work indoors ! 


= Need to increase sensitivity 
— Too Slow !! 


=> Need to get Navigation Messages somewhere else 
— Too Big and Expensive !!! 


— Need to incorporate into existing handset hardware 


a 
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Sensitivity Enhancement 


— CDMA Network time is synchronized with 
GPS time 

— Coherent Integration time can be increased 
up to 20 ms for CDMA/PCS system. 


pp ee cas 


|! Second =1.023.000 Chips! 
=1000 PN Epoch 9 b= =~ 
=50 0, 


fa Bits | 


AAMT PAA 


Navigation Messages 


— Navigation Messages can be received via 
network. 


— Received Navigation Messages can be 
used to prolong coherent integration time 
for sensitivity improvement. 

= Position Determination Service Standard 
for Dual Mode Spread Spectrum Systems 
(TIA/EIA/IS-801 ) 


a 
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TIA/EIA/IS-801 


« -Mobile Terminated Response ¢ Mobile Originated Response 


Message(BS => MS) Messages(MS => BS) 
— Provide BS Capabilities - Provide MS Information 
— Provide GPS Acquisition Assistance ~ Provide Autonomous Measurement 
- Provide GPS Location Assistance Weighting Factors 
- Provide GPS Sensitivity Assistance — Provide Pseudorange Measurement 
— Provide Base Station Almanac — Provide Pilot Phase Measurement 
— Provide GPS Almanac - Provide Location Response 


Provide Time Offset Measurement 


- Provide GPS Ephemeris 
- Provide GPS Navigation Message Bits 
~ Provide Location Response 

Provide GPS Almanac Correction 


a 


_— > % tie 
GPS satettites ~~ % 


a 
| 4-—— GPS signais —_—\, 


Reference GPS 
receiverfservice 
j 


Mobile w/ Position 


Partial GPS £8 OSASHOL ats bart > Determining 
N <£ = : 


\ 


assisting 
information 
about CPS : 
signal S Mobile £ 
Switching | 


Base statinns Center — Mobije 


Positioning i 
Center 
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Qualcomm gpsOne 


¢ Network Aided GPS + Advance Forward Link 
Trilateration(AFLT) = Hybrid Solution 
¢ Signal coverage of wireless network and GPS 
satellite complement each other. 
— Rural & Suburban area : sparse base station coverage, 
more than four GPS satellites 
— Urban & In-Building : few GPS satellites and more 
base station. 
* gpsOne support: MSM3300(IS95), 
MSM5100(CDMA2000 1X) 


¢ gpsOne in 
MSM3300 
—MSM:Mobile 
Station 
Modem 
— CDMA 


Processing 
ASIC + CPU 


+ Peripherals 1 a 
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¢ 1S801 
Message 
Flow _— 
(Example) : 


equest GPS Acquisition Assistance 
equest GPS Sensitivity Assistance 


Bal vs 


: Provide Pseudo-range Measurement ; | (c) 
i : Provide Pilot Phase Measurement | : 


| Request Location Response 
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DGPS #2 2t 


Phone (02) 450-3358 
FAX (02) 444-6670 
Email : younglee@kkucc.konkuk.ac.kr 


a2 _| 
a 


Sw hs SS eed 
National Research Lab. 


20004 128 121 
Ml 7 % GNSS AAG Fees 


@ GPS HAA 

@GPS S21} 2&7 asc 
@ DGPS CA PA Be 

@® RTCM SC-104 #2 7:3 
ez=> AAA DASA 
ee eS OS Wes 
® DGPS #2 

@ DGPS S2 =0} 

@ DGPS #d A2Vze 


=r Sa Atstm Gps aAl2e awa 


National Reasearch Lab 


=o les 


GPS #2 Oat Olol 


eseanrc 
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6 IA] SRL FR OR 


eASAS OF : ASHE ASAT AS OH 

eASNA SH : ASHE ASAI SH OA 

ets Ct :GPSlso As SUA ASA 2? AA Oa 
eckea ot :GPsdiso WA SAN ASAl2t AA On 
eC3ae ox :Gpsiiso OS Sel Bet On 


@esAI CH 3: SES, SHUAY SA, Hae 2aor, 


= 
ay ek 


SN Peps irae er Research Lab. = CH St GPS Al AG Aaa 


HOI A SAS ASLazt 


QF ds 


HOLAA & it RH OF} 


SPs Pa ee APD GPs A\AG APs 


Nationsal Reasearch tab 


—314- 


GPS #8 @ *#(1) 


Filtered UERE 


FAL QRt (VDOP = 2.5) 
BA Lxt (HDOP = 2.0) 


ABs GPs AlAw APH 


GPS BAL A(2) 


rms SAE AP] 10(m) 
Random 
Aga & 
AAA 
mals 
aa 
Clea =e 
FAl7| 


UERE 
Filtered UERE 


tinea eee oe ASWetm GPs AJA ATA 


=a lo 


Fan os 


GPS 2el 2 at (3) 


rms 2X 37] 10(m) 


SA Ala 
Hels 

We 

sae 

+Al7| 

UERE 

Filtered UERE 

SA] QA (VDOP = 2.5) 
BA Lat (HDOP = 2.0) 


Nationa Rotosrcn tae Ahem Gps Az Apa 
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Al+2! Gla 


&¢]3 (Ionosphere) 


43H (Stratosphere) 


— 7 bes 3 a 


National Research Lab. 


Geomagnetic fj } 
Equator” 


== reas ee 


pana © 
National Research Lat. 


2317 = 


150 ~ 1,000km 


Asst GPs A|A8 A> 


Al 
= 


acls2 sd 


@ AA 5OkmO/ Al 1,000kmAOIOA) 2cl ERA RE AN 
j2sel aad 
@ 100MHz O/dt 2) ATE StAtotLt, VHFO| AS] SAH 1OOMHz)= Si} 


(GPS #F ST} : L1 = 1.575GHz, L2 = 1.228GHz) 

© GPS A's SaAQ| SHO AASAS FR Ct Lol 

e ADAS Al2t AAS BS AAt+* (TEC: Total Electron Content)(l 
Heo) 42D 2] ASO Bul ell 


MEA 


mtictat test oorer toe A2Wan GPs AA#! apAé 


e A= 2 DA 
ee ESB 
@ EHet Si} HA (114 471): OS ADA 2000 ~ 20014 
@ X= Aela AS (Geodetic) S2/A& SOE NT 
XtI| (Geomagnetic) 3&/A&0l AAl 


, — Ps a 
AFA National Research Law 


es eo a 


AclSOl Get Al2t AS SJ (2) 


@ 100MHz O/ eS] ANU Aclss 


0 Polarization Rotation 


0 Angular Reflection 
0 Amplitude and Phase Scintillation 
0 Group Delay and Carrier Phase Advance) 

¢ DC (Group Delay): Sit 327 S ANSOl Udlote= 
AS Al2zt AS 
SI} (Carrier Phase Advance) : 2 
Ai2zt tS aa 


ee 


Le i 
National Research Lab. ae e 


® Dual Frequency #4!)| : LID L2VSTS| SA Aj2t HOS 
OSat AANA Ba 


A, = 2rec (2), A,, =P rec () 


fi fi, 
@ Liv} L2 Sit GSOCeS AGAA20/ Cec. 
> L129] I}4AI} SOOS We] SHA <ALD 
A =A = Ay = 03 rec - > TEC 
Sis Si 
_ 40.3 1 _ 40.3 


mee 
—)TEC TEC 
c eT Te € ee ‘fia 


arrival 


TEC a Reve al” Gere meer fy fe ) 
40.3 fi. -Srs 


Sr es 


Nationat 
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© Differential GPS : 7) 221} ASA AOS] SSLMES AD 
@ Klobuchar 72! 0/S 
-AiPASQ WS CS 


SHHAAS Sol C/A DS MASAO AS 
- 2 50% 459 A sit 


SPQ SSeS ee ASUSD GPs AA ap! 


Research Lab. 


Ads SA BS (3) 


Tae = DC+ Avos =) : Zenith Delay 
DC= BFTUS5,10), 42098, 7-H ALZH(50,4002%), p= Fl 


Fe=l+ 2{ 2S.—2) 
90° 


2a(t-T.) 
P 


: slant factor 
Y= 
A : : Daa 
(= is +GPStime : Local Time Ionospheric point 
A=Sa@-, P=SfPO. ,a,R:CPS BAHN AS Says 


®. = ©4116 cos(a,— 291): AL PAPIA S 


4 


b nee S \ x 
‘ome =F %4 ALL -— 4 — \ _ 
T. +|sx10 + ( : +z) for RI<5 


= Fx [5107] for ble 5 


Pe Sea ae ASWSD GPs AA APY 
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Klobuchar Model 


Time Delay (nono sec) 


10 15 
Lecal Time (Hourt 


ee es AShsa Gps A|Agq aaa 


National Resaarchn Lab. 


Ae 2D 


ee : 1995. 3. 21. 


Sea WSUS GPs A|AB] otra! 


National Research Lab. 


=el 


@ A> Galas as & 
Lcd 2ECh. 


e SHWE YAor BC. 


OQ 
Cc 


A e|S (Ionosphere) 


150 ~ 1,000km 


, hie — Ss Oa os 
7 National Research Lat. 


a 


esa Jan Gpstislt ASAHS L2] SLOA 
HOLEAL 2 ALE BAY (ALG & At) 
eet, 25,55 59 SHO Oe} AM Soo HS 
e2CO BST 2E aI AAS es 
@® 15GHZ)}7A| SHS WAU, FHA Fe 

® GPS 4! S0/l Attenuation, Scintillation, Delay S°] Ae 


ZH Se AZOWeD GPs Ag ape 
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1 OF (2) 


i=] 


4 (dry constituents)1} GA 4 =(wet constituents)°O & 
OWE, +2] SHS Dae. 
712+(elevation angle) 0i| [é 


: mapping function 
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Se AA QA (3) 


TO 


A=[A,,m,(E)+A,,m,(E)] | 


7], A,, =dry zenith delay 
A, —=Wwet zenith delay 


m, (E) = dry mapping function (obliquity factor) 
m, (E) = wet mapping function (obliquity factor) 


Note : m,(E) =m, (E)& D2t7} WEBI SEO] BA 


asst Gps Al|A8 Ara 


: me I- . — ha 
as National Rosearch Lab. 


observation ‘SS 


rd ee | = eS 
National Research Lak. 


earth's 
surface 


A=ahst Gps A|\Ae ata 
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CSA (Multipath) 


eAdgs Sst AT VAS SBS Sad Ad aes 
(direct and indirect pathh= *AlI/0] [Soka Ba 
eAA)| #20 Get We OLA UELUSS SSQnS 
AAW AIF] Il BIS 
@ Als 47/7 SAO] Hs 
0 Al2t AS 
1a= 
ES 
OAA Ass 


¢ — > HBuWstD GPs AA ATA 


National Research Lab 
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CSAS OT) S43 (1) 


@ SAO CSoeE ASA AZAES lst offset2] YYOS 
Al SA QAI Bs 

e Arisa On & Ollie! AS LockingO| Ot 3+ US 
e200] Oe SA0i 40] (AAS OAS Fills S4) 


— areas eS 


pom 5 
National Research Lap. i—] = 


CSAS ONY Sa (2) 


e 2CO Ds CHE SPSE6~ 102 FI 

© C/AZEO Ast Oz 
OLA SI! OAS) OF AAI] LOA PE IHS 
0 Ch QUE 1.5chip ASS Ot :1466ns (2¢ 440 m) 
02£1.2~1.4m aS 


@PDCOl Sst OA 
0 S/CHO Tt 1.Schip S&F SO/SS C/ABESE S& B12 
HAT PAO QS! SA 
0 2ICH 1/4 AOIB2S] Si Acm ASE (L129) BF SF 4em)2 Lit 
122 Aad BS ABOU SS lem BE OWA OF 
0 S2é LSA ore 


Asst GPs AlAe ata 
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CSAS Sst JH (1) 


@ +Al)| CEILI BADE AS 
OORILE FA SAS DATA 22) 
OAls Sas Ole 
: right-handed circular polarized signal 
: left-handed circular polarized signal 
°} Gain PatternS 22010 D0] KES WSS SA 


0 52! (choke ring)S NSH 42 toe se2ote As ate 


Ose SA2 SS 


fee Se Se ASS GPs AA aya 


Be aaa AHS GPS AIAG oA 
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CSAS Sst I (3) 


@ #A0/ AZEMOO! Sst oH 


0 Correlation Function’s Leading Edge 

0 Narrow Correlator Technology (‘90 - ‘93) 

0 Correlation-Function Shapes (‘94 - ‘95) 

0 MEDLL (Multipath Estimating Delay-Lock Loop) 
0 Strobe Correlator (‘96) : MEDLL®] Faz ls 

0 Modified Correlator Reference Waveforms 
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ASAZTSAA 


e GPSA Ao ACA AZo] ASE NSA AA LAO AS 
ae 
@eAS DB Al2t Bele AS 
0 53412] 23 (Monitor Station)lA AGS 
OSAMA HES SUA FAHAdA Ald 
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SHA AS Sol ASAI AS 
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20, ASS As 
0 Dual Frequency? 4!9|, MAIS AIA, IIACIOIE SS 
0 30A\2! GOIKS SA (GPSALI SAA WD 6Al2! El 
ee SHOHAAS ACASE WSS ROSS, 4 aa aIOagS 
Ont HS 
@ARSAAN ALOHE AAD FAIS ALAM El (line-of-sight 
vector)0] =Se S| HS OA YA 
eAANAS GSO rms LOS 3.6m ASS CAI} BH 
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@ GPSS 0/22 AASAA 
@ GPSF1 22! AIA 
0 Block I, ITA : AIA}S (29H) + SHICIS (27H) 
OBlock UR : AAS (19H) + SHICIS (2N)Ce 
4H WOU, SUS AAS See a 
@ GPS OCS(Operational Control System)= GPS Time= 
UTCOl 1 msecO|J WIS ASCs FX 
e2 AMO AALS 1 msecS CS YEOS SHAS ONS 
AS: ANNAS SEAARSS ASS MH 


SS Se ASHSD GPS AAG CFA! 
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Selective Availability (SA) JH 2 (1) 


@ 197085 GPS BAAA C/AZEY) Oz 
20~30mafe2 HSU, OSHA] DAS ANC! Hal 
SaQ7 24 

@ Aa} 500m (2drms)% CU, 100mz 2A (198344) 

S29/9] ALAA! VORG} blwet AD 
Random Processing® = 42: WS0| AAS 
AIK| SE 60~80m A&P Oz 
eS SHISAS 
~ AAAS Ae 24 :¢-BenAdDAINSos 
- AM AA A (clock dither) : 5-Z2&AIA > FLAN 


Bre Swe ae ABest GPS A|Ag azal 
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. FAA GPS Product Team, http://gps.faa.gov/ 
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v Network-assisted GPS 
v Principle 
v Case studies 
o Snaptrack 
o SIRF 
o Qualcomm 


Intelligent Vehicle Lab. Wireless Location Technologies 2 


= 301= 


Introduction 


v Personal navigationOll Het AA 2 2p Sa 
v Wide deployment of Wireless network 
o Cellular, PCS, TRS, IMT2000 


o Wireless data communication infra including 
internet access 


v Need for Personal navigation 
o Safety 
— Wireless 911, Anti-kidnapping, etc. 
o Location-based commerce 
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Introduction 


v E-911(Enhanced 911) in U.S. 
Wireless 911 calls represent about 43% of all 911 calls 
More than 25% could not identify their location 


FCC mandate to locate emergency mobile telephone calis 
(CC Docket No. 94-102, Sep/1999, www.fcc.gov) 


Two phase implementation to be completed 
o Phase |: location of cell site/ or sector(by Apr/1998) 
o Phase Il : location of users with an accuracy of the 
followings(by Oct/2001) 
— Network-based solution : 100m for 67%, 300m for 95% 
of calls 


— Handset-based solution : 50m for 67%, 150m for 95% 
of calls 
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! Introduction | | | 


v Other applications of wireless location 
v Location-based commerce 
o Location-based billing 
o Location-based yellow page 
v Traveler information 
o Personal navigation aids 
o Vehicle navigation & route guidance 
v Fleet management 
v Tracking 
o Asset tracking for freight carriers 
o Personal tracking, child tracking 
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Introduction | 


v Wireless location®o] 22a 
v GPS 
o Good positioning accuracy 
o Limited in open sky(outdoor) 
o Long TTFF(time-to-first-fix) 
o Needs an additional equipment in MS(mobile station) 
v Wireless Location 
o Positioning using existing cellular network 
o No additional equipment in MS 
— Efficient solution for personal navigation 
o Indoor positioning capability 
o Accuracy limitation 
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| Wireless Location Technology | 


v Radiolocation 
v Location derived from simple geometric relationships 
v 2 types of location system 
o Reverse link location : remote positioning 
— BSs received signal from MS 
— No additional equipment in MS 
— Complexity added to network 
— Security & Privacy problem 
o Forward link location : self positioning 
— BSs received signa! from MS 
— Minor revision in MS 
— No additional equipment in network 
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| Wireless Location Technology | 


v Radiolocation(cont’d) 
v 3 types of measurements 
o Signal strength 
o AOA(angle of arrival) 
o Time-based 
— TOA(time of arrival) 
— TDOA(time difference of arrival) 
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Wireless Location Technology 


v Signal strength 


v Empirical propagation models describe signal 
attenuation with distance 


Several models for macro-, micro-, and pico-cell 
Signal attenuation : proportional to d? 
Noise sources : shadow fading, multipath fading 


Very inaccurate in multi-path fading with empirical 
models 
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| Wireless Location Technology | 


v AOA(direction finding) 
v Need antenna arrays or directional antennas 
o Calculation of direction 
a=sin'!(d/l) 


Array antenna 
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Wireless Location Technology 


v AOA(cont’d) 
v Location at intersection of straight lines of position 
v Measurement 


Line of position 
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! Wireless Location Technology | 


v AOA(cont’d) 
v Solution 
o Usually 2D solution 
o Needs more than 2BSs for 2D solution 
v Environment dependent : scatter effects 
o Macrocell 
~ Ring of scatters : good accuracy 
o Microcell 
— Elliptical and dual scatter models : bad accuracy 
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Wireless Location Technology 


v Time-based methods 
v Time synchronization of network is required 
v TOA - 
o Line-of position 

— Circle 

o Location 
— Intersection 
Of LOP 


Line of position 
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Wireless Location Technology 


v TOA 

v Measurement 

o Travel time t, 
Cty PAK) VS F? 

v Solution of TOA 
o MS Clock may have bias 
o Needs more than 4 BSs for 3D solution 
o Linearized LS(least square) 
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| Wireless Location Technology | 


v Time-based methods(cont’d) 
v TDOA 
o Differencing between BS ___ Line of position 


o Line of position ae 
~— Hyperbola “ 
o Location i 


a. i ‘ Soe oat re 
- Intersection of LOP A Gee A 


’ 
o 


: CN tae 
rts} 
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v TDOA 

v Measurement 

o Travel time t, 

Ct gig PAL XY HO VDP? AG XY AO Y DFE 

v Solution of TOA 

o MS Clock bias is eliminated 

o Needs more than 4 BSs for 3D solution 

— More than 3 measurement equations 
o Linearized weighted LS 
o Analytic solution for 3 or 4 BSs 
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Wireless Location Technology 


v Time-based methods(cont’d) 


v Linearized LS solutions of TOA and TDOA are 
equivalent 


o Same position estimate 
o Same error covariance 
v Error sources 

o Multipath fading 

o NLOS(non-line-of-sight) propagation 
— Typically 400-700m 

o MAI(multiple-access interference) 
— Near-far effect(by power control) 
— Particularly in CDMA system 
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| Wireless Location Technology | 


v Limitation in wireless location 
v Hearability 
o For wireless location, 3(2D) or 4(3D) BSs should be 
related 
o Hard to hear from 3 or 4 BSs in large portion of 
cellular coverage areas 
v Accuracy 
o By multipath, MAI, and NLOS, positioning error is 
much larger than GPS 


Intelligent Vehicle Lab. Wireless Location Technologies 19 


Wireless Location Technology 


v Alternatives 
v Network Assisted GPS 
o GPS positioning assisted by wireless network 
o Assistance from network : Timing, Doppler, MS clock 
correction, visible SVs PRN No., ephemeris(SV 
position), PR correction, etc. 
o Fast TTFF, Indoor capability 
v Hybrid method 
o AOA+TOA/TDOA 
o TOA/TDOA+GPS(including network assisted) 
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Wireless Location Technology 


v Location Service Categories 


Category 1 
Basic Service Level 


Location of all handsets, including legacy, 
based on cell or improved cell accuracy 


Location of all new handsets with improved 
accuracy at a reasonable cost 


Category 2 
Enhanced Service Level 


Location of new handsets with high accuracy 
and higher (compared to Category 2) cost 
with customer choice 


Category 3 
Extended Service Level 


I I YT SNP EEA ES TTT TE TE SS EE IE OS FF EN 
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Wireless Location Technology 


v Categories Mapping 


System | Category | Location Technology 


4 Improved Ci location method. The improved Cl location method utilizes in addition to 
Cl information e.g. TA and other GSM specific measurement data. This accuracy 
category can be implemented without changes to MS. 


GSM specific E-OTD (MS assisted/based) positioning technology, which is, 
standardised for GSM Release'98. 


The MS Based GPS concept which is standardised for GSM Release'98. 
Cell coverage based method, which is already included in ANSI-136. 
Currently no support in TIA/EIA-136 C standard exists for this category. 


MS Based GPS location concept that is standardised for TIA/EIA-136 C. 
Recognising the commonalties and convergence in TIA/EIA-136 and GSM, the 
TIA/EIA-136 implementation of NW Assisted MS Based GPS location technology is 
harmonised with the GSM GPS concept to the greatest extent practical. 


Ci based method, which is standardised for TIATR45.2 IS-41. 


IS-95 specific AFLT (MS assisted/based) position technology, which is standardised 
for TIA/IS-801. 


MS Based GPS concept which is standardised for TIA/IS-801. 


—_ 
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v Related Companies 


Company Technology 
TruePosition TDOA 
Cell-Loc 
Grayson wireless TDOA 
Lockheed Sanders TDOA 
SnapTrack Cellular Aided GPS 
US Wireless Multipath fingerprint 
SigmaOne AOA 
Cambridge Positioning Systems Cross Correlated TDOA 


Harris Corporation MicroTrax 
Corsair — Phonetrack AOA/TDOA 
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Network Assisted GPS 


v Main concept 
v Limitation in GPS positioning 
o Long TTFF in cold or warm start 
o Positioning only when open sky(outdoor) is provided 
v Assistance from wireless network 
o Time assistance 
— In CDMA cellular, GPS time is used as reference 
o Frequency assistance 


— MS clock drift can be compensated by cellular carrier 
tracking 


o Doppler/position assistance 


— SV Doppler is provided from a location server that is 
attached to cellular network 


— SVID, ephemeris, PR correction from location server 
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Network Assisted GPS 


vy Main concept(cont’d) 
v Solution 
o Fast TTFF : Using assistance data(time, ephemeris, 
frequency, Doppler) 
o Sensitivity 


— Using frequency assistance, MS clock drift can be 
compensated 


— By multiple convolution(multiple PN sequence set), 
sensitivity is increased 


PL ae ee 
WLU ULL Sore 


— By GPS navigation data, # of overlapping PN sequence 
is limited 
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MOBILE OR REMOTE UNIT 
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Frequency 
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Tranernit 
Shifts, and Other Signal Parameters: 
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Network Assisted GPS 


v SnapTrack solution(cont’d) 


Receive Fix Command Power down and 
& Ooppler from Base wait for new command 


Find oscittator offset Find position 
from Modem AFC ot Send Data to Base 


Turn On Receiver 


Callect snapshot 
of GPS signa! 


Turn Of Receiver 


Turn On OSP 
Geurn Sede foe Interpolate Peak Location 
the Satellite 
Ooppler correction ware. Piocese 7c 


Relative timing determination 
using FFT Based Matched 
Filtering 


mn A ED 
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| Network Assisted GPS | 


v SnapTrack solution(cont’d) 


Se etter ee Find Maunliideeauared of Navigation data 
ae demodulation 
Multiply by FFT of PN code Inverse FFT 
Compensate for 
Doppler Time shift 
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v SnapTrack solution(cont’d) 

v Fast TTFF 

o Network assisted Doppler 

o GPS local oscillator calibration using BS carrier 
v Sensitivity 

o Multiple convolution : use 10ms block 

o GPS local oscillator calibration using BS carrier 
v Time assistance 

o Time from BS : Time tag 


o Correlation results to PDE(position determination 
entity) with time tag 


v Low Power Consumption 


Intelligent Vehicle Lab. Wireless Location Technologies 30 


=o 


v SnapTrack solution(cont’d) 
v PDE(position determination entity) 

o Attached to BSC or MSC 

o Find exact position with 
— Ephemeris 
— Correlation result from MS 
— Time index for measurements from MS 
— PR correction if necessary 
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Network Assisted GPS 


v SIRF 


v Wireless Aided GPS(WAG) technology + SIRF’s GPS 
receiver technology 


v Fast TTFF 
o Network assisted time and Doppler information 
o Uncertainty is reduced to about 10us and 150Hz 
o Drastically reduces the size of the receiver’s search 
— Network assisted time and Doppler information 
— Uncertainty is reduced to about 10us and 150Hz 
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| Network Assisted GPS | 


v_ SiRF (cont'd) 
v High sensitivity 
o Provides satellite position information or ephemeris 
o Cell phone clock sharing with GPS 
o Reduce the integration time per dwell 
— If navigation bit transition information available 
> Integration time up to 20ms 
~ If bit phase information available 
> Integration time beyond 20ms 
v Low power consumption 
o SiRF’s TricklePower techniques 
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Network Assisted GPS | 


v_ SiRF (cont'd) 
v Robustness to assistance information inaccuracy 


o SiRFstar architecture feature a large number of 
satellite-searching correlators. 


v Supports various kinds of location mode 


o Network centric, network assisted, and autonomous 
mode 


o Can reduce Wireless Geolocation Server computing 
load 


o Flexible enough to accommodate a “broadcast type” 
of location assistance information 
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| Network Assisted GPS 


v SiRF (cont'd) 
v “Network Centric” mode : remote positioning 
o Location computation is made at the server 


o Acquisition assistance can be high level(almanacs) or low 
level(code phase and carrier frequency) 


o GPS in MS deliver only some variation of pseudorange 
v “Network Assisted” mode : self positioning 

o Location computation is made at MS 

o Stand-alone or DGPS position is sent back to the server 
v “Autonomous” mode : no assistance 


o Location is computed locally at MS, and is sent back to the 
server, with or without DGPS correction 


o Little or no assistance information 
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Network Assisted GPS 


v Qualcomm’s solution(MSM 3300) 
v “gpsOne solution” 
o Similar to Snaptrack’s solution 
o Time domain correlation 
v Implementation 
o Antenna switch 
o GPS correlation using cellular components 
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| Network Assisted GPS | 


v Standards 
v TIA/EIA 
o 1801 : CDMA cellular signaling standard 
o PN3890 : CDMA cellular network entities standard 
v ETSI 
o GSM spec. 
o Being moved to 3GPP 
v 3GPP(3" Generation partnership project) 
o 3G TS 22.071, 3G TS 23.171, 3G TS 25.305 
v T1P1 
v ITU 
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v E911 Network Reference Model(IS-801) 


E 


MSC 
A D, Emergency Services 


Emergency 
Services 
Network 

Entity 


Emergency 
Services | - 
Message 


Entity Legend 
CRDB: Coordinate Routing Database 
MPC: Mobile Position Center 

MSC: Mobile Switching Center 

PDE: Position Determining Entity 
PSAP: Public Safety Answering Point 
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